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I.     INTRODUCTION 

A.  DIRECT  FORM  STRUCTURE  REALIZATION 

A  scheme  for  high  speed  recursive  digital  filter  realization  was  suggested  by  H.  H. 
Loomis  and  B.  Sinha  [Ref.  1].  The  main  aspect  of  the  design  is  in  modifying  the  filter 
difference  equation  so  that  the  recursive  part  requires  only  the  Lth  past  output  instead 
of  the  immediate  past  output.  This  removes  the  tight  recursive  requirement  thus  enabling 
pipelining  of  multipliers  and  adders  in  the  direct  form  implementation  resulting  in  high 
throughput.  One  of  the  difficulties  with  this  method  is  that  the  modification  introduces 
extra  poles  (not  necessarily  stable)  and  zeros  in  the  filter  transfer  function  which  cancel 
each  other  in  the  infinite  arithmetic  implementation.  In  practice  due  to  the  finite  word 
length  effect  this  cancellation  may  not  occur.  This  may  result  in  an  unstable  realization 
even  if  the  original  filter  were  a  stable  one.  However,  if  L  is  chosen  to  be  large  enough, 
then  it  has  been  proven  that  the  new  poles  introduced  will  also  be  within  the  unit  circle 
thus  guaranteeing  stability  even  if  the  pole-zero  cancellation  does  not  occur.  This 
method  is  called  the  direct  form  structure,  or  the  scalar  structure. 

B.  STATE  SPACE  REALIZATION 

Another  method  by  R.  Gnanasekaran  is  an  alternate  technique  which  relaxes  the 
recursive  nature  of  the  state  equation  (state  space  realization)  [Ref.  2],  That  is,  the 
present  state  vector  computation  does  not  require  the  immediate  past  state  vector.  The 
proposed  technique  allows  high  speed  implementation  using  the  state-space  structure 
directly.  High  throughput  is  achieved  by  use  of  pipelined  multiply  -add  modules.  The 
actual  hardware  complexity  will  depend  upon  the  method  of  pipelining,  though  very  high 
speed  is  possible  by  increasing  the  number  of  states  in  the  pipelining.  This  method  im- 
plements AL  instead  of  A.  the  system  matrix.  The  statespace  realization  will  be  stable 
if  the  original  filter  is  stable.  Furthermore,  the  A1  filter  has  improved  tolerance  to  coef- 
ficient error. 

C.  OVERVIEW  OF  THESIS 

Chapter  II  describes  the  IIR  notch  filter.  The  transfer  function  and  the  coefficients 
that  are  appropriate  for  the  IIR  notch  filter  are  developed. 

Chapter  III  presents  the  realization  of  the  high  speed  IIR  notch  filter  by  using  the 
pipeline  technique  of  Loomis  and  Sinha.  Stability  testing,  frequency  response  and  unit 


sample  response  are  described.  System  designs  for  pipelined  multiply-add  modules  are 
included. 

Chapter  IV  presents  Gnanasekaran's  state  space  realization  which  includes  a  modi- 
fied state  space  equation.  System  designs  for  pipelined  multiply-add  modules  are  de- 
scribed. 

Chapter  V  contains  a  comparison  of  the  results  of  the  Loomis/Sinha  pipeline  re- 
cursive design  and  the  Gnanasekaran  state  space  realization  design. 

Appendix  A  contains  the  program  and  data  for  frequency  response  analysis  of  the 
II R  notch  filters  in  Chapter  III. 

Appendix  B  contains  the  program  and  data  for  unit  sample  response  analysis  of  the 
II R  notch  filters  in  Chapter  III. 


II.     IIR  NOTCH  FILTER  DESIGN 

From  Roue,  [Ref.  3:  p.  72]  the  IIR  notch  filter  has  the  transfer  function 

b0z   -\-b]zJrbc) 

H{z)  =        2        ' °-.  (2.1) 

OqZ   +  bxz  +  a2 

Dividing  the  numerator  and  denominator  of  (2.1)  by  a„,  we  have 


bo 
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by 

bo 
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+  -^r 

17.-) 

4] 

ao 

1 

>- 

a2 
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Cl{ 

~  z 
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a0 

H{z)  =  - — —  .  (2.2) 


#p  bx  a-, 

Let  —  =  BQ.  —  =  B,.  and  —  =  A0.   Therefore  the  transfer  function  is 

#0z2  +  fl,z  +  B0 

U{z)  =     '-' ,  *-.  (2.3) 

From  (2.1)  and  (2.3)  we  can  see  that 

1.  B{  in  the  numerator,  must  be  equal  to  Bl  in  the  denominator;  and. 

2.  the  first  B0  in  the  numerator  must  be  equal  to  the  last  50  in  the  numerator. 


Rewriting  (2.3)  as 


Z?0(z2  +  ^-z+l) 


H{z)  =  — u- .  (2.4) 

(z2  +  B,z  +  A0) 

the  transfer  function  H(z)  is  of  the  form 

{z-dl){z-d2) 

Let  the  zeros  be  placed  at  le~8,  and  let  the  poles  be  placed  at  JA0  e-e.    Therefore 
the  transfer  function  is 


B0(z   -\e-jB)(z-\e+jB) 
z-^e~J°){z-^e+JU) 


H^  =  -       — -;^      .-r^  (2"6) 


^0[^-lg+;^-lg'yey2+l] 

\_zl  -  Jlje+J9z  -  jA~e~J9z  +  JA0  ] 


^)°.,2     "V^.   ~^z     —  (2-7) 


Z?0[z2  -2(  ;  )z  +  1] 

^) _      +/     -;* (2-8) 

[z2  -  2V^0  (  '       2  )Z  +  A°'] 

B{)[z2  -2  cos  0 z  +1] 

//(z)  =      ,  0L     _ i—  (2.9) 

[z   -2jA0  cos  Oz  +  Aq] 

B{)z2  -  2B0  cos  6z  +  B0 

//(z)  =  -t = -■  (2.10) 

z   -  2VM0  costfz  +  ^0 

According  to  (2.3)  we  can  see  that 
-2B0  cos  0z  =  -2jl^  cos  #z  (2.11) 

which  means 
B0  =  JA^.  (2.12) 

For  our  example  of  the  1 1 R  notch  filter  we  want  the  zeros  to  be  placed  at 

!  e±/0.44671 10339*  =1<?±#  ^ 

and  we  want  the  poles  to  be  placed  at 

O.99^°-446711O33^  =  v/^e±/0.  (2.14) 
We  can  see  that 

0  =  0.446711033971  (2.15) 

A0  =  0.9801  (2.16) 

50  =  7^  =  0.99.  (2.17) 


This  realization  is  shown  in  Figure  1  on  page  6. 
From  equation  (2.1    ) 

B^z2  -  IB,  cos  Bz  +  Bn 

^)=  2    :  — — fl — 7-  (2'18j 

z    -  2y/A0  cos  0z  +  ^40 

Substituting  6,  B0  and  /_0.  we  then  have 

0.99z2  -  2(0.99)  cos(0.44671 10339rr)z  +  0.99 

H(z)  =  —7 y- - - '- (2.19) 

z2  -  2(0.99)  cos(0.4467110339:r)z  +  0.9801 

0.99z2- 0.32993; +  0.99 

_--0.32993r-0.9S01 

,        0.99-0.32993:-'40.99:-;  ,  . 

1  -  0.32993z_1  +  0.9801z~2 

Because  we  have  the  transfer  function  [Ref.  1:  p.  273] 

dn  +  a,z     4-  flhz 

1  —  b-z      —  b2z 

from  (2.21)  and  (2.22)  we  can  see  that  all  of  coefficients  are 

Gq  =  0.99 

a,  =  -0.32993 

a2  =  0.99  (2.23) 

b}  =0.32993 

b2=  -0.9  S01  . 

*(;-)  =  ||  (2.24, 


F(z)        0.99  -  0.32993z-1  +  0.99 
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(2.25) 


A'uj        j  _  o.32993z_1  +  0.9S01z~2 

F(z)[l  -  0.32993z_1  +  0.9S01z~2]  =  A'(z)[0.99  -  0.32993z_1  +  0.99z"2]  (2.26) 

Yiz)  -  0.32993z_I  Y(z)  +  0.9S01z~2Fiz)  =  0.99X(z)  -  0. 3299  3  z~lX\z)  +  0.99z~2X(z)    (2.2") 
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Figure   1. 


Poles  and  Zeros  Location 


Y{z)  =  0.99A'(z)  -  0.32993z-IA'(z)  +  0.99z~2A'(z)  +  0.32993z_1  Y[z) 
-0.980l2"2}'(z) 


(2.28) 


y{n)  =  0.99.r(«)  -  0.32993.v(w  -  1)  +  0.99^r(«  -  2)  +  0.32993y(«  -  1) 
-  0.980 \y{n-  2). 


(2.29) 


This  is  the  basic  second-order  recursive  filter  (p  =  0).  The  number  of  p  will  be  ex- 
plained later  in  Chapter  III  and  the  realization  is  shown  in  Figure  2  on  page  7. 

We  can  design  this  IIR  notch  filter  to  be  a  high  speed  IIR  notch  filter  by  using 
pipeline  techniques  so  that  it  can  be  operated  at  hithero  impossibly  high  rates.  In 
Chapter  III  we  will  apply  Loomis'Sinha  method  and  in  Chapter  IV  Gnanasekaran's 
state  space  realization.  A  comparison  of  the  results  of  these  two  techniques  will  be 
shown  in  Chapter  V. 


x(n) 


o 


On 


y(n) 


■-cr-a 


Define  Symbols 


x(n) 

y(n) 


Unit  Delay 
Pipeline     Multiplier 

Pipeline  Add  Unit 

Input 
Output 


Figure  2.       Second  Order  II R  Notch  Filter:  p  =  0 


III.     REALIZATION  OF  IIR  NOTCH  FILTER  BY  USING  THE 
LOOMIS/SINHA  METHOD 

A.     PIPELINE  RECURSIVE  DESIGN 

From  Chapter  II  we  have  the  second  order  IIR  notch  filter  that  is  described  by  the 
transfer  function 

_,        (an  +  «iz     +  a->z     )  , 

H(z  ])  =  -^ l— 2—±  (3.1) 

(1-V     ~b2z  2) 

and  the  difference  equation 

y{n)  =  b]y{n  —  1)  +  b-jy{n  —  2)  4-  aQx(n)  4  axx{n  —  1)  +  a2x(n  —  2).  (3.2) 

The  computation  of  r(«)  requires  the  immediate  past  output  y{n  —  1)  .  This  filter  can 
be  developed  into  high  speed  IIR  notch  filter  by  using  the  pipeline  technique  of  Loomis 
and  Sinha  [Ref  1].  The  main  aspect  of  the  design  is  in  modifying  this  filter  difference 
equation  so  that  the  recursive  part  requires  only  from  the  Lth  past  output  instead  of  the 
immediate  past  output.  This  means  the  speed  limits  can  be  overcome  by  modifying  this 
difference  equation  so  that  y{n)  depends  not  on  the  immediate  past  output.  y(n  —  1).  but 
on  y{n  —  L)  for  some  L  >  1.  This  can  be  accomplished  as  follows: 

For  L  =  1  the  difference  equation  is 

y(n)  =  b^in  -  1)  +  b^y{n  -  2)  +  %x{n)  +  +  axx{n  -1)4-  a2x{n  -  2).  (3.3) 

We  can  see  that  this  equation  is  a  2-order  difference  equation. 

From  [Ref.  1:  p.  274]  the  difference  equation  has  {p  4  m)r/,-order  where  m  is  the  order 
of  the  difference  equation  and  p  is  the  pipeline  delay  used  to  overcome  the  speed  limit. 
Therefore  equation  (3.3)  has  the  number  of  p  =  0.  We  can  call  this  equation  the  second 
order  IIR  notch  filter  with  p  =  0. 

Delaying  equation  (3.3)  by  1  yields 

y(n  -  1)  =  btf{n  -  2)  4  b^in  -  3)  4  a0x{n  -  1)  4  axx{n  -  2)  4-  a2x{n  -  3).  (3.4) 

For  L  =  2.  putting  (3.4)  into  (3.3)  results  in  the  difference  equation 


y{n)  =  b^byyin  -  2)  4-  b-#{n  -  3)  4-  a0x{n  -1)4-  a^n  -  2)  4-  a2x(n  -3)] 

+  b,j\n  -  2)  +  a0jc(«)  +  a^(«  -1)4-  a2^(«  -  2)  (3,5) 

>■(«)  =  b2y(n  -  2)  +  />; ^v(n  -  3)  4-  ^.x^  -1)4-  axbxx{n  -  2)  +  c^x^  -  3) 
+  bxv(n  -  2)  +  a0jr(«)  4-  <*,*(«  -  1)  4-  a2;c(/7  -  2) 

>•(«)  =  (^2  +  62)y(«  -  2)  +  &,M«  -  3)  +  a0x{n)  4-  (a06,  4-  ax)x{n  -  1) 

4-  (a,^  4-  a2)x{n  -  2)  +  a2V("  -  3).  l  '  J 

We  see  that  the  computation  of  >'(«)  requires  only  the  second  past  output  y(n  —  2) 
instead  of  immediate  past  output  y{n  —  1),  and  we  can  use  this  technique  for  the 
higher-order  of  L  . 

We  see  that  equation  (3.7)  is  a  3-order  difference  equation  and  has  p=  1.  Thus  we 
can  call  this  filter  the  augmented  IIR  notch  filter  for  p=  1. 

The  z  transform  of  equation  (3.7)  is 

Y  =  (b2  4-  b2)z~2  Y  4-  b,  b2z~2  Y  4-  a0X  +  {aQbx  +  a^z^X 
4-  («]/?]  4-  a2)z~  X  4-  a2bxz~  X 

y  ^  a0  4-  (a0b]  4-  <7,)r~'  +  (q^  4-  o2)z~2  +  a2bxz~2 
X~  l-(b2l+b2)z-2-blb2z-2 

_,         (ars  +  a.z-   +  a->z~  )(l  +  b,z~  ) 
(l-b.z   l-b2z  2)(1+V   ') 

That  means  the  transfer  function  of  the  augmented  IIR  notch  filter,  p=  1.  is  equal 
to  the  transfer  function  of  the  second  order  IIR  notch  filter,  p=  0,  multiplied  and  divided 
by  the  same  polynomials. 

We  can  design  the  higher-order  difference  equation  by  starting  from  a  second-order 
difference  equation.  That  means  we  can  design  the  higher  order,  p-augmented  IIR  notch 
filter  by  starting  from  our  second  order  IIR  notch  filter,  p=0,  as  follows: 

From  Chapter  II  the  second  order  IIR  notch  filter  transfer  function  is  the  original 
transfer  function 

_i         (Or,  4-  a,z      4-  a-yz     )  ,„ 

H(z  l  =  -^ ]— 2-^-.  (3.11) 


Multiply  the  numerator  and  denominator  of  the  original  transfer  function  by  the 
same  polynomial  equation.  This  can  be  thought  of  as  equation  (3.12)  As  developed  in 
[Ref.  1:  p.  2S1].  the  transfer  function  of  the  p-augmented  filter  is  as  shown  in  equation 
(3.13)  We  know  from  other  considerations  that  the  transfer  function  of  the  higher  order 
realization  must  in  fact  be  the  same  as  the  original  and  therefore  must  have  the  same 
poles  and  zeroes  in  the  z  -  plane. 

The  coefficients  of  the  a  polynomial  depend  on  the  original  bx  and  b2,  and  p  is  de- 
termined by  the  number  of  stages  in  the  pipeline  multiply  and  add  units. 


(a0  +  a,z~x  +  a2z-2)(l  +  afV1  +  -  +  ofz~p) 


H'(z~l)  = 


(1  _  blZ~l  -  b2z~2)i\  +  a.fz-]  +  -  +  4\~p) 


(3.12) 


//'(0  = 


a 


(p)„-i 

o  z 


+  -  +  ap+2z 


(p)     -p-2 


,<P)  ,-/>-! 

>p+V 


l-b™*~r-1  -b™n 


(p)     -p-2 

p+2z 


(3.13) 


The  values  of  bf\  o.f\  and  a'f  are  given  in  Tables  1,  2  and  3  [Ref.  1:  pp.  278-279]. 


Table  1.      &'sfor/?  =  0.1..  ..,5 


Coefficient  of  D  ~l 

Coefficient  of  Z>-2 

p 

tyli 

b'^\ 

0 

b. 

b2 

1 

b\  +  b2 

bxb2 

2 

b\  +  2bxb2 

b\b2  +  b\ 

bi  +  lb\b2  +  b\ 

b\b2  +  7bM 

4 

b\  +  Ab\b2  +  lb& 

b\b2  +  lb\b\  +  b\ 

5 

b\  +  >b\b2  +  6b\b\  +  b\ 

b\b2  +  Ab\bl  +  Ibfil 

10 


Table 

2.       a. 

's  for  /> 

=  0.1 5 

P 

af> 

a<f> 

a£5 

y.r 

a.f 

0 

0 

0 

0 

0 

0 

1 

ft, 

0 

0 

0 

0 

-> 

ft. 

&?  +  Z>2 

0 

0 

0 

3 

ft> 

ftf  +  ft2 

b\  +  2bxb2 

0 

0 

4 

^ 

ft?  +  ft2 

b\  +  2b}b2 

b\  +  lb\b2  +  b\ 

0 

5 

ft, 

ft?  +  ft2 

b\  +  2bxb2 

b\  +  lb\b2  +  b\ 

ft?  +  Ab]b2 

+  Zbfil 

Table  ?.       a  's  for  p  =  0.1 3 


p 

^ 

af> 

a'p 

a'f 

<& 

c#> 

0 

flo 

«i 

a: 

0 

0 

0 

1 

«o 

t70&,  -j-  0., 

a]b]  +  a2 

a2b} 

0 

0 

2 

«0 

aub[  -r  a. 

a,(b\  +  b2) 
+  (TTj^,  +  a2 

«,(ft?  +  ft2) 
+  a26, 

a2{b\  +  b2) 

0 

_> 

<*0 

allbl  +  «, 

a0(b:  +  b2) 
+  a-b.  +  a2 

a0(^  +  2ft]ft2) 

+  a-(bl  +  b2) 
+  a2b, 

+  ai(b\  +  2fc,fc2) 
+  a,(ft?  +  ft2) 

+  a2(6?  +  2ftLft2) 

B.     SECOND  ORDER  IIR  NOTCH  FILTER:  P  =  0 

From  equation  (3.13)  in  section  A,  we  let  p  =  0.  This  equation  becomes: 


_,         (aQ  +  a,z      +  a-,z     ) 


-l 


(l-ft1Z-'-ft22    ') 


(3.14) 


H(2     }-A-(z) 


(3.15) 


7(z)        (a0  +  a,z   1  +  fl2z  2) 


A'(~)        (1  -  V"1  "  ft2z"2) 


(3.16) 


Manipulating  (3.16)  will  result  in  the  difference  equation,  (3.20). 


-l 


}  U)[l  -  ftiz     -  ft22"  ]  -  A'(r)[flb  +  fl!2      +  fljz  '] 


(3.17) 
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Y(z)  -  V~~'  Y{z)  -  b2z~2Y(z)  =  a<>X{z)  +  a]Z-lX(z)  +  a2z-2A'(z)  (3.18) 

Y(z)  =  a0X{z)  +  axz-\X{z)  +  a2z~2X(z)  +  b,z~]  Y(z)  +  b2z~2Y(z)  (3.19) 

y(n)  =  a0x(n)  +  axx(n  -  1)  +  a2x(n  -  2)  +  b${n  -  1)  +  btf{n  -  2).  (3.20) 

This  realization  is  shown  in  Figure  3  on  page  13. 

From  chapter  II  we  know  that  for  our  example  notch  filter,  that  the  coefficient 
values  are 


a0  =  0.99 

a,  =  -0.32993 

fl2  =  0.99  (3.21) 

bx  =  0.32993 

b2  =  -0.9801  . 

The  transfer  function  is 
H'(z~1)  =  ^ l— ~—r-.  (3.22) 

(l- v-    -  V~") 

The  difference  equation  is 
>•(«)  =  a0x{n)  +  axx(n  -  1)  +  a2x(n  -  2)  +  bxy(n  -  1)  +  btf{n  -  2).  (3.23) 

C.     AUGMENTED  IIR  NOTCH  FILTER:  P  =   I 

From  equations  (3.12)  and  (3.13)  in  section  A,  we  let  p  =   1.  Those  equations  be- 
come: 

_,         (an  +  tf,z      +  a-)Z        ItR    z     ) 
H'(z   [  =f2 i— 1—* 1—1  (3.24) 

(1-V   '-62z  ^(l  +  a^z   ') 

^      )  -  /     ,(i)-2        WD-3  (3"25> 

1  —  £>Vz      —  p,  z 


-i        Hz) 


12 


Figure  3.       Second  Order  IIR  Notch  Filter:  p  =  0 


vv   \  0)    ,       0)   -1    ,       0)   -2    ,       (1)   -3 

}  (z)        aQ'  +  a]  'z     +  a2  z      +  a]  z 


X(z) 


I  _  *<V2  -  6<V3 


Manipulating  (3.27)  will  result  in  the  difference  equation,  (3.31). 


(1)   -2         .<!)   -3 


(D    .    „<>>,-!    ,    Jl)   -2    ,     J!),-3- 


r(2)[i  -  b^z-  -  b^z-^  =  xwlw  +  fli"2-  +  <^*    +  *^z-J] 
r(r)  -  *<  V2  roo  -  »1V3  }'(z)  =  4»xz)  +  fl{,>2-,jr(2)  +  a<Vlr(z) 

+  a?)z-3X(z) 
Y(z)  =  aJ1)A'(z)  +  aJV^z)  +  afV"2^*)  +  ajV^z)  +  AJV2^) 

+  *<V3r(z) 

>•(«)  =  c&x)x{n)  +  a\])x(n  -  1)  +  fl?}jc(«  -  2)  +  «?>*(«  -  3)  +  b^y{n  -  2) 

+  tf  !k»  -  3). 

This  realization  is  shown  in  Figure  4  on  page  14. 

From  Tables  1,  2  and  3  [Ref.  1:  pp.  278-279],  we  have 


(3.27) 


(3.28) 
(3.29) 

(3.30) 

(3.31) 
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Figure  4.       The  Augmented  IIR  Notch  Filter:  p  =   1 


a0    - 

1 

«?>- 

*i 

=  0.32993 

*0)- 

ao 

=  0.99 

,0) 


a}"  =  00^,  +  ^  =  -0.0032993 


.(i) 


«2    =  «i*i  +  a2  =  0.SS11462 


(3.32) 


,n) 


a],;  =  a2^  =0.3266307 
b™  =  bl  +  b2  =  -0.871246 
^,>  =  /,1^,2  = -0.32336439 


The  transfer  function  is 


l_6(')2-2_6(D2-3 


(3.33) 


The  difference  equation  is 
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>•(«)  =  a{J\x(n)  +  a\l)x(n  -  1)  +  a^x{n  -  2)  4-  a(3l)x{n  -  3)  +  b^yin  -  2) 


D.     AUGMENTED  IIR  NOTCH  FILTER:  P  =  2 

From  equations  (3.12)  and  (3.13)  in  section  A,  we  let  p  =  2.  Those  equations  be- 
come: 

U  (1  -  ^z"]  -  V"2)d  +  afVJ  +  42)z~2)  l       j 

,        <ff  >  +  fljV1  +  aPz-2  +  *<V3  +  afz~" 

H'(z    )  =  ;      -^Z~3 — 7777^ (3.36) 

1  -  b\  'z     -  b\  'z 

H'(z  X)  =  ^  (3.37) 


Y(z)        a{2)  +  aJV1  +  fz-2  +  a<V3  +  a{2)z~4 


W)  1  -  b?z-3  -  &<V 


Y(z)  -  bf )z"3  Y(z)  -  b^z^Yiz)  =  42)X(z)  +  a{2)z~]X(z)  +  afz~2X{z) 

,       (2)   -3  v,   v    .       (2)   -4  v,   v 

+  t?3  ;z     A(z)  +  a}  ;z     A(z) 


Y(z)  =  42)X(z)  +  a{2)z~lX(z)  +  apz-2X(z)  +  a(22)z~2X{z)  +  af  z"4Z(z) 
+  bi22)z-2Y(z)  +  b{2)z-4Y(z) 

y(n)  =  ^xin)  +  a(2)x{n  -  1)  +  of ]x{n  -  2)  +  af^(«  -  3)  +  a(^x(n  -  4) 
+  6fy(«  -  3)  +  b{2)y(n  -  4). 

This  realization  is  shown  in  Figure  5  on  page  16. 

From  Tables  1,  2  and  3  [Ref.  1:  pp.278-279],  we  have 


(3.38) 


Manipulating  (3. 38)  will  result  in  the  difference  equation.  (3.42). 

7(z)[l  -  b?z-3  -  bPz-*l  =  X(z)l42)  +  a{2)z~l  +  afz~2  +  a{2)z~3  +  afz"4]  (3.39) 


(3.40) 


(3.41) 


(3.42) 
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Figure  5.       The  Augmented  IIR  Notch  Filter:  p  =  2 
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(3.43) 


a{2)  -  1 

a!f]  =  bx  =0.32993 

af  =  b\  +  b2  =  -0Sl  12462 

oJ2)  =  flo  =  0.99 

«{2)  =  a^  +  a,  =  -0.0032993 

«22)  =  ao(*f  +  b2)  +  alb]  +  a2  =  0.01861246 

af}  =  a,(62  +  b2)  +  a2bx  =  0.61408096 

a£}  =  a2{b]  +  b2)  =  -0.86253373 

£<2>  =  &f  +  26^2  =  -0.61081465 

^2)  =  £262  +  £22  =  O.S5390S4. 

The  transfer  function  is 

"(r  )= i-*;v>-»«v •     '  (    ' 

The  difference  equation  is 

v(w)  =  oJ2)jc(k)  +  flP)JC(»  -  1)  +  a^xfa  -  2)  +  af^w  -  3) 
+  a?]x{n  -  4)  +  ^2)y(«  -  3)  +  bf\in  -  4). 


(3.45) 


E.     AUGMENTED  IIR  NOTCH  FILTER:  P  =  3 

From  equations  (3.12)  and  (3.13)  in  section  A,  we  let  p  =  3.  Those  equations  be- 
come: 

(        i  -1    ,  -2\/i     ,       (3)   -l    ,       (3)   -2    ,       (3)   -3x 

^2      '  n         i.     -I        L     -2\/i     ,       (3)   -1    ,       (3)   -2    .      (3)   -3,  1 .3.40) 

"'<z  )-- — ! --w^-u^r " —  <3'47> 

1  —  b\  z      —  b\  z 

A\Z) 


vv   \  (3)    ,       O)   -1    ,       (3)   -2    ,       (3)   -3    .       (3)    -4    ,       (3)   -5 

Y(z)        %    +  a\  z      +  a2z      +  a]  z      +  a\  z      +  a\  'z 


XW  l-^-V4-^-5 


(3.49) 
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Manipulating  (3.49)  will  result  in  the  difference  equation,  (3.53). 

r(z)[i  -  tf  >z~4  -  ^V5]  = 

X(z)U2)  +  ajV1  +  afz~2  +  «f  >2"3  +  afz"4  +  a<V5] 


1  -  tf }^"4  "  ^"5 


The  difTerence  equation  is 


(3.50) 


Y(z)  -  b^z -*Y{z)  -  bf}z~5Y(z)  =  aj3)*(z)  +  a^z    X(z)  +  43)z~2X(z) 

+  afz~3A'(z)  +  a^z^Xiz)  (3.51) 


(3.52) 


(3.53) 


+  «<3)  z~5X{z) 

Y(z)  =  43)X{z)  +  a\3)z-]X(z)  +  afz-2X{z)  +  af  }2_3X(z)  +  a^z~^X(z) 
+  a{pz-5X{z)  +  bPz-*Y(z)  +  b{3)z-5Y(z) 

y(n)  =  43)x{n)  +  af]x{n  -  1)  +  a^x{n  -  2)  +  af\{n  -  3)  +  a{^x{n  -  4) 
+  af^n  -  5)  +  bi3)y{n  -  4)  +  bf}y(n  -  5). 

This  realization  is  shown  in  Figure  6  on  page  19. 

From  Tables  1,  2  and  3  [Ref.  1:  pp.  278-279],  we  have 

43)=1 

af }  =  Z»,  =0.32993 

43)  =  £)2  +  62  =  -0.S7 12462 

43)  =  6?  +  2£,£2  =  -0.61081465 

a<3)  =  fl()  =  0.99 

«j3)  =  a0bl  +a1  =  -0.0032993 

a{3)  =  o0{b;  +  b2)  +  tfjfcj  +  a2  =  0.01861246 

a<3>  =  «0(/.3  +  2^2)  +  flj^i  +  62)  +  a2bx  =  0.00937445 

fl<3)  =  Ql(/>3  +  2fej62)  +  a2(£?  +  b2)  =  -0.66100766 

«s3>  =  ai(b\  +  MM  =  -0.6047065 

^3)  =  b*  +  3b2b2  +  b\  =  0.65238232 

b{3)  =  b3b2  +  2bJ>\  =  0.59865944  . 

The  transfer  function  is 

(3)    ,       (3)    -1     ,       (3)   -2    ,       (3)    -3    ,       (3)    -4    ,       (3)    -5 

.        g^;  +  a\  >z      +  a^  'z     +  a\'z      +  a\  >z      +  a\  'z 


(3.54) 
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Figure  6.       The  Augmented  II R  Notch  Filter:  p  =  3 
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y(n)  =  43)x(n)  +  af^rt  -  1)  +  of  }x(n  -  2) 

+  a? \x(n  -  3)  +  a(^x(n  -  4) 

n  (3.56) 

+  af}*(«  -  5)  +  ^3V(«  -  4) 

+  b^y(n-5). 

F.     STABILITY  TESTING 

From  section  A  if  we  start  with  a  stable  filter,  we  would  like  to  be  assured  that  the 
poles  of  the  augmented-order  filter  all  be  outside  the  unit  circle  in  the  z~x  plane.  We  de- 
sired this  because,  even  though  the  added  poles  are  cancelled  by  the  zeroes  of 
(1  +  ocwz_1  +  •••  +  c/.^z-"),  realization  imperfections  will  prevent  exact  cancellation  and  the 
augmented  filter  would  be  unstable.  In  order  to  examine  the  stability  question,  we  will 
use  Jury's  stability  test  as  applied  to  successively  higher-degree  denominator  polynomials 
in  z  as  p  increases  [Ref.  1:  p.  281].   The  results  from  this  paper  are  summarized  below: 

1.     Condition  for  Stability  of  the  Original  IIR  Notch  Filter:  p  =  0 

The  original  IIR  filter  is  described  by  equation  (3.14)  and  to  test  for  stability, 
we  analyze  the  following  characteristic  polynomial  resulting  from  it: 


-> 


F{z)  =  ;"  -b]Z-b2  =  0  (3.57) 

The  conditions  for  stability  are: 

\-bl-b2>0...O 

l  +  bl-b2>0...Q>  (3.58) 

\b2\  <1...© 

These  conditions  are  shown  graphically  in  Figure  7  on  page  21. 

Circled  numbers  in  the  figure  refer  to  the  conditions  of  the  same  number  in  equation 

(3.58) 

Checking  for  stability  of  our  example  IIR  notch  filter;  we  have  from  ®  of  (3.58) 

1  -  bx  -  b2  >  0 

1- (0.32993) -( -0.9801)  >0 
1.65017  >0 

from  ©  of  (3.58) 
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Figure  7.       Condition  for  Stability  of  the  Original  Filter  [Ref.  1:  p.  283] 

1  +  b]  -  b2  >  0 

1  +  (0.32993) -( -0.9801)  >0 

2.31003  >0 

from  ©  of  (3.58) 


1^1  <1 

0.9S01I  <  1 

0.9S01  <  1 


Thus  the  conditions  are  all  satisfied  and  the  original  filter  is  stable. 

2.     Condition  for  Stability  of  the  Augmented  IIR  Notch  Filter:  p=  1 

The  augmented  IIR  notch  filter,  p=  1,  is  described  by  equation  (3.33)  and  to  test 
for  stability,  we  analyze  the  following  characteristic  polynomial  resulting  from  it: 


F(z)  =  2  +  ^=0 


(3.59) 


The  conditions  for  stability  are: 
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1  +  *,  >  0  ...  ® 

^  (3.60) 

These  conditions  are  shown  graphically  in  Figure  8  on  page  23. 
Circled  numbers  in  the  figure  refer  to  the  conditions  of  the  same  number  in  equation 
(3.60) 

Checking  for  stability  of  our  example  IIR  notch  filter;  we  have  from  ©  of  (3. 60) 

1  +  b)  >  0 

1  +  0.32993  >0 

1.32993  >0 

from©  of (3.60) 

I  -  bx  >  0 

0.67007  >  0 

Thus  the  conditions  are  all  satisfied  and  the  augmented  IIR  notch  filter,  p  =  1,  is  stable. 
3.     Condition  for  Stability  of  the  Augmented  IIR  Notch  Filter:  p  =  2 

The  augmented  IIR  notch  filter.  p=  2,  is  described  by  equation  (3.44)  and  to  test 
for  stability,  we  analyze  the  following  characteristic  polynomial  resulting  from  it: 

F(z)  =  z2  +  b1z  +  (b2l+b2)  (3.61) 

The  conditions  for  stability  are: 

1  +  bx  +  (b2  +  b2)  >  0  ...  ® 

1  -  bx  +  {b2  +  b2)  >  0  ...  ©  (3.62) 

\b[  +  b2\  <1  ...  © 

These  conditions  are  shown  graphically  in  Figure  9  on  page  24. 
Circled  numbers  in  the  figure  refer  to  the  conditions  of  the  same  number  in  equation 
(3.62) 

Checking  for  stability  of  our  example  IIR  notch  filter;  we  have  from  ©  of  (3. 62) 

1  +  bx  +  {b2  4-  b2)  >  0 

1  +  (0.32993)  +  [(0.32993)2  +  (  -0.9801)]  >  0 

0.4586838  >0 


22 


\A 

Lb2 

1 

Z 

■ 

i 

k          -*Di 

-i        l 

• 

Figure  8.       Condition  for  Stability  of  the  Augmented  Filter:  p  =   1  [Ref.  1:  p.  284] 
from©  of (3.62) 

1  -  £>,  +  {b\  +  b2)  >  0 

1  -  (0.32993)  +  [(0.32993)2  +  (  -0.9801)]  <  0 

-0.201 1762  <0 

therefore  the  filter  is  unstable,  from  ©  of  (3.62) 

\b2  +  b2\  <1 

I  (0.32993)2  -f-  (  — 0.9S01)  I  <  1 

0.8712462  <  1 

Note  that  ©  is  violated  and  here  the  realization  will  be  unstable. 

4.     Condition  for  Stability  of  the  Augmented  IIR  Notch  Filter:  p  =  3 

The  augmented  IIR  notch  filter,  p=  3,  is  described  by  equation  (3.55)  and  to  test 
for  stability,  we  analyze  the  following  characteristic  polynomial  resulting  from  it: 


F(z)  =  z3  +  b,z2  +  (b2  +  b2)z  +  [b]  +  2bxb2) 


(3.63) 
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Figure  9.       Condition  for  Stability  of  the  Augmented  Filter:  p  =  2  [Ref.  1:  p.  285] 


The  conditions  for  stability  are: 


1  +  6j  +  (62  +  fc2)  +  {b*  +  2^)  >  0  ...  0 
1  -  fcj  +  (£»f  +  b2)  -  {bl  +  2bxb2)  >  0  ...  © 
\b\  +  lblb1\  <  1  ...  © 
I  {b\  +  Ib.bJ2  -  1 1  >  I  b,{b\  +  2£,62)  -  (£2  +  *j)  I 


(3.64) 


© 


These  conditions  are  shown  graphically  in  Figure  10  on  page  25. 
Circled  numbers  in  the  figure  refer  to  the  conditions  of  the  same  number  in  equation 
(3.64) 

Checking  for  stability  of  our  example  IIR  notch  filter;  we  have  from  ©  of  (3. 64) 

1  +  by  +  [b]  +  b2)  +  (bl  +  2bxb2)  >  0 

1  +  0.32993  +  [(0.32993)2  +  (  -0.9801)]  +  [(0.32993)3  4-  2(0.32993)(  -0.9801)]  >  0 

-0.15213085  <0 
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Figure   10.        Condition  for  Stability  of  the  Augmented  Filter:  p  =  3  [Ref.  1:  p.  2S6] 
therefore  the  filter  is  unstable,  from  ©  of  (3.64) 

1  -  6,  +  [b]  +  b2)  -  {b\  +  2bxb2)  >  0 

1  -  0.32993  +  [(0.32993)2  +  (  -0.9801)]  -  [(0.32993)3  +  2(0.32993)(  -0.9801)]  >  0 

0.40963S46  >  0 

from©  of (3.64) 

\b3x+2bxb2\  <1 

I  (0.32993)3  +  2(0.32993)(  -0.9S01)  I  <  1 

0.61081465  <  1 

from  ©of  (3. 64) 

Kif  +  26162)2-l|  >  \bi{b3l+2blb2)-(bl  +  b2)\ 

I  {(0.32993)3  +  2(0.32993)(  -0.9801)}2  -  1 1  < 

I  (0.32993){(0.32993)3  +  2(0.32993)(  -0.9801))  -  {(0.32993)2  +  (  -0.9801)}  I 

0.62690546  <  0.66972012 
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therefore  the  filter  is  unstable. 

Note  that  ©  and  ©  are  violated  and  therefore  the  realization  will  be  unstable. 

G.     FREQUENCY  RESPONSE 

From  the  transfer  functions  (3.22),  (3.33),  (3.44),  and  (3.55),  we  use  the  program  in 
Appendix  A  to  calculate  the  magnitude  and  phase  of  the  frequency  response  of  the  se- 
veral realizations.  We  plot  the  magnitude  of  these  responses  for  the  original  filter  and  for 
the  pipeline  augmented  that  all  of  the  frequency  responses  are  the  same  from  those  fig- 
ures. 

1.  Second  Order  IIR  Notch  Filter:  p  =  0 

The  magnitude  response  is  shown  in  Figure  1 1  on  page  27. 

2.  The  Augmented  IIR  Notch  Filter:  p  =   1 

The  magnitude  response  is  shown  in  Figure  12  on  page  28. 

3.  The  Augmented  IIR  Notch  Filter:  p  =  2 

The  magnitude  response  is  shown  in  Figure  13  on  page  29. 

4.  The  Augmented  IIR  Notch  Filter:  p  =  3 

The  magnitude  response  is  shown  in  Figure  14  on  page  30. 

We  know  from  section  F  that  the  example  augmented  IIR  notch  filter  for  p  = 
2,  and  3  are  unstable  but  we  cannot  see  that  from  these  magnitude  responses.  We  have 
to  find  another  way  to  show  which  one  is  stable  or  unstable.  This  will  be  shown  in  sec- 
tion H. 

H.     UNIT  SAMPLE  RESPONSE 

From  the  difference  equations  (3.23),  (3.34),  (3.45),  and  (3.56),  we  use  the  program 
in  Appendix  B  to  calculate  the  output  for  a  unit  sample  input  to  each  of  the  pipeline 
realizations.  We  plot  the  outputs  in  Figures  15,  16,  17  and  18.  We  can  see  that  the  out- 
puts for  Figure  15  (p  =  0)  and  Figure  16  (p=  1)  are  stable,  however  the  outputs  for  Fig- 
ure 17  (p=2)  and  Figure  18  (p=3)  are  unstable  as  predicted  by  the  stability  testing  in 
section  F. 

1.  Second  Order  IIR  Notch  Filter:  p  =  0 

The  output  is  shown  in  Figure  15  on  page  31. 

2.  The  Augmented  IIR  Notch  Filter:  p  =   1 

The  output  is  shown  in  Figure  16  on  page  32. 

3.  The  Augmented  IIR  Notch  Filter:  p  =  2 

The  output  is  shown  in  Figure  17  on  page  33. 
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Figure  11.       Magnitude  Response  of  Second  Order  IIR  Notch  Filter:  p  =  0 
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Figure  12.       Magnitude  Response  of  the  Augmented  IIR  Notch  Filter:  p  =  1 
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Figure  13.       Magnitude  Response  of  the  Augmented  IIR  Notch  Filter:  p  =  2 
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Figure  14.       Magnitude  Response  of  the  Augmented  I1R  Notch  Filter:  p  =  3 
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Figure  15.       Output  Sequence  of  Second  Order  IIR  Notch  Filter:  p  =  0 
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Figure  16.       Output  Sequence  of  The  Augmented  II R  Notch  Filter:  p  =  1 
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Figure  17.       Output  Sequence  of  The  Augmented  IIR  Notch  Filter:  p  =  2 
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4.     The  Augmented  IIR  Notch  Filter:  p  =  3 

The  output  is  shown  in  Figure  18  on  page  35. 

I.     PIPELINE  IMPLEMENTATION 

All  of  the  IIR  filters  in  sections  B,  C,  D,  and  E  can  be  developed  into  high-speed  IIR 
filters  by  applying  pipeline  techniques.  We  use  Ka-stage  pipeline  add  units  and  Km-stage 
pipeline  multipliers  to  increase  the  speed  of  the  filters.  The  Km  stages  are  involved  in  the 
coefficient  multiplies  and  the  Ka  stages  are  involved  in  the  adds.  The  increase  in  order 
of  the  filter  is  P  =  2Ka  +  Km-2  [Ref  1:  p.  280]. 

1.     Second  Order  IIR  Notch  Filter:  p  =  0 

To  achieve  a  p  =  0  realization,  we  cannot  use  pipeline  modules  with  longer  de- 
lays than  Ka  =  1  and  Km  =  0. 

P  =  2Ka  +  Km-2 

=  2(1) +  0-2 
=  0 

We  have  a  realization  of  the  original  IIR  notch  filter  but  with  an  additional  3 
delays.  The  equation  is 


-3 


Y=z~3[aQX+a]z~]X+a2z~2X  +  b]z~lY+  b2z~2Y'].  (3.65) 


This  realization  is  shown  in  Figure  19  on  page  36. 
2.     The  Augmented  IIR  Notch  Filter:  p  =   I 

To  achieve  a  p=  1  realization,  we  cannot  use  pipeline  modules  with  longer  de- 
lays than  Ka  =  1  and  Km  =  1 . 

P  =  2Ka  +  Km  -  2 
=  2(1)+  1-2 
=  1 

We  have  a  realization  of  the  augmented  IIR  notch  filter,  p  =   1,  but  with  an 
additional  4  delays.  The  equation  is 

z~AY  =  z~\a0X  +  alz~lX  +  a2z~2X+  a3z~3X  +  b2z~2Y  +  b3z~3Y'].  (3.66) 

This  realization  is  shown  in  Figure  20  on  page  37. 
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Figure  18.       Output  Sequence  of  The  Augmented  IIR  Notch  Filter:  p  =  3 
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Figure   19.       Second  Order  II R  Notch  Filter:  p  =  0  with  Ka=  1,  Kra=0 

3.  The  Augmented  IIR  Notch  Filter:  p  =  2 

To  achieve  a  p=2  realization,  we  cannot  use  pipeline  modules  with  longer  de- 
lays than  Ka  =  1  and  Km  =  2. 

P=2Ka+  Km -2 
=  2(l)  +  2-2 

_  t 

—  .*- 

We  have  a  realization  of  the  augmented  IIR  notch  filter,  p  =  2,  but  with  an 
additional  5  delays.  The  equation  is 

2"5  Y  =  z~\oqX  +  a^-]X  +  a2z~2X  +  a3z~3X  +  a^X  +  b3z~3  Y  +  b^  Y].  (3.67) 

This  realization  is  shown  in  Figure  21  on  page  38. 

4.  The  Augmented  IIR  Notch  Filter:  p  =  3 

To  achieve  a  p=  3  realization,  we  cannot  use  pipeline  modules  with  longer  de- 
lays than  Ka  =  1  and  Km  =  3. 
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Figure  20.       The  Augmented  IIR  Notch  Filter:  p  =  1  with  Ka=  1,  Km=  1 
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Figure  21.       The  Augmented  IIR  Notch  Filter:  p  =  2  with  Ka=  1,  Km  =  2 
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P=2Ka+  Km-2 

=  2(1)  + 3 -2 
=  j 


We  have  a  realization  of  the  augmented  IIR  notch  filter,  p 
additional  6  delays.  The  equation  is 


-6V  -6 

z     r  =  z 


OqA'+^z   'X+a2z  2Z+a3z  3X 

+  a.^X  +  a5z~5X  +  64z"4  Y  +  b5z~5  Y 


3,  but  with  an 


(3.6S) 


This  realization  is  shown  in  Figure  22  on  page  40. 


Table  4. 

LOOMIS/SINHA  PIPELINE  RECURSIVE  DESIGN 

p-order 

Ka-add  units 

Km-multipliers 

D  (output  delayed) 

Figure 

0 

1 

0 

3 

19 

1 

1 

1 

4 

20 

2 

1 

2 

5 

21 

3 

1 

6 

22 

J.     SUMMARY 

This  chapter  investigated  a  method  for  applying  pipeline  techniques  to  the  design 
of  a  high-speed  IIR  no-  ii  filter.  Using  the  Loomis  Sinha  method,  IIR  notch  filters  op- 
erating at  rates  hitherto  impossible  can  be  designed.  The  general  structure  of  the  IIR 
notch  filter  and  the  ^ithod  of  calculating  the  multiplier  coefficients  have  been  presented. 
The  stability  of  the  resulting  realization  has  been  investigated  and  a  technique  for  as- 
certaining the  stability  of  the  realization  has  been  presented.  The  delay  and  complexity 
of  each  realization  are  summarized  in  Table  4. 
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Figure  22.       The  Augmented  II R  Notch  Filter:  p  =  3  with  Ka=  I,  Km=  3 
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IV.     GNANASEKARAN'S  STATE  SPACE  REALIZATIONS 

IIR  filter  realizations  are  possible  and  most  of  them  can  be  presented  within  the 
general  time  domain  representation  of  the  system  under  consideration.  In  continuous 
time  it  was  convenient  to  represent  a  linear  constant  coefficient  differential  equation  by 
a  system  of  first-order  differential  equations  for  analysis  purposes.  It  is,  therefore,  not 
surprising  that  it  is  useful  to  represent  a  linear  constant  coefficient  difference  equation 
by  a  system  of  first-order  linear  constant  coefficient  difference  equations.  Such  a  for- 
mulation results  in  a  special  form  of  representation  of  the  difference  equation  called  a 
state  variable  representation. 

A.     STATE  SPACE  REPRESENTATION  OF  A  NTH  ORDER  RECURSIVE 
DIGITAL  FILTER 

We  are  used  to  working  with  a  difference  equation  of  the  form 

y(n)  =  bxy(n  -  1)  +  •••  +  b^v{n  -  A")  +  a0x(n)  +  axx(n  -  1)  +  -  +  aLx{n  -  M)     (4. 1) 
or  in  the  notation 

N  M 

y(n)  =  Yj>&(n  -  k)  +  TJakA'1  -  k)  (4.2) 

k=\  k=0 

that  is  called  the  Nth-order  form.  An  alternative  is  to  describe  an  LTI  system  with  N 
first- order  difference  equations  and  an  output  equation.  This  alternative  form  for  de- 
scribing a  linear  time-invariant  system  is  known  as  the  state-space  description  and  can 
be  written  as  a  matrix  equation: 

v(«)  =  Xv(n  -  1 )  +  Bx{n)         (stale  equation)  (4.3) 

with  the  system  output  given  by 

y(n)  =  CT\(n)  +  dx{n)         (output  equation).  (4.4) 

In  Eqs.  (4.3)  and  (4.4),  \(n)  is  a  vector  with  N  rows  and  1  column  (an  A"  x  1  vector  or 
column  matrix)  of  newly  defined  variables  called  the  states,  x(n)  is  the  system  input,  and 
y(n)  its  output.  In  the  matrix  state  equation  of  Eq.  (4.3)  the  matrix  A  must  be  an 
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A'  x  N  matrix  and  the  B  matrix  an  N  x  1  vector.  The  scalar  output  y(n)  is  defined  by  Eq. 
(4.4)  where  C  is  a  N  x  1  column  matrix  and  i/is  simply  a  real  number  that  multiplies  (or 
scales)  the  input  x(n).  The  elements  in  A,  B,  C,  and  d  are  constants  that  depend  only 
upon  the  system  parameters. 

B.     WRITING  THE  STATE  SPACE  EQUATIONS 

To  illustrate  this  approach  to  modeling  of  LTI  systems  we  return  to  a  version  of  the 
Nth-order  difference  equation  described  as  a  second  order  II R  filter,  namely 

y(n)  =  bty(n  —  1)  +  b^n  —  2)  4-  aQx(n)  +  axx{n  —  1)  +  a2x(n  —  2)  (4.5) 

Y  =  bxz~x  Y  +  b2z~2  Y  +  a0X  +  a^z~]X  +  a2z~2X.  (4.6) 

We  define  the  state  varible  as:   Vu   V2.     Then,  we  can  establish  the  following  re- 
lationships: 


V\  =  y  =  V1  Y  +  b2z~2  Y  +  a0X  +  axz~x  X  +  a2z~2X 
Yl  — 2     v\  —  2     l 

-It/  -2  v 


(4.7) 


Therefore  the  state  equation  is: 


r,  =  b,z~]  l\  +  b,z~]  V0  +  anX+a,z~*X+a,z~2X 


1  T  t'2-        '  2  ^  M(P 


V  -7~lV 

\  2  —  Z        KL 


(4.8) 


We  express  it  in  the  matrix  state  equation: 


v(«)  =  Av(/7  -  1)  +  B.v(w) 


\  =  Az~l\  +  BX 


>r 

> 

v 

-1 
z 

>r 

+ 

_r2_ 

_i 

0_ 

_^2. 

<20  +  «i^        +  ^22 


-2 


X. 


(4.9) 
(4.10) 

(4.11) 


That  means: 


V  = 


(4.12) 
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B  = 


bx     b2 


-l 


an  +  a,z      +  a~,z 


-21 


(4.13) 


(4.14; 


The  output  equation  will  be: 


y{n)  =  Cy \{n)  +  dx{n) 


Y=lcx     c,D 


V, 


+  dX 


Y=c,V,  +ClV,  +  dX. 


1  '  1  ^  u2'  2 


(4.15) 
(4.16) 

(4.17) 


C.     MODIFIED  STATE  SPACE  EQUATION 

From  the  state  equation 

v(«)  =  Av(«  -  1)  +  Bx(n), 


(4.18) 


the  computation  ofv(n)  requires  the  immediate  past  state  vector  v(n-l).  This  recursive 
nature  severely  limits  the  speed  that  can  be  achieved  by  direct  implementation  of  the 
state  equation.  The  speed  limits  can  be  overcome  by  modifying  this  state  equation  so 
that  v(n)  depends  not  on  the  immediate  past  vector,  v(n  —  1),  but  on  v(n  —  L)  for  some 
L  >  1.  This  can  be  accomplished  as  follows: 
For  L  =  1  the  state  equation  is 


\{n)  =  A\{n  -  1)  +  Bx(n) 


-i 


V  =  A2"'V  +  BA\ 


Delaying  equation  (4.19)  by  1,  2,  and  3  yields 

v(«  -  1)  =  Av(«  -  2)  +  Bx(n  -  1) 
v(n  -  2)  =  Av(«  -  3)  +  Bx(n  -  2) 
v(/7  -  3)  =  Av(«  -  4)  +  Bx(n  -  3). 
For  L  =  2  putting  (4.21)  into  (4.22),  the  state  equation,  results  is 


(4.19) 
(4.20) 

(4.21) 
(4.22) 
(4.23) 
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v(«)  =  A[Av(«  -  2)  +  B.r(»  -  1)]  +  B.y(«)  (4.24) 

\(n)  =  A2\{n  -  2)  +  ABx(n  -  1)  +  Bx(n)  (4.25) 

V  =  A2z~2V  +  ABz_,A'  +  BAr.  (4.26) 
For  L  —  3  putting  (4.22)  into  (4.25),  the  state  equation,  results  is 

v(«)  =  A2[Av(«  -  3)  +  B*(w  -  2)]  +  ABjc(«  -  1)  +  Bjc(h)  (4.27) 

v(/i)  =  A3v(«  -  3)  +  A2Bx(«  -  2)  +  ABx[n  -  1)  +  Bjc(/i)  (4.28) 

V  =  A3z~3V  +  A2Bz~2X+ABz~]X+BX.  (4.29) 
For  L  =  4  putting  (4.23)  into  (4.28),  the  state  equation,  results  is 

v(,7)  =  A3[Av(«  -  4)  +  Bx(n  -  3)]  +  A2B*(«  -  2)  +  ABx{n  -  1)  +  Bjc(«)  (4.30) 

v(«)  =  A4v("  -  4)  +  A3B.v(«  -  3)  +  A2Ba'(«  -  2)  +  AB.r(«  -  1)  +  B.y(«)  (4.3 1) 

V  =  A4z~4V  4-  A3Bz"3A'  +  A2Bz"2A'  +  ABz"1^  +  BA.  (4.32) 
Equation  (4.32)  can  be  written  in  the  following  form  for  any  L>  1: 


L-\ 

/   ' 

j=0 


ALz~L\  +  \  A'Bz_/X 


One  of  the  benefits  of  the  state  space  realization  is  that  if  the  original  filter  (L=  1) 
is  stable,  all  higher  order  realizations  are  also  stable,  and  furthermore,  AL  has  improved 
error  properties  over  A.  This  advantage  of  the  state  space  method  is  discussed  in  the 
paper  by  Lu,  Lee,  and  Messerschmitt  [Ref.  4]. 

D.     SECOND  ORDER  IIR  FILTER  STATE  SPACE  REALIZATIONS 

From  the  state  space  equation  of  the  second  order  IIR  filter  in  section  B  and  the 
modified  state  space  equation  in  section  C,  we  can  develop  techniques  for  high  speed 
implementation  of  second  order  IIR  filters  for  several  values  of  L.  This  section  shows 
how  to  apply  these  concepts  for  L  =   1,2,  3.  and  4: 
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1.     L  =   1 

Beginning  with  equation  (4.20), 


V  =  Az     V  +  BX 


and  substituting  for  V,  A,  and  B  from  section  B  yields: 


(4.34) 


v2_ 

= 

"*1 

1 

0_ 

-l 

z 

r           -i  ,      -2"i 

aQ  +  fljZ      4-  a2z 

+ 

0 

(4.35) 


X 


Beginning  with  the  output  equation,  equation  (4.17),  and  referring  to  Figure  23, 
we  note  that  the  outpuf  of  the  multiplier  a  is  dz~2X.  This  results  is  a  revised  output 
equation  of 


y=C]K  +(hV,  +  dz~'X. 


(4.36) 


Substituting  the  state  equation  into  the  output  equation,  results  in 


y«(Vi  +  c2>",^i 


+  b2c}z   '  V2 
+  OqCxX 
+  «jC]Z     A 


-2 


+  (fl2C!  +  d)z      X 


(4-37) 


-l 


r=(Vi +  <>>>> 

+  b2c1z~2Y 

+  aQC\X 


+  a,c,z     A 


-2 


+  (a2q  +  d)z     X 


(4.38) 


To  demonstrate  equality  with  the  Y  of  the  original  II R  filter  as  expressed  in  equation 
(3.23)  for  p  =  0,  let  c,  =  1,  c2  =  0,  and  d  =  0  .   Then  the  output  equation  is: 


-1  ir 


-2V 


-2  v 


r=^z"'r+^z     >'+«o^+«i^    X+a2z  'X. 


(4.39) 


Bv  substituting: 
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*] 

= 

*n 

bi 

= 

*>12 

a0 

= 

au 

ax 

= 

ax2 

a-)  =  a. 


(4.40) 


we  can  write  this  equation  as: 

Y=buz~]  Y+  bX2z~zY  +  axxX+ax2z~'X+ax2z~'X. 


(4.41) 


This  realization  is  shown  in  Figure  23  on  page  47. 
2.     L  =  2 

Beginning  with  the  state  equation  (4.26), 


V  =  A2z~2V  +  XBz~lX+  BX 


and  substituting  for  V,  A,  and  B  from  section  B  yields: 


(4.42) 


>r 

"*i 

W 

"6 

1        *21 

_-2 

K 

r 

1      o_ 

J       0_ 

JX2. 

+ 

~bx     b{ 
1       0_ 

r             -i          -2~i 

aQ  +  axz      +  a2z 

0 

+ 

«o  + 

(7,2 

0 

-2 

X 

[x 


(4.43) 


>r 

~b\ 

_v2_ 

h 

+ 

+ 

\b2 

o 

Z 

>r 

2 

ji. 

by  +  b~,     b,b 


b 


a0bx  +  axbxz      +  a2bxz 
a0  +  axz      +  a2z 

-2 


-2 


lx 


a0  +  axz      +  a2z 
0 


A' 


(4.44) 


>r 

"a? 

_F2_ 

_*1 

+ 

4-  b2     bxb2 


-2 


Fo 


1  __^  T  ' 

aQ  +  (a0bx  +  ax)z      +  (axbx  +  a2)z      +  a2bxz 

-l    .         -2   ,         -3 
OqZ      +  axz      +  a2z 


X. 


(4.45) 
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Figure  23.       Second  Order  IIR  Filter:  L  =  1 


47 


(4.47) 


Beginning  with  the  output  equation,  equation  (4.17),  and  referring  to  Figure  24, 
we  note  that  the  output  of  the  multiplier  d  is  dz~lX.  This  results  is  a  revised  output 
equation  of 

Y  =  Cj  Vy  +c2V2  +  dz~3X.  (4.46) 

Substituting  the  state  equation  into  the  output  equation,  results  in 

Y=[{b2  +  b2)cx+bxc2~]z-2Vx 
+  Lb\b2cx  +  b2c2~\z     V2 

+  C(«(A  +  al)Cl  +  aQc2~]z~  % 

+  (axbxcx  +  a2cx  +  axc2)z    X 
+  {a2bxcx  +  a2c2  +  d)z~  X 

Y=\_{b\  +  b2)cx+b,c2~\z-2Y 

+  LV72C1  +b2c2^z~2y 
+  Oq^X 

+  [(^i  +  a\)c\  +  aoci1z~]X 
+  {axbxcx  +  a2cx  +  axc2)z    X 

+  (a2bxcx  +  fl2c2  +  d)z~  X- 

To  demonstrate  equality  with  the  Y  of  the  augmented  II R  filter  as  expressed  in  equation 
(3.34)  for  p  =   1,  let  c,  =  1,  c2  =  0,  and  d=0.   Then  the  output  equation  is: 


(4.4S) 


By  substitutine: 


}'=  [b1  +  b2)z~2Y+  bxb2z~2Y  +  a0X+(a0bx  +  ax)z~xX 
+  (axbx  +  a2)z~lX  +  a2bxz~  X. 


b\   +  b2  =  b\\ 

bxb2  =  bn 
«o  =  a\  1 


(4.49) 


(4.50) 
a0bx  +  ax  —  an 

axbx  +  a2  =  ax3 

a2bx  =  au 


we  can  write  this  equation  as: 
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-3-,/ 


-1 


Y  =  buz     Y  +  £nz     Y+a„X+  a„z~'X 


'12 

+  al2z-2X+auz-*X. 


w 


(4.51) 


This  is  the  same  as  the  equation  (3.34)  for  p=  1  found  on  page  15.    This  realization  is 
shown  in  Figure  24  on  page  50. 
3.     L  =  3 

Beginning  with  the  state  equation  (4.29), 


V  =  A3z"3V  +  A7Bz~2X  +  ABz~lX  +  BX 


(4.52) 


and  substituting  for  V,  A,  and  B  from  section  B  yields: 


v7 


bl  +  b2  bxb2 


b, 


+ 


+ 


+ 


bx     b2 
1      0 


-3 


b\  +  b2     bxb2 


bx     b2 
1       0 


'i 

*o  ~r  ui 


a0  +  axz      +  a2z 
0 


-T 


-2 


X 


an  +  a,z      +  a0z 


z~]X 


-l 


a0  +  axz      +  a2: 
0 


X 


(4.53) 


W 


L"2j 


bx{b2  +  2b2)     b2(b2  +  b2) 


-3 


b\  +  b2  bxb2 

a{){b2  +  b2)  +  ax{b\  +  b2)z~]  +  a2(b2  +  b2)z 


+ 


+ 


+ 


-2 


a06,  +  Q\b\Z      +  a2bxz 


'■x 


a()bi  +  a}b}z     +  a2bxz 

.  -2 

+  a2z 


-2 


z-X 


-i 


Oq  4-  axz      +  a2z 


-v 


X 


(4.54) 
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iM> 


&4 


i> 


B 


[3 


H 


J! 


[0-6  [0-© 


Figure  24.       The  Augmented  IIR  Filter:  L  =  2 
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1^2  J 


bx(h\  +  7b2)    b2{b]  +  b2) 


b\  +  b2 


b,b 


\u2 


-3 


+ 


a0  +  (a0bx  +  ax)z     +  (a0bx  +  aQb2  +  a]b]  +  a2)z 
+  {a]bl  +  axb2  +  a2bx)z~  +  (a2bx  +  a2b2)z 
a0z~   +  (a0b]  +  ax)z~  +  (axbx  +  a2)z~   +  a2bxz~ 


X. 


(4.55) 


Beginning  with  the  output  equation,  equation  (4.17),  and  referring  to  Figure  25, 
we  note  that  output  of  the  multiplier  d  is  dz^X.  This  result  is  a  revised  output  equation 
of 


-4 


Y=c]V]+c2V2  +  dz~^X, 


Substituting  the  state  equation  into  the  output  equation,  results  in 


(4.56) 


)'=[>?  4-  2bxb2)C]  +  [b]  +  b2)c2]z-3V] 

+  [(b~b2  +  bl)c]+blb2c2]z-2V2 
+  aQcxX 

+  C(fl<A  +  fli)ci  +  uqc2^z~]x 

+  [(aQbx  +  a0b2  +  #,£>,  +  a2)cx  +  (ao^i  +  a\)c2~\z~  % 
4-  \_{a]b]  -f  a]b2  4-  a2b])cx  +  (ai^i  +  a2)c2]z~  % 
+  \_{a2b[  +  a2b2)C]  +  a2bxc2  +  d~\z~  X 


(4.57) 


Y=l(b2l+2b]b2)cl  +  (bl2  +  b2)c2~\z'3Y 

+  \_{b2xb2  +  b\)cx  +  bxb2c2~\z~*Y 

+  a0cxX 

+  [(a06,  +  ax)cx  +  aQc{]z~XX  (4.58) 

+  \_{%b]  +  %b2  +  axbx  +  a2)c]  +  {a0b]  +  ax)c2~\z~~  X 

+  [_(axbx  +  axb2  +  a2bx)cx  +  (axbx  +  a2)c2~]z~  X 

+  \_{a2b]  +  a2b2)cx  +  a2bxc2  +  d~\z    X. 

To  demonstrate  equality  with  the  Y  of  the  augmented  II R  filter  as  expressed  in  equation 
(3.45)  for  p  =  2,  let  Cj  =  1,  c2  =  0.  and  d=0.   Then  the  output  equation  is: 
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Y={b\  +  2bxb2)z~2Y+{b]b2  +  b\)z~*Y  +  %X  +  {aQb{  +  ax)z~lX 

+  \_ad{b]  +  b2)  +  axbx  +  a2']z~2X  +  [a^b2  +  b2)  +  a2b,~\z~2X  (4.59) 

+  a2(b2  +  b2)z~*X. 


Bv  substitutina: 


b\  +  2b\b2  =  &i  j 
b]  b2  +  b2  =  bx2 


%{b2  +  b2)  +  alb]  +a2  =  a13 


ax(b\  +  b2)  +  a2l\  =  au 
a2(b\  +  b2)  =  al5 


we  can  write  this  equation  as: 


Y=buz~3Y+bnz"iY+anX+anz~lX 
+  auz    X  +  a]lXz    A  +  axsz    X. 


(4.61) 


This  is  the  same  as  the  equation  (3.45)  for  p=  2  found  on  page  17.    This  realization  is 
shown  in  Figure  25  on  page  53. 
4.     L  =  4 

Beginning  with  equation  (4.32), 

V  =  A4z_4V  +  A3Bz~2X  +  A2B-T:A'  +  ABz~lX  +  BX  (4.62) 

and  substituting  for  V,  A,  and  B  from  section  B  yields: 
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Figure  25.       The  Augmented  IIR  Filter:  L  =  3 
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>r 

"*1 

{b]  +  2b2)     b2(L 

i  + 

lh) 

~b 

l 

W 

-4 

>r 

V2_ 

2 

bx  +  b2             bxb2 

J     o_ 

.^2. 

+ 

~bx(b]  +  2b2)     b2{b]  +  b2) 

2 

bx  +  b2             bxb2 

r     ,       -i  ,       -2i 

Oq  +  d\Z      +  a2i 
0 

z"lY 

b}  +  b2     bxb2 

Oq  +  axz      +  a2z 

—2  „ 

+ 

bx             b2 

0 

Z    ^ 

+ 

~bx     b2~ 

r            -i  ,       -2i 
Oq  +  axz      +  a2z 

0 

z"1* 

r     i      -1  ,      -2i 

a0  +  Q\Z      +  ^2Z 

+ 

n 

X 

(4.63) 


Ko 


fef  +  3^2  +  ^     bxb2(b2x  +  2b2) 


W  +  2b2) 


b2(bx2  +  b2) 


K, 


+ 


+ 


+ 


+ 


^{b2  +  262)  +  axbx{b2  +  2b2)z~]  +  a2bx(b2  +  2b2)z~2 
a0(b2  +  b2)  +  ax(b2x  +  b2)z']  +  a2{b\  +  b2)z~2 


a0(b2  +  b2)  +  ax{b2  +  b2)z   '  +  a2(b2  +  b2)z  2 
a0bx  +  axbxz      +  a2bxz 


a0bx  +  axbxz     +  a2bxz 

-1     ,  -2 

a0  +  axz      +  a2z 


-l 


X 


<2q  +  axz      +  a2z 
0 


X 


■x 


-3 


X 


(4.64) 


^1 


bx  +  2bxb2  +  b\ 
bx(b2  +  2b2) 


bxb2{b{  +  2b2) 
b2(b2  +  b2) 


V, 


+ 


-2 


Oq  +  (dQbx  +  ax)z      +  {aQb\  +  a^b2  +  axbx  +  a2)z 

2  2  3 

+  {a0bx  -f  2a0bxb2  +  axbx  +  tfj62  +  «2^i)z 
4-  (flj/?!  +  2axbxb2  +  a2b2)z     +  (a2bx  +  2a2bxb2)z~' 
a0z      +  (a0bx  +  ax)z      +  (aQbx  +  aQb2  +  axbx  +  a2)z 
+  (axbx  +  axb2  +  a2bx)z~   +  (a2bx  +  a2b2)z~' 


X. 


(4.65) 
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Beginning  with  the  output  equation,  equation  (4.17),  and  referring  to  Figure  26, 
we  note  that  the  output  of  the  multiplier  d  is  dz~sX.  This  results  is  a  revised  output 
equation  of 

Y=cx  Vl  +  c2  V2  +  dz~5X.  (4.66) 

Substituting  the  state  equation  into  the  output  equation,  results  in 


Y=  [>f  +  lb\b2  +  b\)cx  +  {b\  +  26162)c2>_V1 
+  [(b3b2  +  2bxb22)cx  +  [b\b2  +  bl)c2]z~4V2 


(4.67) 


+  %C\X 

+  L(a0bx  +  «i)c,  +  %c{]z~\X 

+  \_{%b2x  +  a0b2  +  axbx  +  a2)c]  +  (a0bx  +  ax)c2~]z~2X  ^ 

-f  \_{ciQbx  +  2a0bxb2  +  axbx  +  axb2  +  a2bx)cx  +  (a0b2  +  a0b2  +  axbx  +  a2)c2~\z~3X 

+  [_(axbx  +  2axbxb2  +  a2b2  +  a2b2)cx  +  {axb2  +  axb2  +  a^b^c-^z^X 

+  [(fl2^i  +  2a2b\b2)cx  +  {a2bx  +  a2b2)c2  +  d~]z~5X 

Y=  [(£?  +  lb\b2  +  b])cx  +  {b3x  +  2bxb2)c2~]z-4Y 
+  [(bx3b2  +  2bxb\)cx  +  (bxb2  +  6:2)c2]-_5r 
+  «oci^ 

+  [(^-j/'i  +  a1)cl  +  <70c2];-1A' 

+  [(a0bx  +  a0b2  +  axbx  +  a2)cx  +  (a0bx  +  ax)c2~\z~^X 

4-  \_{aobx  +  2a0bxb2  +  axb\  +  axb2  +  a2bx)cx  +  {a0bx  +  Oq^  +  axbx  +  a2)c2\z     X 
+  [_{axbx  +  2axbxb2  +  a2bx  +  a2b2)cx  +  (axbx  +  axb2  +  a2bx)c2~\z~  X 
+  \_(o2bx  +  2a2bxb2)cx  +  {a2bx  +  a2b2)c2  +  d\z~  X. 

To  demonstrate  equality  with  the  Y  of  the  augmented  II R  filter  as  expressed  in  equation 
(3.56)  for  p  =  3,  let  c,—  1,  c2  =  0,  and  d=  0.   Then  the  output  equation  is: 


Y  =  (bx  +  3b2b2  +  bl)z~A Y  +  (b3xb2  +  2bxb\)z~5 Y  +  a0X  +  {a^bx  +  ax)z   XX 
+  \_%{bx  +  b2)  +  axbx  4-  a2~\z~  X 
+  [a0(b3x  +  2bxb2)  +  ax(b2  +  b2)  +  a2bx~]z-3X 
+  [_ax{b\  +  2bxb2)  +  a2(b2  +  b2)']z~~AX  +  a2{b\  +  2bxb2)z~5X. 


(4.69) 


Bv  substituting: 


55 


b}  -f  3b^b2  +  b2  —  bn 
b{b2  +  2b-[b2  =  b]2 

«0^1  +  a\  =  a\2 

a0(b]  +  b2)  +  albl  +  a2  —  a13 
a0(bl  +  2bxb2)  +  ax{bx  +  b2)  +  a2bx  =  a14 

3  2 

a,  (&!  +  2bxb2)  +  a2(b1  +  b2)  =  a]5 
a2{b\  +2b]b2)  =  a]6 


(4.70) 


we  can  write  this  equation  as: 


Y=b]]z~4Y  +  b^Y+a^X+a^z^X 


'\2< 
-2  v   ,  -3  v 

+  al5z~^X  +  al6z~  X. 


mi- 


'12' 


(4.71) 


This  is  the  same  as  the  equation  (3.56)  for  p=3  on  page  20.    This  realization  is  shown 
in  Figure  26  on  page  57. 

From  equation  (4.33)  we  can  see  that  the  feedback,  portion  of  the  realization 
must  have  a  delay  of  L  stages  in  any  loop.  The  general  realization  of  that  portion  is 
shown  in  Figure  27  on  page  58. 

E.     PIPELINED  IMPLEMENTATION 

All  of  the  1 1 R  filters  in  section  D  can  be  developed  into  high-speed  IIR  filters  by 
applying  pipeline  techniques.  We  use  Ka-stage  pipeline  add  units  and  Km-stage  pipeline 
multipliers  to  increase  the  speed  of  the  filters.  The  Km  stages  are  involved  in  the  coeffi- 
cient multiplies  and  the  Ka  stages  are  involved  in  the  adds.  The  net  result  is  a  (Km  + 
Ka)-stage  pipeline  multiply-add  unit.  The  equations  for  the  filter  outputs  for  different 
L's  are  given  below: 

1.     L  =   1 

The  equation  is 


-(8A'a+2A'mh,  _    -(&Ka+2Km) 


buz-{2Ka+Km)Y+bnz-QKa+Km+l)Y 
+  auX+  a122    X  +  anz    X 


(4.72) 


This  realization"  is  shown  in  Figure  28  on  page  59. 
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Figure  26.       The  Augmented  IIR  Filter:  L  =  4 
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Figure  27.       Heart  of  the  Pipeline  Recursive  Filter  Realization 

2.     L  =  2 

The  equation  is 


-(9Ka+2Km)  y  _     -(9Ka+2Km) 


,        -(2Ka+Km)v       ,        -(2Ka+Km+\)v 
°\\z  l  +  D\2Z  l 

+  a,iX+  a,->z    X  +  a,->z  ^X+a^z    X 


'11^  ^  u12' 


'14- 


(4-73) 


This  realization  is  shown  in  Figure  29  on  page  60. 
3.     L  =  3 

The  equation  is 


-(\0Ka+2Km)  y  _     -(\QKa+2Km) 


b„z 


-(2Ka+Km)y  ,,        -{2Ka+Km+\)y 


'11' 

+  a-nX+  a^z    X  +  a,-,z    X 


•if1  t  "12' 

+  tf^z-  >V+  a15z     A' 


(4.74) 


This  realization  is  shown  in  Figure  30  on  page  61. 
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Figure  28.       Second  Order  IIR  Filter:  L  =  1  with  Ka,  Km 
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Figure  29.       The  Augmented  I1R  Filter:  L  =  2  with  Ka,  Km 
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Figure  30.       The  Augmented  IIR  Filter:  L  =  3  with  Ka,  Km 
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4.     L  =  4 

The  equation  is 


-(\\Ka+2Km)v 


Y  = 


-(\\Ka+2Km) 


-(2Ka+Km) 


Y+bnz 


-(2Ka+Km+\) 


-2 


Y 


+  auX  +  anz    X  +  al3z    X 
+  auz~X+al5z    X+a]6z~X 


(4.75) 


This  realization  is  shown  in  Figure  31  on  page  63. 

From  analysis  of  the  feed  back  loop  portion  of  the  state  space  realization,  as 
shown  in  Figure  27  and  equation  (4.33),  we  have  seen  that  the  delay  around  any  loop 
must  be  L.  A  general  realization  for  the  feedback  portion  of  the  pipeline  case  is  shown 
in  Figure  32  on  page  64. 

Therefore  we  can  see  that 

L  =  2Ka+  Km.  (4.76) 

The  following  examples  show  the  analysis  of  the  cases  for  L  =  1,  2,  3  and  4. 
5.     Example  for  L  =   1: 
a.      Ka  =  0,  Km  =  1 

L  =  2Ka  +  Km 
=  2(0)  +  1 
=  1 

The  only  integer,  non-negative  values  for  Ka  and  Km  which  solve  equation  (4.76)  are 
Ka  =  0.  Km  =  1.  This  means  we  cannot  use  pipeline  modules  with  longer  delays  than 
Ka  =  0.  Km  =  1  . 

This  realizes  the  original  11R  filter  for  L  =   1,  but  with  an  additional  2  de- 
lays and  the  equation  is 


z  2Y=  z~\bViz~]Y  +  b]2z~2Y+  anX  +  anz~lX+  a]3z~2X~]. 


(4.77) 


This  realization  is  shown  in  Figure  33  on  page  65. 
6.     Example  for  L  =  2: 

The  only  integer,  non-negative  values  for  Ka  and  Km  which  solve  equation 
(4.76)  are  Ka  =  0,  Km  =  2  or  Ka  =  1,  Km  =  0. 
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Figure  31.       The  Augmented  IIR  Filter.  L  =  4  with  Ka,  Km 
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Figure  32.       Heart  of  More  General  Pipeline  Filter  Realization 

a.  Ka  =  0,  Km  =  2 

L  =  2Ka  +  Km 
=  2(0)  +  2 
=  2 

For  these  values,  we  use  pipeline  modules  with  delays  Ka  =  0,  Km  =  2  . 

This  realizes  the  augmented  IIR  filter  for  L  =  2,  but  with  an  additional  4 
delays.  The  equation  is 

z^Y  =  z~*[bnz-2Y  +  bnz~3  Y  +  auX  +  anz~XX  +  auz~2X  +  a14z~3A'].      (4.78) 

This  realization  is  shown  in  Figure  34  on  page  66. 

b.  Ka  =  1,  Km  =  0 

L  =  2Ka  +  Km 
=  2(1)4-0 

=  2 

For  these  values,  we  us~  pipeline  modules  with  delays  Ka  =  1,  Km  =  0  . 
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Figure  33.       Second  Order  IIR  Filter:  L  =  1  with  Ka  =  0,  Km  =  1 
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Figure  34.       The  Augmented  II R  Filter:  L  =  2  with  Ka=  0,  Km  =  2 
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This  realizes  the  augmented  IIR  filter  for  L  =   2,  but  with  an  additional  9 
delays.  The  equation  is 

z-9Y  =  z~9lbu2~2Y+bnz~3Y+anX+an2~lX+auz~2X+auz~2X'].      (4.79) 

This  realization  is  shown  in  Figure  35  on  page  68. 
7.     Example  for  L  =  3: 

The  only  integer,  non-negative  values  for  Ka  and  Km  which  solve  equation 
(4.76)  are  Ka  =  0,  Km  =  3  or  Ka  =  1,  Km  =  1. 
a.      Ka  =  0,  Km  =  3 

L  =  2Ka  +  Km 
=  2(0)  +  3 
=  3 

For  these  values,  we  use  pipeline  modules  with  delays  Ka  =  0,  Km  =  3  . 

This  realizes  the  augmented  IIR  filter  for  L  =  3,  but  with  an  additional  6 
delays.  The  equation  is 


Y=z~° 


-3 


6i-,z  ~Y+b:->:     Y+  auX  +  a^z    X 


+  al3z  ~X+al4z    X+a^z    X 


5^ 


(4. SO) 


This  realization  is  shown  in  Figure  36  on  page  69. 
b.      Ka  =  1,  Km  =  1 

L  =  2Ka  +  Km 
=  2(1)+1 

=  3 

For  these  values,  we  use  pipeline  modules  with  delays  Ka  =  1,  Km  =  1  . 

This  realizes  the  augmented  IIR  filter  for  L  =  3,  but  with  an  additional  12 
delays.  The  equation  is 


CnY=z~u 


bnz~*Y+bnz~*Y+anX+anz  ]X 
+  a,-,z    X  +  a]Az    X+a^z     X 


'13' 


'14' 


(4.81) 


This  realization  is  shown  in  Figure  3  /  on  page  70. 
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Figure  35.       The  Augmented  IIR  Filter:  L  =  2  with  Ka=  1,  Km=  0 
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Figure  36.       The  Augmented  IIR  Filter:  L  =  3  with  Ka=  0,  Km=  3 
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Figure  37.       The  Augmented  II R  Filter:  L  =  3  with  Ka=  1,  Km=  1 
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8.     Example  for  L  =  4: 

The  only  integer,  non-negative  values  for  Ka  and  Km  which  solve  equation 
(4.76)  are  Ka  =  0,  Km  =  4  and  Ka  =  1 ,  Km  =  2  or  Ka  =  2,  Km  =  0. 
a.      Ka  =  0,  Km  =  4 


L  =  2Ka  +  Km 
=  2(0)  +  4 
=  4 

For  these  values,  we  use  pipeline  modules  with  delays  Ka  =  0,  Km  =  4  . 

This  realizes  the  augmented  IIR  filter  for  L  =  4,  but  with  an  additional  8 
delays.  The  equation  is 


-5  V 

Z        1    = 


buz  4Y+bnz   5Y+  auX+anz   lX 

+  al3z~2X  +  auz~3X  +  al5z~*X  +  al6z~sX 


(4.82) 


This  realization  is  shown  in  Figure  3S  on  page  72. 
b.      Ka  =  1,  Km  =  2 

L  =  2Ka  +  Km 
=  2(1) +  2 

=  4 


For  these  values,  we  use  pipeline  modules  with  delays  Ka=  \,  Km  =  2  . 

This  realizes  the  augmented  IIR  filter  for  L  =  4  but  with  an  additional  15 
delays.  The  equation  is 


r-15r=z-15 


buz~AY+bnz~5Y+  auX+  a]2z~]X 

+  a,iZ~"X  +  a]Az~  X  +  a,cz     X  +  a^z~  X 


(4.83) 


This  realization  is  shown  in  Figure  39  on  page  73. 
c.      Ka  =  2,  Km  =  0 

L  =  2Ka  +  Km 
=  2(2)  +  0 
=  4 

For  these  values,  we  use  pipeline  modules  with  delays  Ka  =  2,  Km  =  0  . 

This  realizes  the  augmented  IIR  filter  for  L  =  4,  but  with  an  additional  22 
delays.  The  equation  is 
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Figure  38.       The  Augmented  IIR  Filter:  L  =  4  with  Ka=  0,  Km  =  4 
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Figure  39.       The  Augmented  IIR  Filter:  L  =  4  with  Ka=  1,  Km  =  2 
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-22 ,,  -22 

2  1  =  Z 


-4 


-5 


-1 


Y+  bnz     Y+auX+anz    X 


-2 


-3 


-4 


-5 


+  auz    X +  auz    X  +  ax5z    X  +  axbz    X 


(4.84) 


This  realization  is  shown  in  Figure  40  on  page  75. 


Table  5. 

STATE  SPACE  REALIZATION  SUMMARY 

L 

Ka-add  units 

Km-multipliers 

D(output  delayed) 

Figure 

1 

0 

1 

2 

31 

2 

0 

2 

4 

32 

2 

1 

0 

9 

33 

-> 

0 

3 

6 

34 

3 

1 

1 

12 

35 

4 

0 

4 

8 

36 

4 

1 

2 

15 

37 

4 

2 

0 

22 

38 

F.     SUMMARY 

This  chapter  develops  an  alternate  method  for  high  speed  implementation  of  state 
space  digital  filters.  The  proposed  method  implements  AL  instead  of  A  as  seen  in 
equation  (4.33).  The  proposed  method  allows  arbitrary  A  in  the  implementation,  and  it 
can  take  full  advantage  of  the  optimum  state  space  filters  along  with  the  benefit  of  A^- 
realization.  The  delay  and  complexity  of  each  realization  are  summarized  in  Table  5. 
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Figure  40.       The  Augmented  IIR  Filter:  L  =  4  with  Ka=2,  Km=0 
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V.     CONCLUSIONS 

A.  SUMMARY 

This  thesis  describes  pipeline  realizations  of  the  1IR  notch  filter.  The  transfer 
function  and  the  coefficients  that  are  appropiate  for  IIR  notch  filter  have  been  pre- 
sented. Realization  of  a  high  speed  IIR  notch  filter  using  a  pipeline  technique  by  Loomis 
and  Sinha  has  been  presented.  An  alternate  technique  in  using  state  space  realization 
by  R.  Gnanasekaran  has  also  been  presented. 

B.  COMPARISON 

Table  4  in  Chapter  III  summarizes  the  Loomis,  Sinha  pipeline  recursive  design  and 
Table  5  in  Chapter  TV  summarizes  Gnanasekaran's  state  space  realization.  It  was  shown 
in  Chapter  III  section  I  that  we  have  only  one  choice  of  Ka  and  Km  for  each  value  of 
p  and  in  Chapter  IV  section  E  that  we  have  1  choice  of  Ka  and  Km  for  L=  1,  2  choices 
of  Ka  and  Km  for  L=  2  and  for  L=  3,  and  3  choices  of  Ka  and  Km  for  L  =  4. 

In  the  Loomis/Sinha  pipeline  recursive  design  we  note  that  for  p=  2,  Ka=  1,  Km=  2, 
output  delay  =  5.  the  filter  is  unstable.  However,  in  Gnanaskaran's  state  space  realization 
for  L  =  4,  Ka=  1.  Km  =2.  output  delay  =  15,  the  filter  is  stable.  Thus  we  see  that  in  cases 
where  Loomis,  Sinha  realization  are  unstable,  a  state  space  realization  always  exists  and 
is  stable.  However,  it  will  require  more  hardware  and  have  a  greater  overall  delay.  All 
of  these  results  are  summarized  in  Table  6.  In  that  table,  the  column  labelled  D  refers 
to  the  overall  delay  and  the  columns  labelled  #Mlt  and  #Add  refer  to  the  number  of 
multipliers  and  number  of  adders  as  shown  in  the  appropriate  figure. 

Note  that  for  a  given  value  of  Ka  and  Km,  the  number  of  adders  and  multipliers 
required  for  the  state  space  realization  is  larger  than  that  required  for  the  scalar.  This 
disadvantage  is  offset  by  the  guaranteed  stability  of  the  state  space  realization. 
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Table 

6.     SUMMARY  OF  REALIZATIONS 

State  Space  Realization 

Scalar  Realization 

Ka 

Km 

L 

D 

Fig 

#Mlt 

#Add 

P 

D 

Fig 

#Mlt 

#Add 

0 

1 

1 

2 

31 

13 

11 

0 

2 

2 

4 

32 

15 

13 

1 

0 

2 

9 

33 

15 

13 

0 

3 

19 

5 

5 

0 

6 

34 

17 

15 

1 

1 

-> 

j 

12 

35 

17 

15 

1 

4 

20 

6 

6 

0 

4 

4 

8 

36 

19 

17 

1 

2 

4 

15 

37 

19 

17 

2 

5 

21 

7 

7 

2 

0 

4 

22 

38 

19 

17 

1 

** 

j 

6 

22 

8 

S 
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APPENDIX  A.     PROGRAM  AND  DATA  OF  FREQUENCY  RESPONSE 

A.  PROGRAM  FOR  FREQUENCY  RESPONSE  CALCULATION 

C  DIGFREQ.  FOR              VERSION:  2/03/88 

C 

C 

C  PURPOSE:   THIS  PROGRAM  COMPUTES  THE  FREQUENCY  RESPONSE  OF 

C  DISCRETE  SYSTEMS.   THE  PROGRAM  CONSISTS  OF  A  MAIN 

C  PROGRAM  THAT  CONTROLS  THE  INPUT/OUTPUT  AND  THE 

C  SUBROUTINES  dfresp  AND  coeff.   SUBROUTINE  dfresp 

C  COMPUTES  THE  FREQUENCY  RESPONSE  OF  EACH  SYSTEM. 

C  SUBROUTINE  coeff  ALLOWS  THE  USER  THE  OPTION  OF 

C  GENERATING  THE  FILTER  COEFFICIENTS  OF  THE  SYSTEMS 

C  TO  BE  ANALYZED  BY  WRITING  THE  APPROPRIATE  EQUATIONS. 

C  IF  THE  USER  ELECTS  TO  GENERATE  THE  COEFFICIENTS  BY 

C  USING  THE  SUBROUTINE  coeff,  THE  EQUATIONS  MUST  BE 

C  WRITTEN  INTO  THE  SUBROUTINE  USING  STANDARD  FORTRAN  77 

C  STATEMENTS.   THE  COEFFICIENTS  MUST  BE  STORED  IN  THE 

C  ARRAYS  b()  AND  c( )  WHICH  CORRESPOND  RESPECTIVELY  TO  THE 

C  NUMERATOR  AND  DENOMINATOR  TERMS  OF  THE  SYSTEM  EQUATION. 

C  THE  USER  CAN  SELECT  ONE  OF  TWO  OPERATING  MODES:  BATCH 

C  OR  TEST.   IN  BATCH  MODE  THE  AMOUNT  OF  INTERFACE  WITH 

C  THE  USER  IS  MINIMIZED  AND  IT  IS  ASSUMED  THAT  THE  INPUT 

C  DATA  HAS  BEEN  STORED  IN  THE  DEFAULT  FILE  ' DIGFREQ. IN' . 

C  IN  TEST  MODE  THE  USER  IS  PROMPTED  FOR  THE  NAME  OF  THE 

C  INPUT  FILE  OR  HAS  THE  OPTION  TO  PERFORM  A  TRIAL  RUN  BY 

C  USING  THE  INPUT  DATA  STORED  IN  THE  FILE  ' DIGFREQ.  TST1 . 

C  IT  IS  RECOMMENDED  THAT  FIRST-TIME  USERS  SELECT  THE  TEST 

C  MODE  AND  MAKE  A  TRIAL  RUN  WITH  THE  PRESTORED  INPUT  DATA. 

C  THE  TEST  MODE  ECHOES  THE  INPUT  DATA  ONTO  THE  MONITOR  TO 

C  ALLOW  VERIFICATION  OF  ITS  ACCURACY.   THIS  PROGRAM  WILL 

C  COMPUTE  THE  FREQUENCY  RESPONSE  OF  UP  TO  THREE  SYSTEMS. 

C  FOR  EACH  SYSTEM,  THE  USER  HAS  THE  OPTION  OF  HAVING  THE 

C  OUTPUT  EXPRESSED  IN  DECIBELS  (db).   THE  OUTPUT  OF  THIS 

C  PROGRAM  IS  STORED  IN  TABULAR  FORM  IN  THE  FILE 

C  'DIGFREQ. OUT'  AND  IN  A  FORM  SUITABLE  FOR  PLOTTING 

C  IN  THE  FILE  'DIGFREQ.DAT'. 

C 

C 

c 

C  THIS  PROGRAM  ASSUMES  THAT  EACH  DISCRETE  SYSTEM  IS  MODELED  BY  THE 

C  EQUATION:   H(z)  =  num/den  WHERE: 

C 

C  num  =  b(0)*z**T.  -f  b(l)*z**(L-l)  +  ...  +  b(L-l)*z  +  b(L) 

C 

C  den  =  c(0)*z**N  +  c(l)*z**(N-l)  +  ...  +  c(N-l)*z  +  c(N) 

C 

C  L  =  A  NON-NEGATIVE  INTEGER,  THE  DEGREE  OF  THE  NUMERATOR 

C  POLYNOMIAL. 

C  N  =  A  NON-NEGATIVE  INTEGER,  THE  DEGREE  OF  THE  DENOMINATOR 
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POLYNOMIAL. 
b(0)...b(L)  =  REAL  COEFFICIENTS  OF  THE  NUMERATOR  TERMS. 
c(0)...c(N)  =  REAL  COEFFICIENTS  OF  THE  DENOMINATOR  TERMS. 

THE  INPUT  PARAMETERS  SHOULD  BE  STORED  IN  A  FILE  NAMED 
'DIGFREQ. IN'.   ALL  OF  THE  READ  STATEMENTS  USED  BY  THIS  PROGRAM 
REQUIRE  FORMATTED  INPUT.   PARTICULAR  ATTENTION  SHOULD  BE  PAID 
TO  THE  FORMATS,  ESPECIALLY  THE  USE  OF  THE  DECIMAL  POINT  TO 
DENOTE  'REAL'  NUMBERS.   THE  INPUT  PARAMETERS  REQUIRED  BY  THE 
PROGRAM  ARE  LISTED  BELOW. 


NAME 

nurasys 

L 

N 

dsorce 

yscal 

thetaO 

dithta 

nurapts 

b() 

c() 

WHERE: 

numsvs 


TYPE 

INTEGER 

INTEGER 

INTEGER 

CHARACTER 

CHARACTER 

REAL 

REAL 

INTEGER 

REAL 

REAL 


RANGE  (ARRAYS) 


0,  1,  2, 
0,  1,  2, 


,    L 

,  N 


RESTRICTIONS 

1  <=  nurasys  <=  3 
0  <=  L  <=  128 
0  <=  N  <=  128 
'F'  OR  'S' 
'STD1  OR  'LOG' 


1  <=  numpts  <=  101 
0  <=  L  <=  128 
0  <=  N  <=  128 


THE  NUMBER  OF  DISTINCT  SYSTEMS  H(z)  TO  BE  ANALYZED. 
THIS  INTEGER  VALUE  MUST  OCCUR  AT  THE  TOP  OF  THE  INPUT 
FILE.  IT  DELINEATES  THE  NUMBER  OF  SYSTEMS  TO  BE  READ  BY 
THE  PROGRAM  AND  ANALYZED.  FOR  EACH  SYSTEM  (1,  ...  ,  nurasys) 
THE  PARAMETERS  BELOW  MUST  APPEAR  IN  THE  INPUT  FILE. 


L  =  AN  INTEGER  VALUE  SPECIFYING  THE  DEGREE  OF  THE  NUMERATOR 
POLYNOMIAL. 

N  =  AN  INTEGER  VALUE  SPECIFYING  THE  DEGREE  OF  THE  DENOMINATOR 
POLYNOMIAL. 

dsorce  =  THE  CHARACTER  STRING  'F'  OR  'S1  DENOTING  WHETHER  THE 

SYSTEM  COEFFICIENTS  ARE  TO  BE  READ  FROM  THE  INPUT  FILE  (F) 
OR  GENERATED  (S)  THROUGH  USE  OF  THE  SUBROUTINE  coeff. 

yscal  =  A  CHARACTER  STRING  SPECIFYING  THE  DESIRED  MAGNITUDE  OPTION: 
'STD'  WILL  PRODUCE  STANDARD  MAGNITUDE  OUTPUT; 
'LOG'  WILL  PRODUCE  MAGNITUDE  EXPRESSED  IN  DECIBELS  (db). 

thetaO  =  THE  STARTING  VALUE  OF  THETA  (RADS)  AS  IN  Z=EXP( J*THETA). 

dithta  =  THE  INCREMENT  OF  THETA  (RADIANS). 

numpts  =  THE  NUMBER  OF  FREQUENCY  POINTS  FOR  WHICH  THE  OUTPUT  IS 
TO  BE  COMPUTED. 

b()  =  THE  NUMERATOR  COEFFICIENTS  IN  ORDER  b(0),  b(l),  ...,  b(L). 
IF  dsorce  =  'F'  IS  SELECTED  THEN  THE  USER  MUST  SUPPLY  THE 
L+l  NUMERATOR  COEFFICIENTS  IN  THE  FILE.   IF  dsorce  =  'S' 
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THEN  THE  USER  HAS  ELECTED  TO  GENERATE  THE  NUMERATOR 
COEFFICIENTS  BY  WRITING  THE  APPROPRIATE  FORTRAN  STATEMENTS 
IN  THE  SPACE  PROVIDED  IN  SUBROUTINE  coeff.   IF  THIS  METHOD 
OF  DATA  GENERATION  IS  ELECTED  THE  PROGRAM  MUST  BE  RECOMPILED 
BEFORE  EXECUTION. 

c()  =  THE  DENOMINATOR  COEFFICIENTS  IN  ORDER  c(0),  c(l),  ...,  c(N). 
IF  dsorce  =  'F'  IS  SELECTED  THEN  THE  USER  MUST  SUPPLY  THE 
N+l  DENOMINATOR  COEFFICIENTS  IN  THE  FILE.   IF  dsorce  =  'S' 
THEN  THE  USER  HAS  ELECTED  TO  GENERATE  THE  DENOMINATOR 
COEFFICIENTS  BY  WRITING  THE  APPROPRIATE  FORTRAN  STATEMENTS 
IN  THE  SPACE  PROVIDED  IN  SUBROUTINE  coeff.   IF  THIS  METHOD 
OF  DATA  GENERATION  IS  ELECTED  THE  PROGRAM  MUST  BE  RECOMPILED 
BEFORE  EXECUTION. 

NOTE:   THE  INPUT  FORMAT  STATEMENTS  OCCUR  IN  THE  MAIN  PROGRAM 

FOLLOWING  THE  CAPTION:  ********   INPUT  FORMAT  ********. 

THE  FORM  OF  THE  INPUT  DATA  FILE  IS: 


LINE  # 


1 
2 
3 

NOTE 
NOTE 
NOTE 


ENTRIES 

numsys 

L,N,dsorce,yscal 

dlthta, theta0,numpts 

b(k),  k=0,l,... ,L 

c(k),  k=0,l,... ,N 


FORMAT 

II 
i3,tll,i3,t21,al,t31,a3 
2fl0. 0,i3 
6fl0.0 
6fl0.0 


WHERE:   NN  =  1  +  (L/6  ROUNDED  DOWN  TO  THE  NEXT  SMALLER  INTEGER). 
ND  =  1  +  (N/6  ROUNDED  DOWN  TO  THE  NEXT  SMALLER  INTEGER), 


NOTES  1. 


2. 


THE  NEXT  NN  LINES  ARE  ONLY  REQUIRED  IF  dsorce  =  'F'.  IF 
dsorce  =  'S'  THEN  THE  USER  HAS  ELECTED  TO  GENERATE  THE 
L+l  NUMERATOR  COEFFICIENTS  IN  THE  SUBROUTINE  coeff. 
THE  USER  MUST  PROVIDE  THE  APPROPRIATE  FORTRAN  STATEMENTS 
IN  SUBROUTINE  coeff  TO  GENERATE  THE  VALUES  FOR  b(). 

THE  NEXT  ND  LINES  ARE  ONLY  REQUIRED  IF  dsorce  =  *F'.  IF 
dsorce  =  'S'  THEN  THE  USER  HAS  ELECTED  TO  GENERATE  THE 
N-rl  DENOMINATOR  COEFFICIENTS  IN  THE  SUBROUTINE  coeff. 
THE  USER  MUST  PROVIDE  THE  APPROPRIATE  FORTRAN  STATEMENTS 
IN  SUBROUTINE  coeff  TO  GENERATE  THE  VALUES  FOR  c(). 

FOR  numsys  >  1  THE  FORMAT  OF  LINES  2. . .  IS  REPEATED. 

THE  FORMAT  f 10. 0  USED  FOR  INPUT  DATA  PERMITS  THE  DECIMAL 
POINT  TO  BE  PLACED  ANYWHERE  IN  THE  FIELD  OF  10  COLUMNS 
AND  ALSO  ALLOWS  THE  EXPONENTIAL  FORMAT  TO  BE  USED  (EG. 
3146.  2  =  3.  1462E+03). 


c****************************  OUTPUT  ******************************* 

THE  OUTPUT  DATA  CREATED  BY  THE  PROGRAM  IS  STORED  IN  TABULAR  FORM 
IN  THE  FILE  'DIGFREQ. OUT' .   ADDITIONALLY,  THE  OUTPUT  DATA  IS 
WRITTEN  INTO  THE  FILE  'DIGFREQ.DAT'  TO  FACILITATE  PLOTTING  BY 
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C  A  SEPARATE,  USER  SUPPLIED  PROGRAM.   THE  FORMAT  OF  THE  DATA  IN 

C  'DIGFREQ.DAT'  IS:   el2.6,  2x,  el2.  6.   THE  FIRST  ENTRY  CORRESPONDS 

C  TO  THE  ORDINATE  VALUE  (THETA)  AND  THE  SECOND  ENTRY  THE  ABSCISSA 

C  VALUE  (MAGNITUDE  OR  PHASE).   ADDITIONAL  HEADER  INFORMATION  IS 

C  WRITTEN  INTO  THE  DATA  FILE  TO  ALLOW  FOR  CONTROL  AND  LABELING  OF 

C  EACH  PLOT. 

C 

C 

Q-kir-fcJe-it-k-kMr-k-k-kit-ir-k-kJv-k-jtlrkitirk-kjrjtJrJrk       EXAMPLE   ****************************** 

c 

C  THE  INPUT  PARAMETERS  FOR  THE  SYSTEM  DESCRIBED  BELOW  ARE  STORED  IN 

C  THE  SAMPLE  INPUT  FILE  ' DIGFREQ. TST'  AND  CAN  BE  USED  FOR  A  TRIAL 

C  RUN  IN  THE  TEST  MODE. 

C 

C 

C   SYSTEM:   H(z)=z/(z-0. 5) 

C 

C   GOAL:     TO  OBTAIN  THE  FREQUENCY  RESPONSE  FOR  THIS  SYSTEM  FROM 

C  THETA  =  0.0  TO  THETA  =  3.14159  (PI  RADIANS)  IN  STEPS 

C  OF  dlthta  =  PI/10.  0 

C 

C   FOR  THE  SYSTEM  DESCRIBED  ABOVE  THE  INPUT  FILE  IS: 

C 

C   1 

C   001       001       F         STD 

C   0.314159   0.0       Oil 

C   1.  0       0.  0 

C   1.0       -0.5 

C 

C 

C  THE  RESULTING  OUTPUT  DATA  FILE:  'DIGFREQ. OUT'  IS: 

C 

C  INPUT  DATA  FOR  SYSTEM  #   1 

C 

C   INPUT  DATA  SOURCEFILE:  DIGFREQ.  TST 

C  DEGREE  OF  NUMERATOR  =   1 

C   DEGREE  OF  DENOMINATOR  =    1 

C  dsorce  =  F 

C  NUMBER  OF  FREQUENCY  POINTS  =   11     MAGNITUDE  OPTION  =  STD 

C   STARTING  VALUE  OF  THETA  =   . O0000OE+00 

C   INCREMENT  OF  THETA  =   . 314159E+00 

C 

C  THE  NUMERATOR  COEFFICIENTS  b(0) ,b( 1). .  .  b(L)  ARE 

C 

C    . 1000E+01     . 0000E+00 

C 

C 

C  THE  DENOMINATOR  COEFFICIENTS  c(0) ,c( 1). . . c(N)  ARE 

C 

C    .  1000E+01    -. 5000E+00 

C 

C 

C  OUTPUT  DATA  FOR  SYSTEM  #   1 

C 

C         THETA         MAGNITUDE         PHASE 

C       (RADIANS)  (DEGREES) 
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OOOOOOE+00 
314159E+00 
628318E+00 
942477E+00 
125664E+01 
157080E+01 
188495E+01 
219911E+01 
251327E+01 
282743E+01 
314159E+01 


200000E+01 
182897E+01 
150588E+01 
122886E+01 
103088E+01 
894428E+00 
800894E+00 
737654E+00 
696900E+00 
674038E+00 
666667E+00 


.  OOOOOOE+00 
-.  164149E+02 
-.262677E+02 
-. 29807E+02 
-. 293546E+02 
-.265651E+02 
-. 223862E+02 
-. 173608E+02 
-. 118186E+02 
-.597793E+01 
-.484184E-04 


END  OF  RUN,  SYSTEM  #  1 


FOR  ILLUSTRATIVE  PURPOSES  THE  COEFFICIENTS  b( )  AND  c( )  COULD 
HAVE  BEEN  GENERATED  BY  SPECIFYING  dsorce  =  'S'  AND  WRITING  THE 
APPROPRIATE  FORTRAN  STATEMENTS  INTO  SUBROUTINE  coeff.   THE 
STATEMENTS  THAT  COULD  BE  USED  TO  ACCOMPLISH  THIS  ARE  WRITTEN  INTO 
THE  SUBROUTINE  BUT  ARE  'COMMENTED  OUT'. 

'rVrVoVycVcVr<Wcy?yryoViW-yr*V«VVoWry?Vr;V}V*      MAIN    PROGRAM      4eMcfeMrlek&Mcirtricirtc$riefcMrMc4tieipli 


character  infile*12,  mode*l,  ylabl*13,  dsorce*l,  yscal*3 
real  mh(101),  ph(101),  thetav(lOl),  c(0:128),  b(0:128) 

C   PROMPT  USER  FOR  MODE:  BATCH  OR  TEST. 

write(*,1115) 

read(* , 1117)  mode 

if  ( (mode.  eq.  '  Y'  ).  or.  (mode.  eq.  '  y'  )  )  then 
mode  =  ' Y ' 
write(*,1118) 
read(*,1119)  infile 

else 
infile  =  'DIGFREQ. IN' 

end  if 

C  UNIT=1  DEFINED  AS  INPUT  FILE.   UNITS=2,3  DEFINED  AS  OUTPUT  FILES. 

open(unit=l , f ile=inf ile,status=' old' , iostat=ierr ,err=999) 
open(unit=2,file=' DIGFREQ. OUT' ) 
open(unit=3,file=' DIGFREQ.  DAT' ) 

C  READ  INPUT  PARAMETERS  AND  CONDUCT  ERROR  CHECKS. 

read( 1,1000)  numsys 

numplts  =  numsys*2 

write(3,2000)  numplts 

if  ( (numsys.  It.  1).  or.  (numsys.  gt.  3))  then 
write(*, 1122)  numsys 
stop  'Error,  numsys  must  be  in  the  range:  1  <=  numsys  <=  3. 

endif 
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do  10  nsys=l,  numsys 
data  mh/101*0. 0/,  ph/101*0. 0/,  thetav/101*0. 0/ 
data  b/129->0.  0/,  c/129*0.  0/ 

read(l,1001)  L,  N,  dsorce,  yscal 
read( 1,1002)  dlthta,  thetaO,  nurapts 

if((L.  lt.0).or.  (L.  gt.  128))  then 

write(,v ,  1124)  nsys,  L 

stop  'Error,  L  must  be  in  the  range:  0  <=  L  <=  128.  ' 
elseif((N.  It.  0).or.  (N.  gt.  128))  then 

write(*,1125)  nsys,  N 

stop  'Error,  N  must  be  in  the  range:  0  <=  N  <=  128. ' 
endif 

if((dsorce.  eq.  'F' ).  or.  (dsorce.  eq.  '  f '  ))  then 

dsorce  =  'F1 
elseif(( dsorce.  eq.  '  &'  ).  or.  (dsorce.  eq.  '  s'  ))  then 
dsorce  =  'S' 

else 
write(*,1018)  decree 

stop  'The  allowed  values  for  dsorce  are:  ''S1'  or  ''F''.' 
endif 

if((numpts. It. 1). or. (numpts. gt. 101))  then 
write(* , 1127)  nsys,  numpts 

stop  'Error,  numpts  must  be  in  the  range:  1  <=  numpts  <=  101. 
endif 

if  ( (yscal.  eq.  '  STD' ).  or.  (yscal.  eq.  '  std' ) )  then 

yscal  =  'STD' 

ylabl  =  '   MAGNITUDE   ' 
e Is eif ((yscal.  eq.  'LOG'  ).  or.  (yscal.  eq.  '  log'  ) )  then 

yscal  =  'LOG' 

ylabl  =  ' MAGNITUDE (db)1 
else 

write(-'-,1128)  yscal 

stop  'Error,  yscal  must  be  the  string:  ''LOG1'  or  ''STD1'. 
endif 

C  FOR  dsorce  =  'F'  READ  THE  COEFFICIENTS  b( )  AND  c( )  FROM  THE  INPUT 
C  FILE.   FOR  dsorce  =  'S'  CALL  coeff  TO  GENERATE  THE  COEFFICIENTS. 

if (dsorce. eq. 'F* )  then 

read(l,1003)  (b(k),k=0,L) 
read(l,1003)  (c(k),k=0,N) 

else 
call  coef f(L,N,nsys ,b,c) 
endif 

C  WRITE  INPUT  DATA  INTO  THE  OUTPUT  FILE:  DIGFREQ.  OUT. 

write(2,1008)  nsys 
write(2,1010)  infile 
write(2,1110)  L 
write(2,llll)  N 


S3 


write 
write 
write 
write 
write 
write 
write 
write 
write 
write 
write 


(2,1019) 
(2,1112) 
(2,1113) 
(2,1114) 
(2,1004) 
(2,1005) 
(2,1006) 
(2,1005) 
(2,1009) 
(2,1126) 
(2,1007) 


dsorce 

numpts , 

thetaO 

dlthta 


yscal 


(b(k),k=0,L) 

(c(k),k=0,N) 

nsys 

ylabl 


C  FOR  TEST  MODE  ECHO  ALL  INPUTS  ONTO  MONITOR  (UNIT  =  *). 


if(mode.  eq.  '  Y' 
write('v,1120 

write(*,1110 
write(*,llll 
write(--,1019 
write(»,1112 
write(*,1113 
write(*,1114 
write(-,1004 
write(»,1005 
write(--,1006 
write(*,1005 
write(-,1123 


then 

nsys 

L 

N 

dsorce 

numpts , 

thetaO 

dlthta 


inf ile 


yscal 


(b(k),k=0,L) 


(c(k),k=0,N) 
nsys 

pause  'END  OF  RUN  STRIKE  <CR>  WHEN  READY  TO  CONTINUE, 
endif 


C   CALL  dfresp  TO  COMPUTE  THE  FREQUENCY  RESPONSE. 

call  dfresp(b,c,mh,ph,L,N, thetaO , dltht a, thetav, numpts , yscal) 

C  WRITE  RESULTS  INTO  OUTPUT  FILE:  DIGFREQ.DAT. 

write(3,2001)  numpts 
write(3,*)  'MAGNITUDE  RESPONSE' 
write(3,»)  'THETA  (RADS)' 
write(3,2003)  ylabl 
do  55  np=l,  numpts 

write(3,2010)  thetav(np),  mh(np) 

55  continue 

write(3,2001)  numpts 
write(3,*)  'PHASE  RESPONSE' 
write(3,*)  'THETA  (RADS)' 
write(3,2003)  '  PHASE  (DEG)  ' 
do  56  np=l,  numpts 

write(3,2010)  thetav(np),  ph(np) 

56  continue 

C  WRITE  RESULTS  INTO  OUTPUT  FILE:  DIGFREQ.  OUT. 
do  150  np=l,  numpts 
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write(2,1013)  thetav(np),  rah(np),  ph(np) 
150     continue 
write(2,1123)  nsys 
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continue 


write(*,1121) 

999  close(unit=l) 
close(unit=2) 
close(unit=3) 

if(ierr.  gt.  0)  then 
write(*, 1116)  ierr 
endif 

C***********   INPUT  FORMAT  **************** 

1000  format(il) 

1001  format( i3,tll, i33t21 ,al ,t31,a3) 

1002  format(2fl0. 0,i3) 

1003  format(6fl0. 0) 


C* 


r  »\-  «s V  * V  ».*r  •£?  iV  -.V  -V  ». V  ■»'-  ^V  »fr  -.V  iV  iV  «»V  ■! 


5*r  "We  ■sWoWc'iV  ftfafe  Vr%V  -sV  -V  Vc  Vc  Vr  >V  Vc 


1004 
1005 
1006 

1007 
1008 
1009 
1010 
1013 
1018 
1019 
1110 
1111 
1112 

1113 
1114 
1115 

1116 

1117 
1118 


1119 
1120 

1121 

1122 


format( t4,' THE  NUMERATOR  COEFFICIENTS  b(0) ,b( 1).  .  .  b(L)  ARE:  ',/) 

format(6(2X,ell. 4),//) 
0format(//,t4,'THE  DENOMINATOR  COEFFICIENTS  c( 0) ,c( 1) . . . c(N) ' , 
1'  ARE:  ' ,/) 

format(t6,' (RADIANS)' , t38 , ' (DEGREES) ' ,/) 

format(tl6,'  INPUT  DATA  FOR  SYSTEM  #   ,il,//) 

format(///,tl6,'  OUTPUT  DATA  FOR  SYSTEM  #  ',il,/) 

format(t4,' INPUT  DATA  SOURCEFILE:  ' ,al2) 

format(t4,3(el2. 6,4x)) 

format( lx, ' dsorce  =  ' ,al,2x, 'Error ,  illegal  value  for  dsorce.  ' ) 

format( t4 , ' dsorce  =  *,al) 

format(t4, 'DEGREE  OF  NUMERATOR  =  '5i3) 

format( t4, 'DEGREE  OF  DENOMINATOR  =  ',i3) 
0format(t4,'NU:'3ER  OF  FREQUENCY  POINTS  =  ',  i3  ,t39  ,'  MAGNITUDE'  , 
1'  OPTION  =  ' ,a3) 

format(t4,' STARTING  VALUE  OF  THETA  =  ',el2.6) 

format(t4,' INCREMENT  OF  THETA  =  ',el2.6,/) 
0format(lx,'D0  YOU  WISH  TO  RUN  THIS  PROGRAM  IN  TEST', 
1'  MODE  ?    (Y/N)  <CR>  :  ' ,  ,) 

0format(///,lx, 'ERROR  OPENING  INPUT  FILE,  PROGRAM  TERMINATED.', 
l//,lx, 'ERROR  CODE: ' ,i4, /////) 

format(al) 
Cformat(/////,lx,'TYPE  THE  NAME  OF  YOUR  DATA  FILE  FOLLOWED', 
1'  BY  <CR>.  '  ,/,'  IF  YOU  DESIRE  TO  MAKE  A  TEST  RUN  USING  THE' , 
2'  SAMPLE  DATA  ALREADY  STORED',/,1  IN  THE  FILE:  DIGFREQ. TST' , 
3  '       TYPE:  DIGFREQ.  TST  <CR>'  ,/,'  FILENAME:  ',  ,) 

forraat( al2) 

format(////,t4,  'SYSTEM  #'  ,il,'     INPUT  DATA  SOURCEFILE:  \al2) 
Of ormat(//,t4, 'TABULAR  OUTPUT  DATA  IS  STORED  IN  FILE:  DIGFREQ. OUT' , 
l/,t4, 'PLOTTING  DATA  IS  STORED  IN  FILE:  DIGFREQ.DAT.  ') 

f ormat( ////// ,t2, 'The  value  of  numsys  is:  ',il,'.') 


S5 


1123  fonnatC/, lx,13( '-*),'   END  OF  RUN,  SYSTEM  #  ' , il,2x, 13( ' -' ) ,//) 

1124  Oformat(////// ,t2, 'The  degree(L)  of  the  numerator  for  system  ', 
l'#  ' ,  il, '  is  :   L  =  '  ,i3,  ' .  '  ) 

1125  0format(//////,t2, 'The  degree(N)  of  the  denominator  for  system' 
1,'  #  '  ,11,'  is   :  N  =  \i3,\  *) 

1126  format(///, t8,'THETA' ,t21 ,al3 ,t40, ' PHASE' ) 

1127  format(////// ,t2 , 'The  value  of  numpts  for  system  ',il,'  is:  '  ,i3) 

1128  format(/////, t2, 'The  value  of  yscal  is:  ',a3,'.') 

2000  format( il) 

2001  formate i3) 
2003   format(al3) 

2010   formateel2. 6,2x,el2.  6) 

end 


SUBROUTINE:  dfresp 


C  PURPOSE:   THIS  SUBROUTINE  COMPUTES  THE  FREQUENCY  RESPONSE  OF 

C  THE  SYSTEM.   ALL  FREQUENCY  CALCULATIONS  ARE  IN  RADIANS, 

C  HOWEVER  THE  OUTPUT  IS  CONVERTED  TO  DEGREES. 

C  THE  OUTPUT  FORMAT  FOR  EACH  FREQUENCY  INCREMENT  IS: 

C  MAGNITUDE(M)     PHASEeP)      AS  IN:   M*EXP(J*P). 


subroutine  dfresp(b,c,mh,ph,L,N, thetaO, dltht a, thetav, numpts , yscal) 

real  mh(numpts),  phenumpts),  thetav(numpts) ,  imz,  rez 

real  b(0:  L) ,  c(0:N) 

character  yscal"*3 

complex  z,  den,  num,  h,  ci 

C   DEFINE  CONSTANTS. 

pi  =  4.  0*atan(l.  0) 

ci  =  (1.  0,0.  0) 

C  ITERATE  FROM  thetaO,  IN  INCREMENTS  OF  dlthta. 

do  100  np=l,  numpts 
num  =  ci'vb(0) 
den  =  ci-''c(0) 

thetav(np)  =  thetaO  +  (np-l)*dlthta 
rez  =  cos(thetav(np)) 
imz  =  sin(thetav(np)) 
z  =  cmplx(rez,imz) 

C  CALCULATE  NUMERATOR  FOR  GIVEN  VALUE  OF  THETA,  IF  L  >  0. 

if(L.  gt.  0)  then 
do  50  k=l,  L 

num  =  z*num  +  ci''-'b(k) 
50       continue 

end  if 
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C  CALCULATE  DENOMINATOR  FOR  GIVEN  VALUE  OF  THETA,  IF  N  >  0. 

if(N.  gt.  0)  then 
do  70  k=l,  N 

den  =  z-den  +  ci-c(k) 
70       continue 
end  if 
h  =  num/den 


CONVERT  COMPLEX  VALUE  V  INTO  MAGNITUDE(rah)  AND  PHASE(ph)  TERMS. 
IF  yscal  =  'LOG'  THEN  CONVERT  MAGNITUDE  TO  DECIBELS  (db). 


C 

c 

C  DIVIDE  BY  ZERO  AVOIDED  BY  'if'  STATEMENTS. 

mh(np)  =  cabs(h) 
if(yscal.eq.  'LOG')  then 

if(mh(np).  gt.  0.  00001)  then 
mh(np)  =  20.0*logl0(mh(np)) 
else 

tnh(np)  =  -100.  0 
endif 
endif 

if(abs(real(h)).  It.  1.0e-15)  then 

if(abs(aimag(h)).  le.  1.  0e-15)  ph(np)=0.  0 
if(aimag(h).  gt.  1.  Oe-15)  ph(np)=90.  0 
if(aimag(h).  It.  -1.  0e-15)  ph(np)=-90.  0 

else 
ph(np)  =  (180.  0/pi)*atan2(aimag(h),real(h)) 
endif 

100   continue 


return 
end 


SUBROUTINE:  coeff 


C  PURPOSE:   THIS  SUBROUTINE  ALLOWS  THE  USER  TO  GENERATE  THE 
C  NUMERATOR  AND  DENOMINATOR  COEFFICIENTS  THAT  DESCRIBE 

C  EACH  SYSTEM  TO  BE  ANALYZED.   IF  dsorce  =  'S'  THEN 

C  THE  MAIN  PROGRAM  WILL  CALL  THIS  SUBROUTINE. 


subroutine  coef f (L,N5nsys ,b,c) 
real  b(0:L),  c(0:N) 

pi  =  4.  0*atan(l.  0) 
el  =  L 
en  =  N 

C  DEVELOP  THE  EQUATIONS  TO  GENERATE  VALUES  FOR  THE  ARRAYS  b( )  AND  c( ) 
C   IN  THIS  SPACE.   THE  STATEMENTS  TYPED  IN  MUST  FOLLOW  STANDARD 
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C  FORTRAN  77  RULES  AND  MAY  USE  FORTRAN  77  INTRINSIC  FUNCTIONS  SUCH  AS: 

C   SIN(),  COSQ,  ABS()...   AN  EXAMPLE  IS  SHOWN  BELOW.   NOTE  THAT  THE 

C  VALUE  nsys  CAN  BE  USED  TO  DISTINGUISH  BETWEEN  SYSTEMS  IF  MORE  THAN 

C  ONE  SYSTEM  (nurasvs  >  1)  IS  TO  BE  ANALYZED. 

C 

C  ***  EXAMPLE  *** 

C 

C      if(nsys.  eq.  1)  then 

C        do  2  i=0,  L 

C  b(i)  =  cos(i*pi/(2.0*el)) 

C  2      continue 

C        do  3  i=0,  N 

C         c(i)  =  cos(2.0*i*pi/(3.0*en)) 

C  3      continue 

C      endif 


return 
end 


B.     MAGNITUDE,  PHASE  OF  SECOND  ORDER  IIR  NOTCH  FILTER:  P  =  0 

INPUT  DATA  FOR  SYSTEM  #  1 

INPUT  DATA  SOURCEFILE:  digfreq.  tst 

DEGREE  OF  NUMERATOR  =   2 

DEGREE  OF  DENOMINATOR  =   2 

dsorce  =  F 

NUMBER  OF  FREQUENCY  POINTS  =  100   MAGNITUDE  OPTION  =  STD 

STARTING  VALUE  OF  THETA  =   . 000000E+00 

INCREMENT  OF  THETA  =   . 280600E-01 

THE  NUMERATOR  COEFFICIENTS  b( 0) ,b( 1) . . . b(L)  ARE: 

. 9900E+00    -.  3299E+00    . 9900E+00 

THE  DENOMINATOR  COEFFICIENTS  c(0) ,c( 1). . . c(N)  ARE: 
. 1000E+01    -.3299E+00     . 9801E+00 


OIJTF'JI  DATA  FOR  SYSTEM  #  1 


THETA 
(RADIANS) 


MAGNITUDE 


PHASE 
(DEGREES) 


OOOOOOE+OO 
280600E-01 
561200E-01 
841800E-01 
112240E+00 
140300E+00 
168360E+00 
196420E+00 
224480E+C0 
252540E+00 
280600E+00 
308660E+00 
336720E+00 
364780E+00 
392340E+00 
420900E+00 
448960E+00 
477020E+00 
505080E+00 
533140E+00 
561200E+00 
589260E+00 
617320E+00 
645380E+00 
673440E+00 
7O1500E+O0 
729560E+00 
757620E+00 
785680E+00 
813740E+00 
841800E+00 
869860E+00 
897920E+00 
925980E+00 
954040E+00 
98210CE+00 
101016E+01 
103822E+01 
106628E+01 
109434E+01 
112240E+01 
, 115046E+01 
, 117852E+01 
, 120658E+01 
, 123464E+01 
. 126270E+01 
. 129076E+C1 


999939E+00 
999939E+00 
999939E+00 
999939E+00 
999938E+00 
999938E+00 
999937E+00 
999936E+00 
999935E+00 
999934E+00 
999933E+00 
999931E+00 
999930E+00 
999928E+00 
999926E+00 
999923E+00 
999920E+00 
999917E+00 
999914E+00 
999910E+00 
999906E+00 
999902E+00 
999896E+00 
999891E+00 
999884E+00 
999876E+00 
999868E+00 
999858E+00 
999847E+00 
999835E+00 
999820E+00 
999803E+00 
999783E+00 
999759E+00 
999731E+00 
999697E+00 
999655E+00 
999603E+00 
999537E+00 
. 999453E+00 
, 999343E+00 
, 999193E+00 
, 998985E+00 
. 998681E+00 
. 998213E+00 
. 997440E+00 
. 996024E+00 


OOOOOOE+00 
193947E-01 
388289E-01 
583436E-01 
779784E-01 
977754E-01 
117779E+00 
138030E+00 
158581E+00 
179474E+00 
200769E+00 
222516E+00 
244780E+00 
267620E+00 
291112E+00 
315332E+00 
340366E+00 
366300E+00 
393251E+00 
421327E+00 
450663E+00 
481405E+00 
513723E+00 
547804E+00 
583872E+00 
622182E+00 
663033E+00 
706764E+00 
753788E+00 
804587E+00 
859743E+00 
919964E+00 
986100E+00 
105921E+01 
114061E+01 
123199E+01 
133548E+01 
145389E+01 
159095E+01 
, 175176E+01 
, 194344E+01 
, 217626E+01 
. 246560E+01 
, 283559E+01 
, 332633E+01 
. 400969E+01 
. 502855E+01 


S9 


131882E+01 
134688E+01 
137494E+01 
140300E+01 
143106E+01 
145912E+01 
148718E+01 
151524E+01 
154330E+01 
157136E+01 
159942E+01 
162748E+01 
165554E+01 
168360E+01 
171166E+01 
173972E+01 
176778E+01 
179584E+01 
182390E+01 
185196E+01 
188002E+01 
190808E+01 
193614E+01 
196420E+01 
199226E+01 
202032E+01 
204838E+01 
207644E+01 
210450E+01 
213256E+01 
216062E+01 
218868E+01 
221674E+01 
224480E+01 
227286E+01 
230092E+01 
232898E+01 
235704E+01 
238510E+01 
241316E+01 
244122E+01 
246928E+01 
249734E+01 
252540E+01 
255346E+01 
258152E+01 
260958E+01 
263764E+01 
266570E+01 
269376E+01 
272182E+01 
274988E+01 
277794E+01 


. 992994E+00 
. 984530E+00 
. 942857E+00 
. 381903E-01 
. 939926E+00 
. 984111E+00 
. 992867E+00 
.995971E+00 
.997413E+00 
. 998198E+00 
. 998672E+00 
. 998979E+00 
.'999190E+00 
. 999341E+00 
. 999452E+00 
. 999537E+00 
. 999603E+00 
. 999655E+00 
.999698E+00 
. 999732E+00 
. 999761E+00 
. 999785E+00 
. 999805E+00 
. 999822E+00 
. 999837E+00 
. 999850E+00 
. 999861E+00 
. 999871E+00 
. 999880E+00 
. 999888E+00 
. 999894E+00 
. 999901E+00 
. 999906E+00 
. 999911E+00 
. 999915E+00 
.999919E+00 
. 999923E+00 
. 999926E+00 
. 999929E+00 
. 999932E+00 
. 999934E+00 
. 999937E+00 
. 999939E+00 
. 999941E+00 
. 999942E+00 
. 999944E+00 
. 999945E+00 
. 999947E+00 
.  999948E+00 
.  999949E+00 
. 999950E+00 
.999951E+00 
. 999952E+00 


671241E+01 
100259E+02 
194059E+02 
877665E+02 
200014E+02 
102595E+02 
687161E+01 
516115E+01 
413001E+01 
344006E+01 
294546E+01 
257310E+01 
228226E+01 
204852E+01 
185629E+01 
169521E+01 
155808E+01 
143978E+01 
133653E+01 
124550E+01 
116455E+01 
109198E+01 
102647E+01 
966960E+00 
912579E+00 
862632E+00 
816531E+00 
773792E+00 
734015E+00 
696843E+00 
661987E+00 
629193E+00 
598246E+00 
568950E+00 
541144E+00 
514683E+00 
489434E+00 
465285E+00 
442146E+00 
419918E+00 
398521E+00 
377889E+00 
357941E+00 
338639E+00 
319915E+00 
301724E+00 
284024E+00 
266771E+00 
249929E+00 
233459E+00 
217332E+00 
201519E+00 
185987E+00 


END  OF  RUN,  SYSTEM  #  1 
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C.     MAGNITUDE,  PHASE  OF  THE  AUGMENTED  IIR  NOTCH  FILTER:  P=  1 

INPUT  DATA  FOR  SYSTEM  #  1 

INPUT  DATA  SOURCEFILE:  digf req.  tst 

DEGREE  OF  NUMERATOR  =   3 

DEGREE  OF  DENOMINATOR  =   3 

dsorce  =  F 

NUMBER  OF  FREQUENCY  POINTS  =  100   MAGNITUDE  OPTION  =  STD 

STARTING  VALUE  OF  THETA  =   . 000000E+00 

INCREMENT  OF  TKETA  =   . 280600E-01 

THE  NUMERATOR  CO:  TICIENTS  b(0)  ,b(  1).  .  .  b(L)  ARE: 

.9900E+00    -. 3299E-02    . 8811E+00    . 3266E+00 

THE  DENOMINATOR  COEFFICIENTS  c(0) ,c( 1). . . c(N)  ARE: 
. 1000E+01     .0000E+00     .8712E+00     . 3234E+00 

OUTPUT  DATA  FOR  SYSTEM  #  1 


THETA 
(RADIANS) 

000000E+00 
280600E-01 
561200E-01 
841800E-01 
112240E+00 
140300E+00 
168360E+00 
196420E+00 
224480E+00 
252540E+00 
280600E+00 
308660E+00 
336720E+00 
364780E+00 
392840E+00 
420900E+00 
448960E+00 
477020E+00 


MAGNITUDE 


999940E+00 
999940E+00 
999939E+00 
999939E+00 
999939E+00 
999938E+00 
999937E+00 
999936E+00 
999935E+00 
999934E+00 
999933E+00 
999932E+00 
999930E+00 
999928E+00 
999926E+00 
999923E+00 
999921E+00 
999918E+00 


PHASE 
(DEGREES) 

000000E+00 
193950E-01 
388293E-01 
583432E-01 
779798E-01 
977766E-01 
117780E+00 
138031E+00 
158579E+00 
179475E+00 
200768E+00 
222515E+00 
244780E+00 
267618E+00 
291112E+00 
315332E+00 
340364E+00 
366303E+00 


91 


505080E+00 
533140E+00 
561200E+00 
589260E+00 
617320E+00 
645380E+00 
673440E+00 
701500E+00 
729560E+00 
757620E+00 
785680E+00 
813740E+00 
841800E+00 
869860E+00 
897920E+00 
925980E+00 
954040E+00 
982100E+00 
1C1016E+01 
103822E+01 
106628E+01 
109434E+01 
112240E+01 
115046E+01 
117852E+01 
120658E+01 
123464E+01 
126270E+01 
129076E+01 
131882E+01 
134688E+01 
137494E+01 
140300E+01 
143106E+01 
145912E+01 
14871SE+01 
151524E+01 
154330E+01 
157136E+01 
159942E+01 
162748E+01 
165554E+01 
168360E+0I 
171166E+01 
173972E+01 
176778E+01 
179584E+01 
182390E+01 
185196E+01 
188002E+01 
190808E+01 
193614E+01 
196420E+01 
199226E+01 
202032E+01 
204838E+01 


999914E+00 
999911E+00 
999906E+00 
999902E+00 
999897E+00 
999891E+00 
999884E+00 
999877E+00 
999868E+00 
999859E+00 
999848E+00 
999835E+00 
999820E+00 
999803E+00 
999783E+00 
999759E+00 
999731E+00 
999697E+00 
999655E+00 
999603E+00 
999538E+00 
999454E+00 
999343E+00 
999194E+00 
998985E+00 
998681E+00 
998214E+00 
997441E+00 
996025E+00 
992995E+00 
984532E+00 
942860E+00 
382032E-01 
939922E+00 
984110E+00 
992866E+00 
995970E+00 
997413E+00 
998197E+00 
998671E+00 
998979E+00 
999189E+00 
999340E+00 
999452E+00 
999537E+00 
999603E+00 
999655E+00 
999697E+00 
999732E+00 
999761E+00 
999785E+00 
999805E+00 
, 999822E+00 
, 999837E+00 
, 999850E+00 
, 999861E+00 


393250E+00 

421327E-J-00 

450665E+00 

481406E+00 

513721E+00 

547805E+00 

583875E+00 

622183E+00 

663032E+00 

706764E+00 

753788E+00 

804590E+00 

859746E+00 

919964E+00 

986101E+00 

105921E+01 

114061E+01 

123199E+01 

133548E+01 

145389E+01 

159095E+01 

175176E+01 

194344E+01 

217625E+01 

246560E+01 

283559E+01 

332632E+01 

400968E+01 

502853E+01 

671239E+01 

100259E+02 

194058E+02 

877594E+02 

, 200013E+02 

,  102595E+02 

, 687158E+01 

, 516112E+01 

, 412999E+01 

, 344004E+01 

, 294544E+01 

. 257309E+01 

, 228225E+01 

, 204850E+01 

, 185627E+01 

, 169520E+01 

,  155807E+01 

, 143977E+01 

.  133652E+01 

. 124549E+01 

, 116454E+01 

, 109197E+01 

, 102646E+01 

. 966947E+00 

. 912566E+00 

. 862621E+00 

. 816522E+00 
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207644E+01 
210450E+01 
213256E+01 
216062E+01 
218868E+01 
221674E+01 
224480E+01 
227286E+01 
230092E+01 
232898E+01 
235704E+01 
238510E+01 
241316E+01 
244122E+01 
246928E+01 
249734E+01 
252540E+01 
255346E+01 
258152E+01 
260958E+01 
263764E+01 
266570E+01 
269376E+01 
272182E+01 
274988E+01 
277794E+01 


999871E+00 
999880E+00 
999887E+00 
999894E+00 
999900E+00 
999906E+00 
999911E+00 
999915E+00 
999919E+00 
999923E+00 
999926E+00 
999929E+00 
999932E+00 
999934E+00 
999936E+00 
999939E+00 
999940E+00 
999942E+00 
999944E+00 
999945E+00 
999947E+00 
999948E+00 
999949E+00 
999950E+00 
999951E+00 
999952E+00 


. 773783E+00 
. 734003E+00 
. 696833E+00 
. 661978E+00 
. 629187E+00 
.598235E+00 
.568938E+00 
.541135E+00 
.514669E+00 
.489423E+00 
.465277E+00 
.442138E+00 
.419907E+00 
. 398514E+00 
. 377878E+00 
.357937E+00 
. 338629E+00 
. 319905E+00 
. 301716E+00 
. 284016E+00 
. 266763E+00 
. 249920E+00 
. 233452E+00 
.217326E+00 
.201509E+00 
. 185980E+00 


END  OF  RUN,  SYSTEM  #  1 


D.     MAGNITUDE.  PHASE  OF  THE  AUGMENTED  IIR  NOTCH  FILTER:  P  =  2 


INPUT  DATA  FOR  SYSTEM  #  1 


INPUT  DATA  SOURCr.FILE:  digfreq.  tst 

DEGREE  OF  NUMERAL :R  =   4 

DEGREE  OF  DENOMINATOR  =   4 

dsorce  =  F 

NUMEER  OF  FREQUENCY  POINTS  =  100   MAGNITUDE  OPTION  =  STD 

STARTING  VALUE  OF  THETA  =  . 000000E+00 

INCREMENT  OF  THETA  =  . 280600E-01 


THE  NUMERATOR  COEFFICIENTS  b(0) ,b( 1). . . b(L)  ARE: 
. 9900E+00    -.3299E-02    . 1861E-01     . 6141E+00 


•.  8625E+00 


THE  DENOMINATOR  COEFFICIENTS  c(0) ,c( 1). . . c(N)  ARE: 
. 1000E+01    . 0000E+00    . 0000E+00    . 6108E+00 


8539E+00 
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OUTPUT  DATA  FOR  SYSTEM  #  1 


THETA 
(RADIANS) 


MAGNITUDE 


PHASE 
(DEGREES) 


OOOOOOE+OO 
280600E-01 
561200E-01 
841800E-01 
112240E+00 
140300E+00 
168360E+00 
196420E+00 
224480E+00 
252540E+00 
280600E+00 
308660E+00 
336720E+00 
364780E+00 
392840E+00 
420900E+00 
448960E+00 
477020E+00 
505080E+00 
533140E+00 
561200E+00 
589260E+00 
617320E+00 
645380E+00 
673440E+00 
701500E+00 
729560E+00 
757620E+00 
785680E+00 
813740E+00 
841800E+00 
869860E+00 
897920E+00 
925980E+00 
954040E+00 
982100E+00 
101016E+01 
103822E+01 
106628E+01 
109434E+01 
112240E+01 
115046E+01 
117852E+01 
120658E+01 
123464E+01 


999940E+00 
999940E+00 
999940E+00 
999940E+00 
999939E+00 
999938E+00 
999938E+00 
999937E+00 
999936E+00 
999934E+00 
999933E+00 
999932E+00 
999930E+00 
999928E+00 
999926E+00 
999923E+00 
999920E+00 
999917E+00 
999914E+00 
999910E+00 
999906E+00 
999901E+00 
999896E+00 
999891E+00 
999834E+00 
999876E+00 
999868E+00 
999858E+00 
999847E+00 
999835E+00 
999820E+00 
999803E+00 
999783E+00 
999759E+00 
999731E+00 
999697E+00 
999655E+00 
999603E+00 
999537E+00 
999453E+00 
999343E+00 
999193E+00 
998985E+00 
998681E+00 
998213E+00 


000000E+00 
194012E-01 
388423E-01 
583613E-01 
780022E-01 
978035E-01 
117802E+00 
138056E+00 
158606E+00 
179500E+00 
200794E+00 
222545E+00 
244805E+00 
267638E+00 
291136E+00 
315351E+00 
340386E+00 
366321E+00 
393270E+00 
421345E+00 
450679E+00 
481422E+00 
513734E+00 
547815E+00 
583885E+00 
622197E+00 
663044E+00 
706773E+00 
753796E+00 
804595E+00 
859751E+00 
919967E+00 
986102E+00 
105921E+01 
114062E+01 
123199E+01 
133548E+01 
145388E+01 
159094E+01 
175175E+01 
194343E+01 
217625E+01 
246558E+01 
283557E+01 
332630E+01 
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126270E+01 
129076E+01 
131882E+01 
134688E+01 
137494E+01 
140300E+01 
143106E+01 
145912E+01 
148718E+01 
151524E+01 
154330E+01 
157136E+01 
159942E+01 
162748E+01 
165554E+01 
168360E+01 
171166E+01 
173972E+01 
176778E+01 
179584E+01 
182390E+01 
185196E+01 
188002E+01 
190808E+01 
193614E+01 
196420E+01 
199226E+01 
202032E+01 
204838E+01 
207644E+01 
210450E+01 
213256E+01 
216062E+01 
218868E+01 
221674E+01 
224480E+01 
227286E+01 
230092E+01 
232898E+01 
235704E+01 
238510E+01 
241316E+01 
244122E+01 
246928E+01 
249734E+01 
252540E+01 
255346E+01 
258152E+01 
260958E+01 
263764E+01 
266570E+01 
269376E+01 
272182E+01 
274988E+01 
277794E+01 


997440E+00 
996024E+00 
992994E+00 
984530E+00 
942856E+00 
3S1887E-01 
939924E+00 
984111E+00 
992867E+00 
995970E+00 
997413E+00 
998198E+00 
998671E+00 
998979E+00 
999189E+00 
999340E+00 
999452E+00 
999537E+00 
999603E+00 
999655E+00 
999697E+00 
999732E+00 
999761E+00 
999785E+00 
999805E+00 
999822E+00 
999837E+00 
999850E+00 
999861E+00 
999871E+00 
999880E+00 
999887E+00 
999894E+00 
999900E+00 
999906E+00 
999911E+00 
999915E+00 
999919E+00 
999923E+00 
999926E+00 
999929E+00 
999932E+00 
999934E+00 
999936E+00 
999938E+00 
999940E+00 
999942E+00 
999944E+00 
999945E+00 
999947E+00 
999948E+00 
999949E+00 
999950E+00 
999951E+00 
999952E+00 


400966E+01 
502849E+01 
671234E+01 
100258E+02 
194057E+02 
877335E+02 
200011E+02 
102594E+02 
687151E+01 
516108E+01 
412995E+01 
344001E+01 
294542E+01 
257306E+01 
228223E+01 
204849E+01 
185626E+01 
169518E+01 
155806E+01 
143976E+01 
133650E+01 
124548E+01 
116454E+01 
109196E+01 
102646E+01 
966945E+00 
912560E+00 
862616E+00 
816516E+00 
773779E+00 
734001E+00 
696832E+00 
661977E+00 
629182E+00 
598231E+00 
568938E+00 
541134E+00 
514669E+00 
489424E+00 
465279E+00 
442136E+00 
419910E+00 
398511E+00 
377878E+00 
357935E+00 
338632E+00 
319908E+00 
301717E+00 
284017E+00 
266763E+00 
249920E+00 
233450E+00 
217325E+00 
, 201511E+00 
, 185983E+00 
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END  OF  RUN,    SYSTEM  #   1 


E.     MAGNITUDE,  PHASE  OF  THE  AUGMENTED  IIR  NOTCH  FILTER:  P=  3 


INPUT  DATA  FOR  SYSTEM  #   1 


INPUT  DATA  SOURCEFILE:  digfreq. tst 

DEGREE  OF  NUMERATOR  =   5 

DEGREE  OF  DENOMINATOR  =   5 

dsorce  =  F 

NUMBER  OF  FREQUENCY  POINTS  =  100   MAGNITUDE  OPTION  =  STD 

STARTING  VALUE  OF  THETA  =   . 000000E+00 

INCREMENT  OF  THETA  =   . 280600E-01 


THE  NUMERATOR  COEFFICIENTS  b(0) ,b( 1). . . b(L)  ARE: 

. 9900E+00    -.3299E-02     . 1861E-01     . 9374E-02 
-.6047E+00 


-.6610E+00 


THE  DENOMINATOR  COEFFICIENTS  c( 0) , c( 1). . . c(N)  ARE: 

.  1000E+01    .  0000E+00    . 0000E+00    . 0000E+00    -. 6524E+00 
-.5987E+00 


OUTPUT  DATA  FOR  SYSTEM  #  1 


THETA 

MAGNITUDE 

PHASE 

(RADIANS) 

(DEGREES) 

000000E+00 

. 999937E+00 

. 000000E+00 

280600E-01 

. 999938E+00 

-. 194494E-01 

561200E-01 

. 999938E+00 

-. 388970E-01 

841800E-01 

. 999938E+00 

-.584018E-01 

112240E+00 

. 999938E+00 

-. 780285E-01 

140300E+00 

. 999938E+00 

-.978157E-01 

168360E+00 

. 999937E+00 

-. 117817E+00 

196420E+00 

. 999936E+00 

-. 138065E+00 

224480E+00 

. 999935E+00 

-. 158610E+00 

252540E+00 

. 999934E+00 

-. 179497E+00 
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280600E+00 

308660E+00 

33&720E+00 

364780E+00 

392840E+00 

420900E+00 

448960E+00 

477020E+00 

505080E+00 

533140E+00 

561200E+00 

589260E+00 

617320E+00 

645380E+00 

673440E+00 

701500E+00 

729560E+00 

757620E+00 

785680E+00 

813740E+00 

841800E+00 

869860E+00 

897920E+00 

925980E+00 

954040E+00 

98210CE+00 

101016E+01 

103822E+01 

106628E+0I 

109434E+01 

112240E+01 

115046E+01 

117852E+01 

120658E+01 

123464E+01 

126270E+01 

129076E+01 

131882E+01 

134688E+01 

137494E+01 

, 140300E+01 

, 143106E+01 

, 145912E+01 

, 148718E+01 

. 151524E+01 

, 154330E+01 

, 157136E+01 

. 159942E+01 

. 162748E+01 

, 165554E+01 

. 168360E+01 

. 171166E+01 

. 173972E+01 

. 176778E+01 

. 179584E+01 

. 182390E+01 


999933E+00 

999931E+00 

999929E+00 

999927E+00 

999925E+00 

999923E+00 

999920E+00 

999917E+00 

999914E+00 

999910E+00 

999906E+00 

999902E+00 

999896E+00 

999891E+00 

999884E+00 

999876E+00 

999868E+00 

999858E+00 

999847E+00 

999835E+00 

999820E+00 

999803E+00 

999783E+00 

999759E+00 

999731E+00 

999697E+00 

999655E+00 

999603E+00 

999537E+00 

999453E+00 

999343E+00 

999193E+00 

998985E+00 

998681E+00 

9982I3E+00 

997440E+00 

996024E+00 

, 992994E+00 

, 984530E+00 

, 942356E+00 

, 381933E-01 

, 939924E+00 

. 984110E+00 

. 992867E+00 

. 995970E+00 

. 997413E+00 

. 998198E+00 

. 998671E+00 

. 998979E+00 

. 999189E+00 

. 999340E+00 

. 999452E+00 

. 999537E+00 

. 999603E+00 

. 999655E+00 

. 999697E+00 


200787E+00 

222538E+00 

244798E+00 

267635E+00 

291127E+00 

315344E+00 

340375E+00 

366312E+00 

393263E+00 

421337E+00 

450671E+00 

481411E+00 

513727E+00 

547810E+00 

583876E+00 

622189E+00 

663036E+00 

706766E+00 

753794E+00 

804590E+00 

859744E+00 

919960E+00 

986096E+00 

105921E+0J. 

114061E+01 

123199E+01 

133547E+01 

145388E+01 

159094E+01 

175175E+01 

194343E+01 

217624E+01 

246558E+01 

283557E+01 

, 332630E+01 

, 400966E+01 

, 502848E+01 

, 671234E+01 

. 100258E+02 

. 194057E+02 

. 877273E+02 

, 200011E+02 

. 102594E+02 

. 687151E+01 

. 516108E+01 

.412995E+01 

. 344001E+01 

. 294542E+01 

. 257306E+01 

. 228223E+01 

. 204848E+01 

. 185626E+01 

. 169518E+01 

. 155805E+01 

. 143975E+01 

. 133650E+01 
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185196E+01 
188002E+01 
190808E+01 
193614E+01 
196420E+01 
199226E+01 
202032E+01 
204838E+01 
207644E+01 
210450E+01 
213256E+01 
216062E+01 
218868E+01 
221674E+01 
224480E+01 
227286E+01 
230092E+01 
232898E+01 
235704E+01 
238510E+01 
241316E+01 
244122E+01 
246928E+01 
249734E+01 
252540E+01 
255346E+01 
258152E+01 
260958E+01 
263764E+01 
266570E+01 
269376E+01 
272182E+01 
274988E+01 
277794E+01 


999732E+00 
999761E+00 
999785E+00 
999805E+00 
999822E+00 
999837E+00 
999850E+00 
999861E+00 
999871E+00 
999880E+00 
999887E+00 
999894E+00 
999900E+00 
999906E+00 
999911E+00 
999915E+00 
999919E+00 
999922E+00 
999926E+00 
999929E+00 
999932E+00 
999934E+00 
999936E+00 
999938E+00 
999940E+00 
999942E+00 
999943E+00 
999945E+00 
999946E+00 
999948E+00 
999949E+00 
999950E+00 
999951E+00 
999952E+00 


. 124548E+01 
. 116453E+01 
. 109196E+01 
. 102645E+01 
.966937E+00 
.912553E+00 
. 862608E+00 
. 816508E+00 
.  773770E+00 
. 733993E+00 
.696823E+00 
. 661968E+00 
.629177E+00 
.598225E+00 
.568932E+00 
.541125E+00 
.514664E+00 
.489416E+00 
.465272E+00 
.442129E+00 
.419903E+00 
. 398506E+00 
. 377872E+00 
. 357931E+00 
. 338631E+00 
. 319911E+00 
. 301721E+00 
. 284018E+00 
. 266771E+00 
. 249936E+00 
. 233465E+00 
. 217347E+00 
.201533E+00 
. 186002E+00 


END  OF  RUN,  SYSTEM  #  1 
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APPENDIX  B.     PROGRAM  AND  DATA  OF  DIFFERENCE  EQUATION 

A.  PROGRAM  FOR  DIFFERENCE  EQUATION  CALCULATION 

C  DIFFEQ.  FOR          VERSION:  2/03/88 

C 

C 

C  PURPOSE:   THIS  PROGRAM  COMPUTES  THE  ITERATIVE  SOLUTION  TO  A 

C  LINEAR,  TIME -INVARIANT  (LTI)  DIFFERENCE  EQUATION. 

C  THE  DIFFERENCE  EQUATION  MUST  BE  IN  THE  FORM: 

C  y(ns)  =  a(l)*y(ns-l)  +  ...  +  a(N)*y(ns-N)  + 

C  b(0)*x(ns)  +  b(l)*x(ns-l)  +  ...  +  b(L)*x(ns-L). 

C  THE  PROGRAM  CONSISTS  OF  A  MAIN  PROGRAM  AND  TOO 

C  SUBROUTINES.   SUBROUTINE  diffeq  IS  CALLED  BY  THE  MAIN 

C  PROGRAM  TO  ITERATIVELY  SOLVE  THE  DIFFERENCE  EQUATIONS  AND 

C  SUBROUTINE  xgen  ALLOWS  THE  USER  THE  OPTION  OF  GENERATING 

C  THE  INPUT  SEQUENCE  x( )  BY  WRITING  THE  APPROPRIATE 

C  EQUATIONS.   IF  THE  USER  ELECTS  TO  GENERATE  THE  SEQUENCE 

C  x()  BY  "JSING  xgen  THEN  THE  PROGRAM  MUST  BE  COMPILED 

C  AGAIN  BEFORE  EXECUTION.   THE  USER  HAS  THE  OPTION  OF 

C  SELECTING  ONE  OF  TWO  OPERATING  MODES:  BATCH  OR  TEST. 

C  IN  BATCH  MODE  THE  AMOUNT  OF  INTERFACE  WITH  THE  USER 

C  IS  MINIMIZED  AND  IT  IS  ASSUMED  THAT  THE  INPUT  DATA 

C  HAS  BEEN  STORED  IN  THE  DEFAULT  FILE  'DIFFEQ.  IN'.   IN 

C  TEST  MODE  THE  USER  IS  PROMPTED  FOR  THE  NAME  OF  THE 

C  INPUT  FILE  OR  HAS  THE  OPTION  OF  PERFORMING  A  TEST  RUN 

C  USING  THE  INPUT  DATA  STORED  IN  THE  FILE  ' DIFFEQ. TST* . 

C  IT  IS  RECOMMENDED  THAT  FIRST-TIME  USERS  SELECT  THE 

C  TEST  MODE  AND  MAKE  A  TRIAL  RUN  WITH  THE  PRESTORED 

C  INPUT  DATA.   THE  TEST  MODE  ECHOES  PORTIONS  OF  THE 

C  INPUT  DATA  ONTO  THE  MONITOR  TO  ALLOW  VERIFICATION  OF 

C  ITS  ACCURACY.   BOTH  BATCH  AND  TEST  MODES  ALLOW  THE 

C  USER  TO  SOLVE  UP  TO  FOUR  DIFFERENCE  EQUATIONS  IN  A 

C  SINGLE  PROGRAM  EXECUTION.   THE  OUTPUT  OF  THE  PROGRAM 

C  'DIFFEQ. FOR'  IS  STORED  IN  THE  ARRAY  y().   THE  OUTPUT  IS 

C  STORED  IN  TABULAR  FORM  IN  THE  OUTPUT  FILE  'DIFFEQ. FOR' 

C  AND  IN  A  FORM  SUITABLE  FOR  PLOTTING  IN  THE  FILE 

C  'DIFFEQ.DAT'. 

C 

C 

c 
c 

C  THIS  PROGRAM  ASSUMES  THAT  EACH  DIFFERENCE  EQUATION  IS  IN  THE 

C  CANONICAL  FORM: 

C 

C  y(ns)  =  a(l)*y(ns-l)  +  ...  +  a(N)*y(ns-N)  + 

C  b(0)*x(ns)  +  b(l)*x(ns-l)  +  ...  +  b(L)*x(ns-L) 

C 

C  L  =  A  NON-NEGATIVE  INTEGER,  THE  NUMBER  OF  INPUT  DELAYS. 

C  N  =  A  NON-NEGATIVE  INTEGER,  THE  NUMBER  OF  OUTPUT  DELAYS. 

C 
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a(l)...a(N)  =  REAL  COEFFICIENTS  OF  THE  OUTPUT  TERMS. 
b(0)...b(L)  =  REAL  COEFFICIENTS  OF  THE  INPUT  TERMS. 

THE  INPUT  PARAMETERS  SHOULD  BE  STORED  IN  A  FILE  NAMED 
'DIFFEQ. IN'.   ALL  OF  THE  READ  STATEMENTS  USED  BY  THIS  PROGRAM 
REQUIRE  FORMATTED  INPUT.   PARTICULAR  ATTENTION  SHOULD  BE  PAID 
TO  THE  FORMATS,  ESPECIALLY  THE  USE  OF  THE  DECIMAL  POINT  TO 
DENOTE  'REAL'  NUMBERS.   THE  INPUT  PARAMETERS  REQUIRED  BY  THE 
PROGRAM  ARE  LISTED  BELOW. 


NAME 

numsys 

L 

N 

nstop 

xsorce 

b(k) 

a(k) 

Y(k) 

x(ns) 

WHERE: 

numsys  = 


TYPE 


RANGE  (ARRAYS) 


INTEGER 

INTEGER 

INTEGER 

INTEGER 

CHARACTER 

REAL 

0,1, 

...  ,  L 

REAL 

1,2, 

.  .  .  ,  N 

REAL 

-N,  . 

..  ,  -1 

REAL 

0 

,1,  •• 

.  ,  nstop 

RESTRICTIONS 

1  <=  numsys  <=  4 

0  <=  L  <=  128 

0  <=  N  <=  128 

0  <=  nstop  <=  999 

'F'  OR  'S' 

0  <=  L  <=  128 

0  <=  N  <=  128 

1  <=  N  <=  128 

0  <=  nstop  <=  999 


THE  NUMBER  OF  SYSTEMS  TO  BE  EVALUATED. 
THIS  INTEGER  VALUE  MUST  OCCUR  AT  THE  TOP  OF  THE  INPUT 
FILE.  IT  DELINEATES  THE  NUMBER  OF  SYSTEMS  TO  BE  READ  BY 
THE  PROGRAM  AND  ANALYZED.   FOR  EACH  SYSTEM  (1... numsys) 
THE  PARAMETERS  BELOW  MUST  APPEAR  IN  THE  INPUT  FILE. 


L  =  AN  INTEGER  VALUE  THAT  SPECIFIES  THE  MAXIMUM  NUMBER  OF  DELAYS 
IN  THE  INPUT  SEQUENCE. 

N  =  AN  INTEGER  VALUE  THAT  SPECIFIES  THE  MAXIMUM  NUMBER  OF  DELAYS 
IN  THE  OUTPUT  SEQUENCE. 

nstop  =  AN  INTEGER  VALUE  THAT  SPECIFIES  THE  LARGEST  TIME  INDEX 
(ns)  FOR  WHICH  THE  DIFFERENCE  EQUATION  IS  TO  BE  SOLVED. 

xsorce  =  A  CHARACTER  VALUE  OF  'F'  OR  'S'  DENOTING  WHETHER  THE 

INPUT  SEQUENCE  x( )  IS  TO  BE  READ  FROM  THE  INPUT  FILE  (F) 
OR  TO  BE  GENERATED  (S)  USING  THE  SUBROUTINE  xgen.   THIS 
LATER  OPTION  IS  ATTRACTIVE  WHEN  nstop  IS  A  LARGE  NUMBER 
AND  THE  INPUT  SEQUENCE  x( )  CAN  BE  READILY  DESCRIBED  BY  AN 
ANALYTICAL  EXPRESSION.   IF  xsorce  =  'S*  THE  USER  MUST 
PROVIDE  THE  APPROPRIATE  FORTRAN  STATEMENTS  IN  THE  SPACE 
PROVIDED  IN  SUBROUTINE  xgen  AND  THE  PROGRAM  MUST  BE 
RECOMPILED  BEFORE  EXECUTION. 

b(k)  =  REAL  COEFFICIENTS  OF  THE  INPUT  SEQUENCE  x(ns-k)  IN  THE 
ORDER:   b(0),  b(l),  .  .  .  ,  b(L). 

a(k)  =  REAL  COEFFICIENTS  OF  THE  OUTPUT  SEQUENCE  y(ns-k)  IN  THE 

ORDER:   a(l),  a(2),  ... ,  a(N).   IF  N  =  0  THEN  THE  EQUATION 
IS  NON-RECURSIVE  AND  NO  a(k)  COEFFICIENTS  SHOULD  BE  IN 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c* 

c 

c 

c 

c 

c 


THE  INPUT  FILE. 

y(k)  =  THE  INITIAL  CONDITIONS  FOR  THE  OUTPUT  SEQUENCE  IF  THE 
DIFFERENCE  EQUATION  IS  RECURSIVE,  I.E.  ,  N  >  0.   THIS 
PROGRAM  CALCULATES  THE  SOLUTION  TO  THE  DIFFERENCE 
EQUATION  FROM  ns  =  0  TO  ns  =  nstop  THEREFORE  THE  INITIAL 
CONDITIONS  y(-N)  TO  y(-l)  MUST  BE  PROVIDED  IN  THE  INPUT 
FILE  IN  THE  ORDER:  y(-N),  y(-N+l),  ...,  y(-l).   IF  N  =  0 
THEN  THE  EQUATION  IS  NON-RECURSIVE  AND  NO  INITIAL 
CONDITIONS  SHOULD  BE  GIVEN  IN  THE  INPUT  FILE. 

x(ns)  =  THE  INPUT  SEQUENCE.   IF  xsorce  =  'FT  THEN  THE  INPUT 

SEQUENCE  x(0),  ... ,  x( nstop)  MUST  BE  PROVIDED  BY  THE  USER 
IN  THE  INPUT  FILE.   IF  xsorce  =  'S*  THEN  THE  USER  HAS 
ELECTED  TO  GENERATE  THE  INPUT  SEQUENCE  BY  PROVIDING  THE 
APPROPRIATE  FORTRAN  STATEMENTS  IN  THE  SUBROUTINE  xgen. 

NOTE:   THE  INPUT  FORMAT  STATEMENTS  OCCUR  IN  THE  MAIN  PROGRAM 

FOLLOWING  THE  CAPTION:  ********   INPUT  FORMAT  ********. 
THE  FORM  OF  THE  INPUT  DATA  FILE  IS: 


LINE  # 

1 

2 
NEXT  NB  LINES 
NEXT  NA  LINES 
NEXT  NY  LINES 
NEXT  NX  LINES 


ENTRIES 

numsys 

L,N, nstop, xsorce 

b(k),  k=0,l,... ,L 

a(k),  k=l,... ,N 

y(k),  k=  -N,... ,-1 

x(ns) ,  ns=  0, . . . , nstop 


FORMAT 

il 

i3,tll,i3,t21,i3,t31,al 
6fl0.0 
6fl0.  0 
6fl0.  0 
6fl0.  0 


WHERE:  NB  =  1  +  (L/6  ROUNDED  DOWN  TO  THE  NEXT  SMALLER  INTEGER) 

NA  =  0  IF  N  =  0  OR 

NA  =  1  +  ((N-l)/6  ROUNDED  DOWN  TO  THE  NEXT  SMALLER  INTEGER) 

NY  =  0  IF  N  =  0  OR 

NY  =  1  +  ((N-l)/6  ROUNDED  DOWN  TO  THE  NEXT  SMALLER  INTEGER) 

NX  =  0  IF  xsorce  =  'S'  OR 

NX  =  1  +  (nstop/ 6  ROUNDED  DOWN  TO  THE  NEXT  SMALLER  INTEGER) 
IF  xsorce  =  'F* 


"NOTE:   FOR  numsys  >  1  THE  FORMAT  OF  LINES  2. 


IS  REPEATED. 


THE  FORM.  T  f 10. 0  USED  FOR  INPUT  DATA  PERMITS  THE  DECIMAL 
POINT  TO  BE  PLACED  ANYWHERE  IN  THE  FIELD  OF  10  COLUMNS 
AND  ALSO  ALLOWS  THE  EXPONENTIAL  FORMAT  TO  BE  USED  (EG. 
3146.2  =  ?<.  1462E+03). 


r****************************  OUTPUT   ****************************** 


THE  INPUT  DATA  AS  WELL  AS  THE  OUTPUT  DATA  ARE  STORED  IN  TABULAR 
FORM  IN  THE  FILE  ' DIFFEQ. OUT' .   ADDITIONALLY,  THE  INPUT  SEQUENCE 
AND  THE  OUTPUT  SEQUENCE  ARE  WRITTEN  INTO  THE  FILE  'DIFFEQ.DAT'  TO 
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FACILITATE  PLOTTING  BY  A  SEPARATE,  USER  SUPPLIED  PROGRAM.   THE 
FORMAT  OF  THE  DATA  IN  'DIFFEQ.DAT'  IS:  el2. 6,  2X,  el2. 6.  THE  FIRST 
ENTRY  CORRESPONDS  TO  THE  ORDINATE  VALUE,  AND  THE  SECOND  ENTRY,  THE 
ABSCISSA  VALUE.   ADDITIONAL  HEADER  INFORMATION  IS  WRITTEN  INTO 
'DIFFEQ.DAT'  TO  ALLOW  FOR  CONTROL  AND  LABELING  OF  EACH  PLOT. 


■fc'ferkirkicic'ffie'icitic'kirfciciticir/ciriticirkirie'ieicie      EXAMPLE      ^MeiricifMrMrMr^itifirkiticitirkirkirkie'irfrk 

THE  INPUT  PARAMETERS  FOR  THE  SYSTEM  DESCRIBED  BELOW  ARE  STORED  IN 
THE  SAMPLE  INPUT  FILE  'DIFFEQ. TST'  AND  CAN  BE  USED  FOR  A  TRIAL 
RUN  IN  TEST  MODE. 

DIFFERENCE  EQUATION: 

y(ns)  =  1.2*y(ns-l)  +  1. 5*x(ns) 

GOAL:     TO  OBTAIN  THE  SOLUTION  TO  THIS  DIFFERENCE  EQUATION  FOR 
ns  =  0  TO  ns  =  10,  GIVEN:  x( 0) . . . x( 10)  =  100.0  AND 
THE  INITIAL  CONDITION  y(-l)  =  25.0. 

THE  INPUT  FILE  IS: 


1 

000 
1.  5 
1.  2 
25.  0 
100.  0 
100.  0 


001 


100.  0 
100.  0 


010 


100.0 
100.  0 


100.  0 
100.  0 


100.0 
100.0 


100.  0 


THE  RESULTING  OUTPUT  FILE  'DIFFEQ. OUT'  IS: 

INPUT  DATA  FOR  PROBLEM  #   1 

PROBLEM  #  1    INPUT  DATA  SOURCEFILE:  DIFFEQ. TST 

THE  NUMBER  OF  INPUT  DELAYS:  L  =   0 

THE  NUMBER  OF  OUTPUT  DELAYS:  N  =   1 

THE  VALUE  OF  ns top  IS:   10 

THE  COEFFICIENTS  b(0),  b(l),  ...,  b(L)  ARE: 

. 150000E+01 

THE  COEFFICIENTS  a(l),  ...,  a(N)  ARE: 
. 120000E+01 

OUTPUT  DATA  FOR  PROBLEM  #   1 


ns 
-1 


x(ns) 
OOOOCOE+00 


y(ns) 
250000E+02 
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C  0      .  100000E+03        .  180000E+03 

C  1      .  100000E+03        . 366000E+03 

C  2      . 100000E+03        .589200E+03 

C  3      .  100000E+03        .857040E+03 

C  4      .  100000E+03        .  117845E+04 

C  5      . 100000E+03        . 156414E+04 

C  6      . 100000E+03        .202697E+04 

C  7      . 100000E+03        . 258236E+04 

C  8      . 100000E+03        .324883E+04 

C  9      .  100000E+03        .404860E+04 

C  10      .  100000E+03        .500832E+04 
C 

C  END  OF  PROBLEM  #  1   

C 

C 

C  FOR  ILLUSTRATIVE  PURPOSES  THE  INPUT  SEQUENCE  x( )  COULD  HAVE  BEEN 

C  GENERATED  BY  SPECIFYING  xsorce  =  'S*  AND  WRITING  THE  APPROPRIATE 

C  FORTRAN  STATEMENTS  INTO  SUBROUTINE  xgen.   THE  STATEMENTS  THAT 

C  COULD  BE  USED  TO  ACCOMPLISH  THIS  ARE  WRITTEN  INTO  THE  SUBROUTINE 

C  BUT  ARE  'COMMENTED  OUT'. 

C 

C 

Q***************************   MAIN  PROGRAM   **************************** 


character  inf ile*12,  mode,vl,  xsorce*l 

real  a(l:128),  b(0:128),  y(-128:850),  x(-128:850),  ii 

C   PROMPT  USER  FOR  MODE:  BATCH  OR  TEST. 

write(»,1115) 

read(* , 1117)  mode 

if  ( (mode.  eq.  '  Y'  ).  or.  (mode.  eq.  'y '  ))  then 
mode  =  'Y' 
write(-'',1118) 
read(-'-,1119)  infile 

else 
infile  =  'DIFFEQ. IN' 

endif 

C  UNIT=1  DEFINED  AS  INPUT  FILE.   UNITS=2,3  DEFINED  AS  OUTPUT  FILES. 

open(unit=l j  f ile=inf ile,status=' old' , iostat=ierr ,err=999) 
open(unit=2 , f ile=' DIFFEQ.  OUT* ) 
open(unit=3,file=' DIFFEQ.  DAT') 

C  READ  INPUT  PARAMETERS  AND  CONDUCT  ERROR  CHECKS. 

read( 1,1000)  numsys 

if((numsys.  le.  0).  or.  (numsys.  gt.  4))  then 
write( " , 1126)  numsys 
stop  'The  allowed  values  for  numsys  are  1  <=  numsys  <=  4. 

endif 

numplts   =  2-numsys 

wrire(3,2000)    numplts 
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do  10  nprob=l,  numsys 

read(l,1001)  L,  N,  nstop,  xsorce 

if((L.  lt.0).or.  (L.  gt.  128))  then 
write(*,1124)  nprob,  'L',  L 

stop  'The  allowed  values  for  "L1'  are:   0  <=  L  <=  128.' 
endif 

if((N.  It.  0).or.  (N.  gt.  128))  then 
write(*,1124)  nprob,  'N' ,  N 

stop  'The  allowed  values  for  "N"  are:   0  <=  N  <=  128.' 
endif 

if ((nstop. It. 0).  or.  (nstop.  gt.  850))  then 
write( " , 1127)  nprob,  nstop 

stop  'The  allowed  values  for  nstop  are:  0  <=  nstop  <=  850. 
endif 

if  (  (xsorce.  eq.  '  F'  ).  or.  (xsorce.  eq.  '  f '  )  )  then 

xsorce  =  'F' 
elseif  ((xsorce.  eq.  '  S'  ).  or.  (xsorce.  eq.  '  s'  ))  then 

xsorce  =  'S' 
else 

write(-, 1128)  nprob,  xsorce 

l™        ,-,       j       -,         j-       It  It  t  I  -p,  I  t         I  Ipl  I    I 

stop  The  allowed  values  for   xsorce   are:    F   or   S 
endif 

C   INITIALIZE  EACH  ARRAY  TO  ZERO  BEFORE  EACH  RUN. 

data  a/128*0. 0/,  b/129*0. 0/ 
data  y/979*0. 0/,  x/979*0. 0/ 

C  READ  THE  COEFFICIENTS  b( ) ,  a( )  AND  THE  INITIAL  CONDITIONS 
C  y(-N)...y(-l). 

read(l,1002)  (b(k),  k=0,L) 
if(N.  gt.  0)  then 

icstart  =  -N 

read(l,1002)  (a(k),  k=l,N) 

read( 1,1002)  (y(k),  k=icstart , -1) 
endif 

C  FOR  xsorce  =  'F*  READ  THE  ARRAY  x( )  FROM  THE  INPUT  FILE. 
C  FOR  xsorce  =  'S'  CALL  xgen  TO  GENERATE  THE  ARRAY  x(). 

if(xsorce.  eq.  'F' )  then 

read( 1,1002)  (x(k),  k=0, nstop) 

else 
call  xgen(x, nstop, nprob) 
endif 

C  FOR  TEST  MODE  ECHO  INPUT  PARAMETERS  ONTO  MONITOR  (UNIT  =  *)• 

if  (mode.  eq.  '  Y' )  then 
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write(*,1007) 

write(*, 1120)  nprob,  infile 
write(*,1110)  'INPUT*   'L'   L 
write(*> 1110)  'OUTUT  ,  V,  N 
write(*,1112)  nstcp 
write(*,1004) 

write(*,1005)  (b(k),k=0,L) 
if(N.  eq.  0)  then 

write(*,113i) 
else 

write(*,1006) 

write(*,1005)  (a(k),k=l,N) 

end  if 
write(*,1123)  nprob 
pause  'END  OF  RUN,  STRIKE  <CR>  WHEN  READY  TO  CONTINUE.  ' 

endif 

C  WRITE    INPUT  DATA  INTO  FILE:  DIFFEQ.  OUT. 

write(2,1008)  'INPUT',  nprob 
write(2,1120)  nprob,  infile 
write(2,1110)  'INPUT'   'L'   L 
write(2,1110)  'OUTPUT*,  V,  N 
write(2,1112)  nstop 
write(2,1004) 

write(2,1005)  (b(k),k=0,L) 
if(N.  eq.  0)  then 

write(2,1131) 
else 

vrite(2,1006) 

write(2,1005)  (a(k),k=l,N) 
endif 

C  WRITE  THE  INPUT  SEQUENCE  INTO  FILE:  DIFFEQ.DAT. 

write(3,2001)  nstop  +  1 
write(3, •-•-)  'INPUT  SEQUENCE  x(ns)' 
write(3,*)  'SAMPLE  #  (ns)' 
write(3,")  'x(ns) ' 
do  54  ns=0,  nstop 

ii  =  ns 

write(3,2010)  ii,  x(ns) 
54      continue 

C  CALL  diffeq  TO  COMPUTE  THE  SOLUTION  TO  THE  DIFFERENCE  EQUATION. 

call  dif feq(N,L,a,b,x,y,nstop) 

C  WRITE  RESULTS  INTO  FILE:  DIFFEQ.DAT. 

write(3,2001)  N  +  nstop  +  1 
write(3,*)  'OUTPUT  SEQUENCE  Y(ns)' 
write(3,*)  'SAMPLE  #  (ns)' 
write(3,")  'y(ns)' 
do  55  ns=  -N,  nstop 
ii  =  ns 
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write(3,2010)  ii,  y(ns) 
55      continue 

C  WRITE  RESULTS  INTO  FILE:  DIFFEQ.  OUT. 

write(2,1008)  'OUTPUT' ,  nprob 

write(2,1129) 

do  102  ns=  -N,  nstop 

write(2, 1130)  ns ,  x(ns),  y(ns) 
102     continue 
write(2,1123)  nprob 

10    continue 

write(*,1121)~. 

999  close(unit=l) 
close(unit=2) 
close(unit=3) 

if  (  ierr.  gt.  0)  then 
write(* ,1116)  ierr 
endif 

C********      INPUT  FORMAT  ******** 

1000  format( il) 

1001  format(i3,tll,i3,t21,i3,t31,al) 

1002  forraat(6fl0.  0) 


1004  format(t4,'THE   COEFFICIENTS  b(0),    b(l),    ...,    b(L)   ARE:',/) 

1005  format(6(lx,el2.  6)) 

1006  format(//,t4,'THE  COEFFICIENTS  a(l),  ...,  a(N)  ARE:',/) 

1007  format(//////////) 

1008  formatr///, tl6,a6,'  DATA  FOR  PROBLEM  #  ',il,//) 

1110   format( t4,' THE  NUMBER  OF  ',a6,'  DELAYS:  ',al,'  =  ',i3) 
1112   formatet4,'THE  VALUE  OF  nstop  IS:  ',i3) 

1115  0format(lx,'D0  YOU  WISH  TO  RUN  THIS  PROGRAM  IN  TEST' , 
1'  MODE  ?     (Y/N)  <CR>  :  ' ,  ,) 

1116  0format(///,lx, 'ERROR  OPENING  INPUT  FILE,  PROGRAM  TERMINATED.', 
l//,lx, 'ERROR  CODE: ' ,i4, /////) 

1117  formateal) 

1118  0format(/////,lx,'TYPE  THE  NAME  OF  YOUR  DATA  FILE  FOLLOWED', 
1'  BY  <CR>.  '  ,/,'  IF  YOU  DESIRE  TO  MAKE  A  TEST  RUN  USING  THE' , 
2'  SAMPLE  DATA  ALREADY  STORED',/,'  IN  THE  FILE:  DIFFEQ. TST', 
3  '       TYPE:  DIFFEQ.  TST  <CR>' ,/,'  FILENAME:  ',  ,) 

1119  formateal2) 

1120  f ormate///,t4,' PROBLEM  #' ,il,'    INPUT  DATA  S0URCEF I LE:  * ,al2) 

1121  Oformate//,'  TABULAR  OUTPUT  DATA  IS  STORED  IN  FILE:  DIFFEQ.  OUT.  ' 
1,/,'  PLOTTING  DATA  IS  STORED  IN  FILE:  DIFFEQ.DAT.') 

1122  format( i3) 

1123  format(/, lx,16( '-'),'  END  OF  PROBLEM  #* ,i2,2x,16( ' -' ) ,//) 

1124  Oformat(//,'  For  problem  #',i2,'  the  value  for  ',al,'  is:  ',i3, 
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l'.   This  value  is  not  allowed.') 

1126  format(//,'  numsys  =  ',i4,'.   This  value  is  not  allowed.') 

1127  Oformat(//,  '  For  problem  #',i2,'  the  value  for  ''nstop''  is: 
li3,'.   This  value  is  not  allowed.') 

1128  Oformat(//,*  For  problem  #',i2,'  the  value  for  "xsorce"  is: 
lal,'.   This  value  is  not  allowed.') 

1129  format (t6,'ns' ,tl6,'x(ns)' , t35 , ' y(ns) ' ) 

1130  format(t4,i4,tll,el4. 6,t30,el4. 6) 

1131  format(/,'  THIS  SYSTEM  IS  NON -RECURSIVE,  I.E.  ,  N  =  0.  ' ) 

2000  format( il) 

2001  forraat( i3) 

2010   formateel2. 6,2x,el2. 6) 


end 


SUBROUTINE:  diffeq 


C  PURPOSE:   THIS  SUBROUTINE  COMPUTES  THE  SOLUTION  TO  A  DIFFERENCE 
C  EQUATION.   ALL  PARAMETERS  DESCRIBING  THE  EQUATION,  AND 

C  THE  INPUT  AND  OUTPUT  SEQUENCES  xe )  AND  ye)  ARE  PASSED 

C  TO  THE  SUBROUTINE  BY  THE  MAIN  PROGRAM. 

subroutine  dif feq(N,L,a,b,x,y,nstop) 

real  x(  -128:  nr  top)  ,  y(  -128:  nstop) ,  a(l:N),  b(0:L) 

do  500  ns=0,  nstop 
y(ns)  =  0.  0 
do  501  k=0,  max(N,Lj 

y(ns)  =  y(ns)  +  a(k)*y(ns-k)  +  b(k)*x(ns-k) 
501     continue 
500   continue 

return 
end 


SUBROUTINE:  xgen 


C  PURPOSE:   THIS  SUBROUTINE  ALLOWS  THE  USER  TO  GENERATE  VALUES  FOR 
C  THE  ARRAY  x().   IF  xsorce  =  'S'  THE  MAIN  PROGRAM  WILL 

C  CALL  THIS  SUBROUTINE.   IF  xsorce  =  'F'  THIS  SUBROUTINE 

C  WILL  NOT  BE  CALLED  BY  THE  MAIN  PROGRAM. 


subroutine  xgenex, nstop, nprob) 
real  x( -128: nstop) 


C  DEVELOP  THE  ALGORITHM  FOR  GENERATING  VALUES  OF  x( )  IN  THIS  SPACE. 

C  THE  STATEMENTS  TYPED  IN  MUST  FOLLOW  STANDARD  FORTRAN  77  RULES  AND 

C  MAY  USE  FORTRAN  77  INTRINSIC  FUNCTIONS  SUCH  AS:  SIN(),  C0S( ) ,  ... 

C  NOTE  THAT  THE  VALUE  nprob  CAN  BE  USED  IN  A  LOGICAL  'IF'  STATEMENT 

C  TO  MATCH  THE  GENERATING  FUNCTIONS  TO  THE  CORRESPONDING  SYSTEM 
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C  EQUATION  READ  FROM  THE  INPUT  FILE  IF  MORE  THAN  ONE  SYSTEM  OF 

C  EQUATIONS  EXIST.   AN  EXAMPLE  OF  AN  ALGORITHM  GENERATING  VALUES 

C  FOR  x()  IS: 
C 

c 

C***  EXAMPLE  *** 

C 

C     if(nprob.  eq.  1)  then 

C       do  1  k=0,  nstop 

C         x(k)  =  100.  0 

C  1     continue 

C     endif 

x(0)  =  1.  0 

do  1  k=l, nstop 
x(k)  =  0.0 
1     continue 


return 
end 


B.     INPUT,  OUTPUT  DATA  OF  SECOND  ORDER  IIR  NOTCH  FILTER:  P  =  0 


INPUT  DATA  FOR  PROBLEM  #  1 


PROBLEM  #  1   INPUT  DATA  SOURCEFILE:  dif f eq.  tst 

THE  NUMBER  OF   INPUT  DELAYS:  L  =   2 

THE  NUMBER  OF  OUTPUT  DELAYS:  N  =   2 

THE  VALUE  OF  nstop  IS:  400 

THE  COEFFICIENTS  b(0),  b(l),  ...,  b(L)  ARE: 

990000E+00  -.329930E+00   . 990000E+00 
THE  COEFFICIENTS  a(l),  ...,  a(N)  ARE: 
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329930E+00  -. 980100E+00 


OUTPUT  DATA  FOR  PROBLEM  #  1 


ns 

x(ns) 

y(ns) 

-2 

. 000000E+00 

. OOOOOOE+OO 

-1 

.  000000E+00 

. OOOOOOE+OO 

0 

. 100000E+01 

. 990000E+00 

1 

. 000000E+00 

-. 329930E-02 

2 

. 000000E+00 

.  186125E-01 

3 

. 000000E+00 

. 937446E-02 

4 

. 000000E+00 

-.  151492E-01 

5 

. 000000E+00 

-. 141861E-01 

6 

. 000000E+00 

. 101673E-01 

7 

. OOOOOOE+OO 

. 172583E-01 

8 

. 000000E+00 

-.427095E-02 

9 

. 000000E+00 

-. 183239E-01 

10 

. 000000E+00 

-. 185966E-02 

11 

. 000000E+00 

. 173457E-01 

12 

. 000C00E+00 

. 754553E-02 

13 

. 000000E+00 

-.  145111E-01 

14 

. OCOOOOE+00 

-.  121830E-01 

15 

. OOOOOOE+OO 

. 102028E-01 

16 

. OOOOOOE+OO 

. 153068E-01 

17 

. OOOOOOE+OO 

-.494956E-02 

18 

. OOOOOOE+OO 

-. 166352E-01 

19 

. OOOOOOE+OO 

-. 637380E-03 

20 

. OOOOOOE+OO 

.  160938E-01 

21 

.OOOOOOE+OO 

. 593453E-02 

22 

. OOOCOOE+00 

-. 138156E-01 

23 

. OOCOOOE+00 

-. 103746E-01 

24 

. OOOOOOE+OO 

.  101178E-01 

25 

. OOOOOOE+OO 

. 135063E-01 

26 

. OOOOOOE+OO 

-. 546028E-02 

27 

. OOOOOOE+OO 

-. 150390E-01 

28 

. OOOOOOE+OO 

. 389786E-03 

29 

. OOOOOOE+OO 

. 148684E-01 

30 

. OOOOOOE+OO 

.452349E-02 

31 

. OOOOOOE+OO 

-. 130801E-01 

32 

.OOOOOOE+OO 

-. 874898E-02 

33 

. OOOOOOE+OO 

.993321E-02 

34 

. OOOOOOE+OO 

. 118521E-01 

35 

. OOOOOOE+OO 

-.582516E-02 

36 

. OOOOOCE+00 

-.  135382E-01 

37 

. OOOOO.I+OO 

. 124259E-02 

38 

. OOOOOOE+OO 

. 136787E-01 

39 

. OOOOOOE+OO 

. 329516E-02 

40 

. oooonoE+oo 

-. 123194E-01 

41 

. OOOOOOE+OO 

-.  729411E-02 

42 

. OOOOOOE+OO 

. 966765E-02 

43 

. OOOOOOE+OO 

.  103386E-01 

44 

. OOOOOOE+OO 

-. 606425E-02 

45 

. OOOOOOE+OO 

-. 121336E-01 
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46 

. 000000E+00 

.  194032E-02 

47 

. 000000E+00 

. 125324E-01 

48 

. 000000E+00 

. 223309E-02 

49 

. 000000E+00 

-. 115462E-01 

50 

. 000000E+00 

-. 599809E-02 

51 

. 000000E+00 

. 933748E-02 

52 

. 000000E+00 

. 895944E-02 

53 

. 000000E+00 

-.619567E-02 

54 

. 000000E+00 

-. 108253E-01 

55 

. 000000E+00 

. 250079E-02 

56 

. 000000E+00 

. 114349E-01 

57 

. 000000E+00 

. 132171E-02 

58 

. 000000E+00 

-. 107713E-01 

59 

. 000000E+00 

-.484919E-02 

60 

. 000000E+00 

.895708E-02 

61 

. 000000E+00 

. 770790E-02 

62 

. 000000E+00 

-.623577E-02 

63 

. 000000E+00 

-. 961188E-02 

64 

. 000000E+00 

. 294043E-02 

65 

. 000000E+00 

. 103907E-01 

66 

. 000000E+00 

.546304E-03 

67 

. 000000E+00 

-. 100037E-01 

68 

. 000000E+00 

-. 383596E-02 

69 

. 000000E+00 

. 853905E-02 

70 

. 000000E+00 

. 657691E-02 

71 

. 00O0OOE+O0 

-.619920E-02 

72 

. 000000E+00 

-. 849133E-02 

73 

. 000000E+00 

. 327429E-02 

74 

. 00000OE+OO 

. 940264E-02 

75 

. 000000E+00 

-.  106918E-03 

76 

. OOOOOOE+OO 

-. 925080E-02 

77 

. OOOOOOE+OO 

-. 294733E-02 

78 

. OOOOOOE+OO 

. 809430E-02 

79 

. OOOOOOE+OO 

. 555923E-02 

80 

. OOOOOOE+OO 

-. 609907E-02 

81 

. OOOOOOE+OO 

-. 746086E-02 

82 

. OOOOOOE+OO 

. 351613E-02 

83 

. OOOOOOE+OO 

. 847247E-02 

84 

. OOOOOOE+OO 

-. 650839E-03 

85 

. OOOOOOE+OO 

-.  851860E-02 

86 

. OOOOOOE+OO 

-.  217265E-02 

87 

. OOOOOOE+OO 

. 763226E-02 

88 

. OOOOOOE+OO 

.464753E-02 

89 

. OOOOOOE+OO 

-.594701E-02 

90 

. OOOOOOE+OO 

-.651714E-02 

91 

. OOOOOOE+OO 

. 367847E-02 

92 

. OOOOOOE+OO 

. 760109E-02 

93 

. OOOOOOE+OO 

-. 109744E-02 

94 

. OOOOOOE+OO 

-. 781190E-02 

95 

. OOOOOOE+OO 

-. 150178E-02 

96 

. OOOOOOE+OO 

. 716096E-02 

97 

.OOOOOOE+OO 

. 383451E-02 

98 

. OOOOOOE+OO 

-.575334E-02 

99 

.OOOOOOE+OO 

-.565640E-02 

100 

. OOOOOOE+OO 

.377263E-02 

101 

. OOOOOOE+OO 

. 678855E-02 
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102 

. 000000E+00 

-. 145781E-02 

103 

. 000000E+00 

-. 713443E-02 

104 

. 000000E+O0 

-. 925062E-03 

105 

. OOOOOOE+00 

.668725E-02 

106 

. 000000E+00 

. 311298E-02 

107 

. OOOOOOE+OO 

-.552711E-02 

108 

.OOOOOOE+OO 

-.487459E-02 

109 

.  OOOOOOE+OO 

. 380885E-02 

110 

. OOOOOOE+OO 

. 603424E-02 

111 

.OOOOOOE+OO 

-. 174217E-02 

112 

. OOOOOOE+OO 

-. 648895E-02 

113 

.OOOOOOE+OO 

-. 433395E-03 

114 

. OOOOOOE+OO 

. 621683E-02 

115 

. OOOOOOE+OO 

. 247589E-02 

116 

.  OOOOOOE+OO 

-.527624E-02 

117 

.  OOOOOOE+OO 

-.416741E-02 

118 

. OOOOOOE+OO 

. 379629E-02 

119 

. OOOOOOE+OO 

.533699E-02 

120 

.OOOOOOE+OO 

-. 195991E-02 

121 

.  OOOOOOE+OO 

-.587742E-02 

122 

. OOOOOOE+OO 

-. 182246E-04 

123 

. OOOOOOE+OO 

. 575444E-02 

124 

.  OOOOOOE+OO 

. 191643E-02 

125 

. OOOOOOE+OO 

-. 500764E-02 

126 

. OOOOOOE+OO 

-. 353046E-02 

127 

. OOOOOOE+OO 

. 374319E-02 

128 

.  OOOOOOE+OO 

.469519E-02 

129 

. OOOOOOE+OO 

-.  211961E-02 

130 

. OOOOOOE+OO 

-. 530108E-02 

131 

.  OOOOOOE+OO 

. 328445E-03 

132 

. OOOOOOE+OO 

.530396E-02 

133 

. OOOOOOE+OO 

. 142802E-02 

134 

. OOOOOOE+OO 

-.472726E-02 

135 

. OOOOOOE+OO 

-. 295927E-02 

136 

. OOOOOOE+OO 

. 365683E-02 

157 

. OOOOOOE+OO 

.410688E-02 

138 

. OOOOCOE+00 

-. 222908E-02 

139 

. OOOOOOE+OO 

-.476059E-02 

140 

. OOOOCOE+00 

. 614057E-03 

141 

. OOOOOOE+OO 

.486845E-02 

142 

. OOOOOOE+OO 

. 100441E-02 

143 

. OOOOOOE+OO 

-.444019E-02 

144 

. OOOOOOE+OO 

-. 244937E-02 

145 

. OOOOOOE+OO 

. 354370E-02 

146 

.  OOOOOOE+OO 

. 356981E-02 

147 

. OOOOOOE+OO 

-. 229540E-02 

148 

. OOOOOOE+OO 

-.425609E-02 

149 

. OOOOOOE+OO 

. 845509E-03 

150 

. OOOOOOE+OO 

.445035E-02 

151 

. OOOOOOE+OO 

.639621E-03 

152 

. OOOOOOE+OO 

-.415076E-02 

153 

. OOOOOOE+OO 

-. 199635E-02 

154 

. OOOOOOE+OO 

. 340950E-02 

155 

. OOOOOOE+OO 

. 308152E-02 

156 

. OOOOOOE+OO 

-. 232497E-02 

157 

. OOOOOOE+OO 

-. 378727E-02 
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158 

. O00000E+00 

.  102916E- 

-02 

159 

. 000000E+00 

.405146E- 

■02 

160 

. 000000E+00 

. 328015E- 

■03 

161 

.O00O0OE+0O 

-. 386261E- 

■02 

162 

. 000000E+00 

-.  159588E- 

■02 

163 

. OOOOOOE+00 

. 325922E- 

■02 

164 

. OOOOOOE+00 

. 263944E- 

-02 

165 

. OOOOOOE+00 

-. 232353E- 

-02 

166 

. OOOOOOE+00 

-.  335351E- 

•02 

167 

. OOOOOOE+00 

.  117087E- 

-02 

168 

. OOOOOOE+00 

.  367308E- 

-02 

169 

. OOOOOOE+OO 

. 642924E- 

■04 

170 

. OOOOOOE+00 

-.357878E- 

■02 

171 

. OOOOOOE+OO 

-. 124376E- 

■02 

172 

. OOOOOOE+00 

.309721E- 

■02 

173 

. OOOOOOE+00 

.224087E- 

-02 

174 

. OOOOOOE+00 

-.  229624E- 

■02 

175 

. OOOOOOE+OO 

-. 295387E- 

■02 

176 

. OOOOOOE+00 

.  127597E- 

■02 

177 

. OOOOOOE+OO 

. 331607E- 

-02 

178 

. OOOOOOE+00 

-.  156510E- 

■03 

179 

. OOOOOOE+OO 

-. 330172E- 

■02 

180 

.OOOOOOE+00 

-.  935942E- 

-03 

181 

. OOOOOOE+00 

. 292722E- 

-02 

182 

. OOOOOOE+OO 

.  188310E- 

-02 

183 

. OOOOOOE+00 

-.  224768E- 

■02 

184 

. OOOOOOE+00 

-. 258720E- 

-02 

185 

. OOOOOOE+OO 

.  134936E- 

-02 

186 

. OOOOOOE+00 

. 298091E- 

■02 

187 

. OOOOOOE+OO 

-.  339014E- 

■03 

188 

. OOOOOOE+00 

-.  303344E- 

-02 

189 

. OOOOOOE+00 

-. 668554E- 

-03 

190 

. OOOOOOE+OO 

. 275250E- 

-02 

191 

. OOOOOOE+00 

. 156338E- 

■02 

192 

. OOOOOOE+00 

-.  218192E- 

-02 

193 

. OOOOOOE+OO 

-.225215E- 

-02 

194 

. OOOOOOE+00 

.  139544E- 

-02 

195 

. OOOOOOE+00 

. 266773E- 

-02 

196 

. OOOOOOE+OO 

-.487510E- 

-03 

197 

. OOOOOOE+00 

-. 277549E- 

-02 

198 

. OOOOOOE+OO 

-.437908E- 

-03 

199 

. OOOOOOE+00 

. 257578E- 

-02 

200 

. OOOOOOE+OO 

.  127902E- 

■02 

201 

. OOOOOOE+00 

-.210253E- 

-02 

202 

. OOOOOOE+00 

-. 194725E- 

■02 

203 

. OOOOOOE+OO 

. 141823E- 

-02 

204 

.OOOOOOE+00 

.237642E- 

-02 

205 

. OOOOOOE+00 

-.  605957E- 

■03 

206 

. OOOOOOE+00 

-. 252905E- 

■02 

207 

.OOOOOOE+00 

-.240512E- 

■03 

208 

.OOOOOOE+00 

.239937E- 

■02 

209 

. OOOOOOE+OO 

. 102735E- 

■02 

210 

. OOOOOOE+00 

-.201267E- 

■02 

211 

. OOOOOOE+00 

-. 167095E- 

■02 

212 

. OOOOOOE+00 

. 142132E- 

■02 

213 

. OOOOOOE+OO 

. 210663E- 

-02 
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214 

.  000000E+00 

-. 697998E- 

■03 

215 

. 000000E+00 

-. 229500E- 

•02 

216 

. O00000E+0O 

-. 730823E- 

■04 

217 

.  000000E+00 

. 222522E- 

■02 

218 

. 000000E+00 

. 805795E- 

■03 

219 

. 000000E+00 

-.  191508E- 

•02 

220 

. OOOOOOE+OO 

-. 142160E- 

■02 

221 

. ooooocz+oo 

. 140794E- 

•02 

222 

. 000000E+00 

.  185783E- 

■02 

223 

. 000000E+00 

-. 766969E- 

■03 

224 

. 000OQ0E+00 

-.207391E- 

-02 

225 

. 000000E+00 

.674610E- 

■04 

226 

. 000000E+00 

. 205490E- 

■02 

227 

. OOOOOOE+OO 

. 611853E- 

■03 

228 

. 000000E+00 

-. 181214E- 

•02 

229 

. OOOOOOE+OO 

-.  119756E- 

■02 

230 

. OOOOOOE+OO 

.  138096E- 

■02 

231 

. OOOOOOE+OO 

.  162935E- 

■02 

232 

.  OOOOOOE+OO 

-.  815913E- 

■03 

233 

. OOOOOOE+OO 

-. 186612E- 

■02 

234 

. OOOOOOE+OO 

. 183989E- 

■03 

235 

. OOOOOOE+OO 

. 188968E- 

■02 

236 

.  OOOOOOE+OO 

.443136E- 

■03 

237 

. OOOOOOE+OO 

-. 170587E- 

■02 

238 

. OOOOOOE+OO 

-. 997137E- 

■03 

239 

.  OOOOOOE+OO 

. 134294E- 

■02 

240 

. OOOOOOE+OO 

. 142037E- 

•02 

241 

. OOOOOOE+OO 

-.847595E- 

■03 

242 

. OOOOOOE+OO 

-. 167175E- 

■02 

243 

. OOOOOOE+OO 

.279167E- 

■03 

244 

. OOOOOOE+OO 

. 173059E- 

■02 

245 

. OOOOOOE+OO 

. 297362E- 

■03 

246 

. OOOOOOE+OO 

-.  159804E- 

■02 

247 

. OOOOOOE+OO 

-. 818687E- 

■03 

248 

. OOOOOOE+OO 

.  129613E- 

■02 

249 

. OOOOOOE+OO 

. 123003E- 

■02 

250 

. OOOOOOE+OO 

-. 864516E- 

-03 

251 

.  OOOOOOE+OO 

-. 149078E- 

■02 

252 

. OOOCOOE+00 

. 355459E- 

-03 

253 

. OOOOOOE+OO 

. 157839E- 

-02 

254 

. OOOOOOE+OO 

. 172373E- 

-03 

255 

.OOOOOOE+OO 

-.  149011E- 

-02 

256 

. OOOOOOE+OO 

-.660574E- 

-03 

257 

. OOOOOOE+OO 

. 124251E- 

-02 

258 

. OOOOOOE+OO 

. 105737E- 

-02 

259 

. OOOOOOE+OO 

-. 868928E- 

-03 

260 

. OOOOOOE+OO 

-.  132301E' 

-02 

261 

. OOOOOOE+OO 

.415134E- 

-03 

262 

.OOOOOOE+OO 

.  143365E' 

-02 

263 

. OOOOOOE+OO 

.  661320E' 

-04 

264 

. OOOOOOE+OO 

-. 138330E' 

-02 

265 

. OOOOOOE+OO 

-.521209E' 

-03 

266 

. OOOOOOE+OO 

.  118381E' 

-02 

267 

. OOOOOOE+OO 

.901413E- 

-03 

268 

. OOOOOOE+OO 

-.  862852E' 

-03 

269 

. OOOOOOE+OO 

-.  116816E 

-02 
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270 

. 000000E+00 

.460272E- 

•03 

271 

. 00O0O0E+O0 

. 129677E- 

■02 

272 

. 000000E+00 

-.  232701E- 

■04 

273 

. 000000E+00 

-. 127864E- 

■02 

274 

. 000000E+00 

-.399054E- 

•03 

275 

. 000000E+00 

. 112153E- 

■02 

276 

. 000000E+00 

.  761140E- 

■03 

277 

. O000O0E+OO 

-. 848092E- 

■03 

278 

.000000E+00 

-. 102580E- 

■02 

279 

. 000000E+00 

.492771E- 

■03 

280 

. 000000E+00 

. 116797E- 

•02 

281 

. 000000E+00 

-.976162E- 

■04 

282 

. 000000E+00 

-. 117693E- 

•02 

283 

. 000000E+00 

-. 292632E- 

■03 

284 

. 000000E+00 

. 105697E- 

■02 

285 

. 000000E+00 

.635534E- 

■03 

286 

. 000000E+00 

-.  826250E- 

■03 

287 

. 00000-E+00 

-. 895491E- 

■03 

288 

. ooooouz+oo 

.514358E- 

■03 

289 

. 000000E+00 

. 104737E- 

■02 

290 

. 000000E+00 

-. 158563E- 

-03 

291 

. OOnnQOL+00 

-. 107885E- 

■02 

292 

. 000000E+00 

-.  200536E- 

■03 

293 

. 000000E+00 

. 991213E- 

■03 

294 

. 000000E+00 

.523576E- 

■03 

295 

. 000000E+00 

-. 798744E- 

■03 

296 

. 000000E+00 

-. 776687E- 

-03 

297 

. 000000E+00 

.526597E- 

■03 

298 

. 000000E+00 

.934971E- 

•03 

299 

. 000000E+00 

-. 207643E- 

-03 

300 

. 000000E+00 

-. 984873E- 

-03 

301 

. 000000E+00 

-.  121428E- 

-03 

302 

. 000000E+00 

. 925211E- 

-03 

303 

. 000000E+00 

.424267E- 

-03 

304 

. 0O0OOOE+00 

-.  766821E- 

-03 

305 

. 000000E+00 

-.  668821E- 

-03 

306 

. 000000E+00 

.530897E- 

-03 

307 

. 000000E+00 

. 830670E- 

-03 

308 

.000000E+00 

-. 246269E' 

-03 

309 

. 000000E+00 

-.  895391E' 

-03 

310 

. 000OO0E+00 

-.540483E' 

-04 

311 

. 000000E+00 

.859741E' 

-03 

312 

. O000O0E+OO 

. 336627E' 

-03 

313 

. 000000E+00 

-.  731569E' 

-03 

314 

. 000000E+00 

-.571295E' 

-03 

315 

. 000000E+00 

.528523E' 

-03 

316 

. 000000E+00 

. 734302E' 

-03 

317 

. 000000E+00 

-.  275737E' 

-03 

318 

. 000000E+OO 

-.  810663E' 

-03 

319 

.000000E+00 

.278823E' 

-05 

320 

. 000000E+00 

. 795451E' 

-03 

321 

. 000000E+00 

.259710E' 

-03 

322 

. 000000E+00 

-.  693935E' 

-03 

323 

.  000O0OE+O0 

-.483492E' 

-03 

324 

.  000000E+00 

.520607E' 

-03 

325 

.  000000E+00 

.645634E' 

-03 
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326 

. 000000E+00 

-. 297233E-03 

327 

. 000000E+00 

-. 730852E-03 

328 

. 000OOOE+00 

.501878E-04 

329 

. 000000E+00 

. 732867E-03 

330 

. 000000E+00 

. 192606E-03 

331 

. 000000E+00 

-. 654736E-03 

332 

. 000000E+00 

-.404790E-03 

333 

. 000000E+00 

.508155E-03 

334 

.000000E+00 

.564390E-03 

335 

. 000000E+00 

-. 311833E-03 

336 

. 000000E+00 

-.656042E-03 

337 

. 000000E+00 

. 891797E-04 

338 

.000000E+00 

.672410E-03 

339 

.  000000E+00 

. 134443E-03 

340 

. O0OO0OE+O0 

-. 614672E-03 

341 

. 000000E+00 

-. 334566E-03 

342 

.  000000E+00 

.492056E-03 

343 

. 000000E+00 

.490253E-03 

344 

. OOOOOOE+OO 

-. 320515E-03 

345 

. 000000E+00 

-.586244E-03 

346 

. OOOOOOE+OO 

.  120718E-03 

347 

. OOOOOOE+OO 

.614406E-03 

348 

. OOOOOOE+OO 

.843958E-04 

349 

.  OOOOOOE+OO 

-.574335E-03 

350 

. OOOOOOE+OO 

-. 272207E-03 

351 

. OOOOOOE+OO 

.473097E-03 

352 

. OOOOOOE+OO 

.422879E-03 

353 

. OOOOOOE+OO 

-. 324162E-03 

354 

. OOOOOOE+OO 

-.521414E-03 

355 

. OOOOOCE+00 

. 145681E-03 

356 

. OOOOOOE+OO 

.559102E-03 

357 

.  OOOOOOE+OO 

.416829E-04 

358 

. 0nnc00E+00 

-.534224E-03 

359 

. OOOOOOE+OO 

-. 217110E-03 

360 

. OOOOOOE+OO 

.451961E-03 

361 

. OOOOOOE+OO 

. 361905E-03 

362 

. OOOOOOE+OO 

-. 323564E-03 

363 

.  OOOOOOE+OO 

-.461457E-03 

364 

. OOOOOOE+OO 

. 164877E-03 

365 

. OOOOOOE+OO 

.506671E-03 

366 

. OOOOOOE+OO 

.557033E-05 

367 

. OOOOOOE+OO 

-.494751E-03 

368 

.OOOOOOE+OO 

-. 168693E-03 

369 

. OOOOOOE+OO 

.429248E-03 

370 

. OOOOOOE+OO 

. 306958E-03 

371 

. OOOOOOE+OO 

-.319432E-03 

372 

. OOOOOOE+OO 

-.406239E-03 

373 

.  OOOOOOE+OO 

.  179045E-03 

374 

.OOOOOOE+OO 

.457227E-03 

375 

.  OOOOOOE+OO 

-.246287E-04 

376 

. OOOOOOE+OO 

-.456254E-03 

377 

. OOOOOOE+OO 

-. 126393E-03 

378 

.OOOOOOE+OO 

.405474E-03 

379 

. OOOOOOE+OO 

. 257656E-03 

380 

. OOOOOOE+OO 

-. 312396E-03 

381 

. OOOOOOE+OO 

-.355598E-03 
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382 

. 000000E+00 

. 188857E-03 

383 

. 0O0O00E+O0 

.410831E-03 

384 

. 000000E+00 

-.495535E-04 

385 

. 000000E+00 

-.419005E-03 

386 

. OO0O0OE+00 

-. 896748E-04 

387 

. 000000E+00 

. 381080E-03 

388 

. 000000E+00 

. 213620E-03 

389 

. 000000E+00 

-. 303017E-03 

390 

. 000000E+00 

-. 309343E-03 

391 

. 000000E+00 

. 194925E-03 

392 

. OOOOOOE+00 

. 367499E-03 

393 

.000000E+00 

-. 697972E-04 

394 

.OOOOOOE+00 

-.383214E-03 

395 

. 000000E+00 

-.580256E-04 

396 

. 000000E+00 

. 356444E-03 

397 

. 000000E+00 

. 174472E-03 

398 

. 000000E+00 

-. 291787E-03 

399 

. 000000E+00 

-. 267270E-03 

400 

. 000000E+00 

.  197800E-03 

1RTTTM  A     1    

C.     INPUT,  OUTPUT  DATA  OF  THE  AUGMENTED  IIR  NOTCH  FILTER:  P=  1 


INPUT  DATA  FOR  PROBLEM  #  1 


PROBLEM  #  1   INI'UT  DATA  SOURCEFILE:  dif feq.  tst 

THE  NUMBER  OF   INPUT  DELAYS:  L  =   3 

THE  NUMBER  OF  OUTPUT  DELAYS:  N  =   3 

THE  VALUE  OF  nsfcop  IS:  400 

THE  COEFFICIENTS  b(0),  b(l),  ...,  b(L)  ARE: 

990000E+00  -. 329900E-02   . 881146E+00   . 326631E+00 


THE  COEFFICIENTS  a(l),  ...,  a(N)  ARE: 
000000E+00  -.871246E+00  -. 323364E+00 

OUTPUT  DATA  FOR  PROBLEM  #  1 


116 


ns 

x(ns) 

y(ns) 

-3 

. 000000E+00 

.  OOOOOOE+00 

-2 

. 000000E+00 

.000000E+00 

-1 

. 000000E+00 

. OOOOOOE+00 

0 

. 100000E+01 

. 990000E+00 

1 

. 000000E+00 

-. 329900E-02 

2 

. 000000E+00 

.  186127E-01 

3 

. 00O0OOE+0O 

. 937460E-02 

4 

. 000000E+00 

-.  151494E-01 

5 

. 000000E+00 

-. 141862E-01 

6 

.000000E+00 

. 101675E-01 

7 

. 000000E+00 

. 172585E-01 

8 

. O000OOE+OO 

-.427105E-02 

9 

.000000E+00 

-. 183242E-01 

10 

. 000000E+00 

-. 185964E-02 

11 

.O0O000E+0O 

. 173460E-01 

12 

. 000000E+00 

. 754559E-02 

13 

. 00O000E+0O 

-. 145113E-01 

14 

. 000000E+00 

-. 121831E-01 

15 

. 000000E+00 

. 102029E-01 

16 

. 000000E+00 

. 153069E-01 

17 

. 000000E+00 

-.494966E-02 

18 

. 000000E+00 

-. 166353E-01 

19 

. 000000E+00 

-. 637331E-03 

20 

. 000000E+00 

. 160940E-01 

21 

. 000000E+00 

. 593455E-02 

22 

. 00C000E+00 

-.  138158E-01 

23 

. 000000E+00 

-. 103747E-01 

24 

. 000000E+00 

. 101179E-01 

25 

. 000O0OE+0O 

. 135064E-01 

26 

. 000000E+00 

-.546039E-02 

27 

. OOOOOOE+00 

-.  150392E-01 

28 

. 000000E+00 

. 389857E-03 

29 

. 000000E+00 

. 148685E-01 

30 

. 000000E+00 

.452347E-02 

31 

. 000000E+00 

-. 130802E-01 

32 

. 000000E+00 

-. 874900E-02 

JO 

. OOOOOOE+00 

. 993335E-02 

34 

. OOOOOOE+00 

. 118522E-01 

35 

. OOOOOOE+00 

-.582528E-02 

36 

. OOOOOOE+00 

-. 135383E-01 

"5  7 

. OOOOOOE+00 

. 124268E-02 

38 

. OOOOOOE+00 

. 136788E-01 

39 

. OOOOOOE+00 

.329511E-02 

40 

. OOOOOOE+00 

-. 123195E-01 

41 

. OOOOOOE+00 

-. 729410E-02 

42 

. OOOOOOE+00 

. 966777E-02 

43 

. OOOOOOE+00 

. 103386E-01 

44 

. 00000^E+00 

-. 606436E-02 

45 

. OOOOOCE+00 

-. 121337E-01 

46 

. OOOOOOE+00 

. 194041E-02 

47 

. OOOOOOE+00 

. 125324E-01 

48 

. ooonooE+oo 

. 223302E-02 

49 

. OOOOOOE+00 

-. 115463E-01 

50 

. 000000E+00 

-.599805E-02 
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51 

. 0000OOE+0O 

.933758E-02 

52 

. 000000E+00 

. 895943E-02 

53 

. 000000E+00 

-. 619577E-02 

54 

.  000000E+00 

-.  108253E-01 

55 

.000000E+00 

. 250088E-02 

56 

.  000000E+00 

. 114350E-01 

57 

. 000000E+00 

. 132163E-02 

58 

. 000000E+00 

-. 107714E-01 

59 

. 000000E+00 

-.484913E-02 

60 

.000000E+00 

. 895716E-02 

61 

.  000000E+00 

. 770786E-02 

62 

. 000000E+00 

-.623586E-02 

63 

. 000000E+00 

-.961187E-02 

64 

.000000E+00 

. 294052E-02 

65 

. 000000E+00 

. 103908E-01 

66 

. 000000E+00 

.546215E-03 

67 

. O0OOO0E+00 

-. 100038E-01 

68 

. 000000E+00 

-. 383588E-02 

69 

. 000000E+00 

. 853911E-02 

70 

. 000000E+00 

. 657685E-02 

71 

. 000000E+00 

-. 619928E-02 

72 

.  000000E+00 

-. 849130E-02 

73 

. 000000E+00 

. 327438E-02 

74 

. 000000E+00 

. 940263E-02 

75 

. 0O000OE+00 

-. 107009E-03 

76 

.  000000E+00 

-. 925082E-02 

77 

. 000000E+00 

-.294724E-02 

78 

. 000000E+00 

. 809434E-02 

79 

.  OOOOOOE+00 

. 555915E-02 

80 

. 000000E+00 

-.609913E-02 

81 

.  O0000OE+0O 

-. 746081E-02 

82 

.  000000E+00 

. 351621E-02 

83 

. 000000E+00 

. 847244E-02 

84 

. 000000E+00 

-. 650931E-03 

85 

.  000000E+00 

-. 851860E-02 

86 

.  000000E+00 

-. 217256E-02 

87 

. 000000E+00 

. 763228E-02 

88 

. 000000E+00 

.464744E-02 

89 

. 000000E+00 

-.594707E-02 

90 

. 000C00E+00 

-.651707E-02 

91 

.  000000E+00 

. 367854E-02 

92 

.000000E+00 

. 760104E-02 

93 

.000000E+00 

-. 109753E-02 

94 

. 000000E+00 

-. 781188E-02 

95 

. 000000E+00 

-.  150168E-02 

96 

. 000000E+00 

. 716097E-02 

97 

.000000E+00 

. 383442E-02 

98 

. 000000E+00 

-.575338E-02 

99 

.000000E+00 

-.565632E-02 

100 

. 000000E+00 

. 377269E-02 

101 

. 000000E+00 

. 678848E-02 

102 

.O0OOO0E+00 

-. 145789E-02 

103 

.000000E+00 

-. 713439E-02 

104 

. 000000E+00 

-. 924967E-03 

105 

. 000000E+00 

. 668724E-02 

106 

. 000000E+00 

. 311288E-02 
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107 

. 000000E+00 

-.552713E-02 

108 

. 000000E+00 

-.487450E-02 

109 

. 000000E+00 

. 380890E-02 

110 

. 000000E+00 

.603416E-02 

111 

. 000000E+00 

-. 174225E-02 

112 

. 000000E+00 

-. 648890E-02 

113 

. 000000E+00 

-.433302E-03 

114 

. 000000E+00 

. 621681E-02 

115 

.000000E+00 

.247579E-02 

116 

. 000000E+00 

-. 527625E-02 

117 

. 000000E+00 

-.416731E-02 

118 

. 000000E+00 

.379633E-02 

119 

.000000E+00 

.533691E-02 

120 

. 000000E+00 

-. 195998E-02 

121 

. 000000E+00 

-.587736E-02 

122 

.000000E+00 

-. 181364E-04 

123 

. 000000E+00 

.575441E-02 

124 

. 000000E+00 

. 191633E-02 

125 

.  000000E+C0 

-. 500764E-02 

126 

. OOOOOOE+00 

-. 353036E-02 

127 

. OOOOOOE+00 

. 374322E-02 

128 

. OOOOOOE+00 

.469510E-02 

129 

. OOOOOOE+00 

-. 211967E-02 

130 

. OOOOOOE+00 

-.530101E-02 

131 

. OOOOOOE+00 

. 328528E-03 

132 

. OOOOOOE+00 

.530391E-02 

133 

. OOOOOOE+OO 

. 142793E-02 

134 

. OOOOOOE+00 

-.472725E-02 

135 

. OOOOOOE+00 

-. 295917E-02 

136 

. OOOOOOE+OO 

. 365685E-02 

137 

. OOOOOOE+00 

.410679E-02 

138 

.  OOOOOOE+OO 

-. 222913E-02 

139 

. OOOOOOE+00 

-.476052E-02 

140 

. OOOOOOE+00 

. 614133E-03 

141 

. OOOOOOE+00 

.486840E-02 

142 

. OOOOOOE+00 

. 100432E-02 

143 

. OOOOOOE+00 

-.444016E-02 

144 

. OOOOOOE+OO 

-. 244927E-02 

145 

. OOOOOOE+00 

. 354371E-02 

146 

. OOOOOOE+OO 

. 356971E-02 

147 

. OOOOOOE+00 

-. 229544E-02 

148 

. OOOOOOE+00 

-.425600E-02 

149 

. OOOOOOE+OO 

. 845577E-03 

150 

. OOOOOOE+00 

.445029E-02 

151 

. OOOOOOE+00 

. 639532E-03 

152 

. OOOOOOE+00 

-.415072E-02 

153 

.  OOOOOOE+00 

-. 199625E-02 

154 

.OOOOOOE+00 

. 340950E-02 

155 

.  OOOOOOE+00 

. 308142E-02 

156 

. OOOOOOE+00 

-.232500E-02 

157 

. OOOOOOE+00 

-. 378719E-02 

158 

. OOOOOOE+00 

. 102922E-02 

159 

. OOOOOOE+00 

.405139E-02 

160 

. OOOOOOE+00 

.327932E-03 

161 

. OOOOOOE+OO 

-. 386257E-02 

162 

. OOOOOOE+00 

-.  159578E-02 
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163 

. 000000E+00 

. 325921E- 

■02 

164 

. O0O000E+O0 

.263934E- 

■02 

165 

. 00On00E+00 

-.232355E- 

■02 

166 

. 000000E+00 

-.  335342E- 

■02 

167 

. 000000E+00 

. 117092E- 

■02 

168 

.000000E+00 

. 367301E- 

•02 

169 

. 000000E+00 

.642165E- 

■04 

170 

. 000000E+00 

-.  357873E- 

■02 

171 

. 000000E+00 

-. 124367E- 

■02 

172 

. 000000E+00 

.  309719E- 

■02 

173 

. 000000E+00 

. 224077E- 

■02 

174 

.000000E+00 

-.229626E- 

■02 

175 

.000000E+00 

-.295378E- 

■02 

176 

.000000E+00 

.  127602E- 

■02 

177 

.000000E+00 

.331600E- 

•02 

178 

.000000E+00 

-. 156579E- 

■03 

179 

.O0000OE+0O 

-.  330167E- 

■02 

180 

. 0O0000E+OO 

-.  935856E- 

■03 

181 

. 000000E+O0 

. 292720E- 

■02 

182 

. 000000E+00 

. 188300E- 

■02 

183 

. OOOOOOE+00 

-. 224769E- 

■02 

184 

. 000000E+00 

-.258711E- 

■02 

185 

. 000000E+00 

.  134939E- 

■02 

186 

. 000000E+00 

. 298083E- 

■02 

187 

. 000000E+00 

-.339076E- 

■03 

188 

. 000O00E+OO 

-. 303338E- 

-02 

189 

. 000000E+00 

-.668474E- 

-03 

190 

. 000000E+00 

. 275246E- 

-02 

191 

. 000000E+00 

.  156329E- 

-02 

192 

. 000000E+00 

-.  218191E- 

-02 

193 

. 000000E+00 

-. 225206E- 

-02 

194 

. 000000E+00 

.  139547E- 

-02 

195 

. 000000E+00 

.  266765E- 

-02 

196 

. 000000E+00 

-.487564E- 

-03 

197 

. 000000E+00 

-.  277542E- 

-02 

198 

. 000000E+00 

-.437834E- 

-03 

199 

. 000000E+00 

. 257574E- 

-02 

200 

. 000000E+00 

. 127893E- 

-02 

201 

. 000000E+00 

-.210252E- 

-02 

202 

. 000000E+00 

-.  194717E- 

-02 

203 

. 000000E+00 

.  141825E' 

-02 

204 

. OOOOOOE+00 

. 237634E- 

-02 

205 

. OOOOOOE+00 

-. 606003E- 

-03 

206 

. 000000E+00 

-. 252899E' 

-02 

207 

. 000000E+00 

-.240445E' 

-03 

208 

. 000000E+00 

.239933E- 

-02 

209 

. 000000E+00 

. 102727E- 

-02 

210 

. 000000E+O0 

-. 201266E- 

-02 

211 

. 000000E+00 

-. 167086E- 

-02 

212 

. 000000E+00 

. 142134E- 

-02 

213 

. 000000E+00 

.210655E- 

-02 

214 

. 000000E+00 

-.698037E- 

-03 

215 

. 000000E+00 

-.229493E- 

-02 

216 

. 000000E+00 

-. 730215E- 

-04 

217 

. 000000E+00 

. 222517E- 

-02 

218 

. 000000E+00 

. 805719E' 

-03 
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219 

. 000000E+00 

-.  191506E- 

•02 

220 

. 000000E+00 

-. 142152E- 

•02 

221 

. O0000OE+O0 

. 140795E- 

•02 

222 

. 000000E+00 

. 185775E- 

•02 

223 

. 00000CE+00 

-. 767001E- 

•03 

224 

. 000000E+00 

-. 207384E- 

•02 

225 

. 000000E+00 

.675152E- 

•04 

226 

. 000000E+00 

. 205485E- 

•02 

227 

.  000000E+00 

. 611783E- 

■03 

228 

. 00000CE+00 

-. 181211E- 

■02 

229 

.  ooooocz+oo 

-. 119748E- 

■02 

230 

. 000000E+00 

. 138096E- 

•02 

231 

. 0OO000E+0O 

. 162927E- 

■02 

232 

. onnoocE+oo 

-.815938E- 

■03 

233 

. 000000E+00 

-. 186605E- 

■02 

234 

. 000000E+00 

. 184036E- 

■03 

235 

. O0O0O0E+OO 

. 188963E- 

■02 

236 

. 000000E+00 

.443071E- 

■03 

237 

. 000000E+00 

-. 170584E- 

■02 

238 

. 0000G0E+00 

-. 997063E- 

■03 

239 

. 00OO00E+00 

. 134294E- 

•02 

240 

. 00000CE+00 

. 142030E- 

■02 

241 

. 000000E+00 

-.  847613E- 

■03 

242 

. 000000E+00 

-.  167168E- 

■02 

243 

.  OOOOOOE+OO 

. 279207E- 

■03 

244 

. 000000E+00 

. 173054E- 

■02 

245 

. OOOOOOE+OO 

. 297304E- 

■03 

246 

. OOOOOOE+OO 

-.  159801E- 

-02 

247 

. OOOOOOE+OO 

-.  818618E- 

-03 

248 

. OOOOOOE+OO 

. 129612E- 

-02 

249 

. OOOOOOE+OO 

. 122996E- 

-02 

250 

. OOOOOOE+OO 

-. 864528E- 

-03 

251 

. OOOOOOE+OO 

-.  149071E- 

-02 

252 

. OOOOOOE+OO 

. 355493E- 

-03 

253 

. OOOOOOE+OO 

. 157833E- 

-02 

254 

. OOOOOOE+OO 

. 172321E- 

-03 

255 

. OOOOOOE+OO 

-. 149007E- 

-02 

256 

. OOOCOOE+00 

-.  660510E- 

-03 

257 

. OOOOOOE+OO 

. 124250E- 

-02 

258 

. OOOOOOE+OO 

. 105730E- 

-02 

259 

. OOOOOOE+OO 

-.  868935E' 

-03 

260 

. OOOOOOE+OO 

-. 132295E- 

-02 

261 

. OOOOOOE+OO 

. 415163E- 

-03 

262 

. OOOOOOE+OO 

. 143360E' 

-02 

263 

. OOOOOOE+OO 

. 660854E' 

■04 

2b4 

. OOOOOOE+OO 

-. 138326E' 

-02 

265 

. OOOOOOE+OO 

-.521150E- 

-03 

266 

. OOOOOOE+OO 

.  118379E' 

-02 

267 

. OOOOOOE+OO 

.901347E' 

-03 

268 

. OOOOOOE+OO 

-. 862854E' 

-03 

269 

. OOOOOOE+OO 

-.  116809E' 

-02 

270 

.OOOOOOE+OO 

.460295E' 

-03 

271 

. OOOOOOE+OO 

.  129671E' 

-02 

272 

. OOOOCOE+00 

-.233112E 

-04 

273 

. OOOOOOE+OO 

-. 127860E' 

-02 

274 

. OOOOOOE+OO 

-.  399000E' 

-03 
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275 

. 000000E+00 

. 112151E- 

•02 

276 

. 000000E+00 

. 761079E- 

■03 

277 

. 000000E+00 

-. 848089E- 

•03 

278 

. OO0O0OE+0O 

-. 102574E- 

•02 

279 

. 0000O0E+O0 

.492789E- 

■03 

280 

. 000000E+00 

. 116792E- 

■02 

281 

. 000000E+0O 

-.976517E- 

■04 

282 

.  000000E+00 

-. 117689E- 

■02 

283 

.  OOOOOOE+00 

-.292583E- 

■03 

284 

.  000000E+00 

. 105694E- 

■02 

285 

.  000000E+00 

. 635477E- 

■03 

286 

.  000000E+00 

-.826244E- 

■03 

287 

.  000000E+00 

-.895433E- 

■03 

288 

.  000000E+00 

.514371E- 

■03 

289 

.  000000E+00 

. 104732E- 

■02 

290 

. 000000E+00 

-.  158593E- 

■03 

291 

.  000000E+00 

-. 107880E- 

■02 

292 

. 000000E+00 

-. 200492E- 

■03 

293 

. 000000E+C0 

. 991185E- 

■03 

294 

. OOOOOOE+00 

.523523E- 

■03 

295 

.  00000''E+00 

-.  798734E- 

-03 

296 

. OOOOOCE+OO 

-. 776631E- 

■03 

297 

. OOOOOOE+00 

.526606E- 

■03 

298 

. OOOOOOE+00 

. 934919E- 

-03 

299 

. OOnocoE+00 

-.207668E- 

-03 

300 

.  OOOOOOE+00 

-. 984830E- 

-03 

301 

.  OOOOOOE+00 

-. 121389E- 

-03 

302 

.  OOOOOOE+00 

. 925181E- 

-03 

303 

. OOOOOOE+OO 

.424218E- 

-03 

304 

. OOOOOOE+00 

-. 766808E' 

-03 

305 

. OOOOOOE+OO 

-. 668769E' 

-03 

306 

. OOOOOOE+00 

.530901E' 

-03 

307 

. OOOOOOE+00 

. 830620E' 

-03 

308 

. OOOOOOE+OO 

-.246290E' 

-03 

309 

. OOOOOOE+00 

-.  895349E' 

-03 

310 

. OOOOOOE+OO 

-.540135E' 

-04 

311 

. OOOOOOE+00 

. 859710E' 

-03 

312 

. OOOOOOE+00 

.336583E' 

-03 

313 

. OOOOOOE+00 

-.  731553E' 

-03 

314 

. OOOOOOE+OO 

-.571245E' 

-03 

315 

. OOOOOOE+00 

.528524E' 

-03 

316 

. OOOOOOE+00 

. 734253E' 

-03 

317 

. OOOOOOE+00 

-.  275754E' 

-03 

318 

. OOOOOOE+OO 

-. 810621E' 

-03 

319 

. OOOOOOE+00 

. 281864E' 

-05 

320 

. OOOOOOE+00 

. 795419E- 

-03 

321 

. OOOOOOE+OO 

.259670E' 

-03 

322 

. OOOOOOE+00 

-. 693917E- 

-03 

323 

. OOOOOOE+00 

-.483446E' 

-03 

324 

. OOOOOOE+00 

.520605E' 

-03 

325 

. OOOOOOE+00 

. 645588E' 

-03 

326 

.  OOOOOOE+00 

-.297246E' 

-03 

327 

.  OOOOOOE+00 

-.  730811E' 

-03 

328 

.  OOOOOOE+00 

.502140E' 

-04 

329 

. OOOOOOE+OO 

. 732835E' 

-03 

330 

. OOOOOOE+00 

. 192569E' 

-03 

122 


331 

. 000000E+00 

-. 654717E-03 

332 

. 000000E+00 

-.404747E-03 

333 

. 0C0000E+00 

. 508149E-03 

334 

. 000000E+00 

.564346E-03 

335 

. 000000E+00 

-. 311842E-03 

336 

. 000000E+00 

-.656002E-03 

337 

. 000000E+00 

. 892021E-04 

338 

. OOOOOOE+OO 

. 672377E-03 

339 

.  000000E+00 

. 134410E-03 

340 

. OOOOOOE+OO 

-. 614651E-03 

341 

.  OOOOOOE+OO 

-. 334527E-03 

342 

. OOOOOOE+OO 

.492049E-03 

343 

. OOOOOOE+OO 

.490211E-03 

344 

.  OOOOOOE+OO 

-. 320521E-03 

345 

. OOOOOOE+OO 

-.586206E-03 

346 

.  OOOOOOE+OO 

. 120736E-03 

347 

. OOOOOOE+OO 

.614374E-03 

348 

. OOOOOOE+OO 

. 843669E-04 

349 

. OOOOOOE+OO 

-. 574313E-03 

350 

. OOOOOOE+OO 

-.  272171E-03 

351 

. OOOOOOE+OO 

.473087E-03 

352 

. OOOOOOE+OO 

.422840E-03 

353 

. OOOOOOE+OO 

-. 324164E-03 

354 

. OOOOOOE+OO 

-.521377E-03 

355 

. OOOOOOE+OO 

. 145696E-03 

356 

. OOOOOOE+OO 

.559070E-03 

357 

. OOOOOOE+OO 

.416575E-04 

358 

. OOOOOOE+OO 

-. 534201E-03 

359 

. OOOOOOE+OO 

-. 217077E-03 

360 

. OOOOOOE+OO 

.451950E-03 

361 

. OOOOOOE+OO 

. 361869E-03 

362 

.  OOOOO'E+00 

-. 323564E-03 

363 

. OOOOOOE+OO 

-.461421E-03 

364 

. OOOOOOE+OO 

. 164889E-03 

365 

. OOOOOOE+OO 

. 506640E-03 

366 

. OOOCGOE+00 

. 554827E-05 

367 

. OOOOOOE+OO 

-.494727E-03 

368 

. OOOOOOE+OO 

-.  168663E-03 

369 

. OOOOOOE+OO 

.429235E-03 

370 

. OOOOOOE+OO 

. 306924E-03 

371 

. OOOOOOE+OO 

-. 319430E-03 

372 

. OOOOOOE+OO 

-. 406206E-03 

373 

. OOOOOOE+OO 

. 179054E-03 

374 

. OOOOOOE+OO 

.457197E-03 

375 

. OOOOOOE+OO 

-. 246476E-04 

576 

. OOOOOOE+OO 

-.456231E-03 

377 

. OOOOOOE+OO 

-. 126367E-03 

378 

. OOOOOOE+OO 

.405459E-03 

379 

.  OOOOOOE+OO 

. 257625E-03 

380 

. OOOOOOE+OO 

-. 312392E-03 

381 

. OOOOOOE+OO 

-. 355566E-03 

382 

. OOOOOOE+OO 

. 188864E-03 

383 

. OOOOOOE+OO 

.410802E-03 

384 

. OOOOOOE+OO 

-.495696E-04 

385 

. OOOOOOE+OO 

-.418981E-03 

386 

. OOOOOOE+OO 

-. 896512E-04 
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387 

. 000000E+00 

. 381065E-03 

388 

. 000000E+00 

. 213592E-03 

389 

. 000000E+00 

-.  303011E-03 

390 

. OOO00CE+0O 

-. 309313E-03 

391 

. 000000E+00 

. 194929E-03 

392 

. 000000E+00 

. 367471E-03 

393 

. 000000E+00 

-. 698106E-04 

394 

. 000000E+00 

-. 383191E-03 

395 

. 000000E+00 

-.580047E-04 

396 

.000000E+00 

. 356428E-03 

397 

. 000000E+00 

. 174446E-03 

398 

. 000000E+00 

-.291779E-03 

399 

.000000E+00 

-.267242E-03 

400 

. 000000E+00 

. 197802E-03 

D.     INPUT,  OUTPUT  DATA  OF  THE  AUGMENTED  IIR  NOTCH  FILTER:  P  =  2 


INPUT  DATA  FOR  PROBLEM  #  1 


PROBLEM  #  1    INPUT  DATA  SOURCEFILE:  dif feq. tst 

THE  NUMBER  OF   INPUT  DELAYS:  L  =   4 

THE  NUMBER  OF  OUTPUT  DELAYS:  N  =   4 

THE  VALUE  OF  nsrop  IS:  400 

THE  COEFFICIENTS  b(0),  b(l),  ...,  b(L)  ARE: 

990000E+00  -. 329900E-02   . 186124E-01   . 614081E+00  -. 862533E+00 


THE  COEFFICIENTS  a(l),  ...,  a(N)  ARE: 

O00000E+00   . 0n0C0JE+00  -. 610814E+00   . 853908E+00 


OUTPUT  DATA  FOR  PROBLEM  #  1 


ns 

x(ns) 

y(ns) 

-4 

. 000000E+00 

. 00OO00E+00 

-3 

. 000000E+00 

. 000000E+00 

-2 

. 000000E+00 

. 000000E+00 
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-1 

. 000000E+00 

. OOOOOOE+OO 

0 

. 100000E+01 

. 990000E+00 

1 

. 000000E+00 

-. 329900E-02 

. 000000E+00 

. 186124E-01 

3 

. 000000E+00 

. 937506E-02 

4 

. 000000E+00 

-. 151490E-01 

5 

. 000000E+00 

-. 141858E-01 

6 

. 000000E+00 

. 101669E-01 

7 

. 000000E+00 

. 172586E-01 

8 

. 000000E+00 

-.427097E-02 

9 

. 000000E+00 

-. 183234E-01 

10 

. 000000E+00 

-. 186026E-02 

11 

. 000000E+00 

. 173461E-01 

12 

. 000000E+00 

. 754517E-02 

13 

. 000000E+00 

-. 145102E-01 

14 

. 000000E+00 

-. 121837E-01 

15 

. 000000E+00 

. 102032E-01 

16 

. 000000E+00 

. 153059E-01 

17 

. 000000E+00 

-.  494841E-02 

18 

. 000000E+00 

-. 166360E-01 

19 

. 000000E+00 

-. 636442E-03 

20 

. 000C00E+00 

.  160924E-01 

21 

.000000E+00 

.593604E-02 

22 

. 0000OOE+OO 

-. 138169E-01 

23 

. 000000E+00 

-. 103729E-01 

24 

. 000000E+00 

. 101156E-01 

25 

. 000000E+00 

.  135084E-01 

26 

. 000000E+00 

-.546244E-02 

27 

. 000000E+00 

-.  150363E-01 

28 

. 000000E+00 

. 386692E-03 

29 

. 000000E+00 

. 148715E-01 

30 

. 00O0OOE+0O 

.451996E-02 

31 

. 000000E+00 

-. 130758E-01 

32 

. OOOOOOE+OO 

-. 875350E-02 

33 

. OOOOOOE+OO 

. 993800E-02 

34 

. OOOOOOE+OO 

. 118465E-01 

35 

. OOOOOOE+OO 

-.581878E-02 

36 

. OOOOOOE+OO 

-. 135450E-01 

37 

. OOOOOOE+OO 

. 125013E-02 

38 

. OOOOOOE+OO 

. 136700E-01 

39 

. OOOOOOE+OO 

. 330475E-02 

40 

. OOOOOOE+OO 

-. 123297E-01 

41 

. OOOOOOE+OO 

-. 728235E-02 

42 

. OOOOOOE+OO 

.965436E-02 

43 

. OOOOOOE+OO 

.  103531E-01 

44 

. OOOOOOE+OO 

-.608030E-02 

45 

. OOOOOOE+OO 

-. 121155E-01 

46 

. OOOOOOE+OO 

. 192010E-02 

47 

. OOOOOOE+OO 

. 125545E-01 

48 

. OOOOOOE+OO 

. 220828E-02 

49 

. OOOOOOE+OO 

-.  115183E-01 

50 

. OOOOOOE+OO 

-. 602890E-02 

51 

. OOOOOOE+OO 

. 937158E-02 

52 

. OOOOOCE+00 

. 892123E-02 

53 

. OOOOOOE+OO 

-. 615305E-02 

54 

. OOOOOOE+OO 

-. 108724E-01 
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55 

. OOOOOOE+00 

. 255325E-02 

56 

. 000000E+00 

. 113763E-01 

57 

. 000000E+00 

. 138688E-02 

58 

. 000000E+00 

-. 108436E-01 

59 

. 000000E+OO 

-.476854E-02 

60 

. 000000E+00 

.886717E-02 

61 

. 000000E+00 

. 780770E-02 

62 

. 000000E+00 

-. 634675E-02 

63 

. 000000E+00 

-.948809E-02 

64 

. OOOOOOE+00 

. 280269E-02 

65 

. 000000E+00 

. 105437E-01 

66 

.000000E+00 

.375915E-03 

67 

. 000000E+00 

-.981388E-02 

68 

. 000000E+00 

-.404702E-02 

69 

. OOOOOOE+00 

. 877377E-02 

70 

.OOOOOOE+00 

.631545E-02 

71 

. 000000E+00 

-.590817E-02 

72 

. 000000E+00 

-. 881493E-02 

73 

. 000000E+00 

. 363443E-02 

74 

. 000000E+00 

. 900161E-02 

75 

. 000000E+00 

. 339247E-03 

76 

. O0000OE+00 

-.974710E-02 

77 

. 000000E+00 

-.  239484E-02 

78 

. 000000E+00 

. 747933E-02 

79 

. 000000E+00 

. 624335E-02 

80 

. 000000E+00 

-. 686032E-02 

81 

. 000000E+00 

-.661345E-02 

82 

. 000000E+00 

. 257313E-02 

83 

. 000000E+00 

. 952163E-02 

84 

. 000000E+00 

-. 181849E-02 

85 

. 000000E+00 

-. 721898E-02 

86 

. 000000E+00 

-.  361872E-02 

87 

. 000000E+00 

. 924135E-02 

88 

. 000000E+00 

. 285663E-02 

89 

. 000000E+00 

-. 395398E-02 

90 

. 000000E+00 

-.873480E-02 

91 

. 000000E+00 

. 614640E-02 

92 

. 000000E+00 

. 485444E-02 

93 

. 000000E+00 

. 195900E-02 

94 

. 000000E+00 

-. 112130E-01 

95 

. 000000E+00 

. 228329E-02 

96 

. 000000E+00 

.294866E-02 

97 

. 000000E+00 

. 852188E-02 

98 

. 000000E+00 

-. 109696E-01 

99 

. 000000E+00 

. 148640E-03 

100 

. 000000E+00 

-.268740E-02 

101 

. 000000E+00 

. 139773E-01 

102 

. 000000E+00 

-.945779E-02 

103 

.000000E+00 

. 176843E-02 

104 

. OOOOOOE+00 

-. 108323E-01 

105 

. 000000E+00 

. 177122E-01 

106 

.  000000E+00. 

-.915626E-02 

107 

.000000E+00 

.812659E-02 

108 

. 000000E+00 

-.200687E-01 

109 

. OOOOOOE+00 

. 207174E-01 

110 

. 000000E+00 

-.  127824E-01 

126 


Ill 

. 000000E+00 

. 191976E-01 

112 

. 000000E+00 

-. 297913E-01 

113 

. 000000E+00 

.254984E-01 

114 

. 0OOOO0E+O0 

-. 226412E-01 

115 

. 000000E+00 

. 345899E-01 

116 

.000000E+00 

-.410138E-01 

117 

. 000000E+00 

. 356029E-01 

118 

. 000000E+00 

-.404615E-01 

119 

.  000000E+00 

. 545884E-01 

120 

. 000000E+00 

-.567688E-01 

121 

.000000E+00 

.551160E-01 

122 

.000000E+00 

-. 678938E-01 

123 

.000000E+00 

. 812886E-01 

124 

. 000000E+00 

-. 821409E-01 

125 

. 000000E+00 

. 885345E-01 

126 

. 000000E+00 

-. 107627E+00 

127 

. 000000E+00 

. 119586E+00 

128 

. 000000E+00 

-. 124219E+00 

129 

. 000000E+00 

. 141341E+00 

130 

. 000000E+00 

-. 164948E+00 

131 

. 000000E+00 

. 177990E+00 

132 

. 000000E+00 

-. 192404E+00 

133 

. 000000E+00 

. 221445E+00 

134 

. 000000E+00 

-. 249570E+00 

135 

.  000000E+00 

. 269510E+00 

136 

. 000000E+00 

-. 299557E+00 

137 

. 000000E+00 

. 341534E+00 

138 

. 000000E+00 

-. 377730E+00 

139 

. 000000E+00 

. 413111E+00 

140 

. 000000E+00 

-.464408E+00 

141 

. 000000E+00 

. 522361E+00 

142 

. 000000E+00 

-.574880E+00 

143 

. 000000E+00 

.636425E+00 

144 

. 000000E+00 

-.  715627E+00 

145 

. 000000E+00 

. 797194E+00 

146 

. 000000E+00 

-. 879632E+00 

147 

. 000000E+00 

. 980564E+00 

148 

. 000000E+00 

-. 109802E+01 

149 

. 000000E+00 

. 121802E+01 

150 

. 000000E+00 

-. 135007E+01 

151 

. 000000E+00 

. 150800E+01 

152 

. 000000E+00 

-. 168159E+01 

153 

. 000000E+00 

. 186472E+01 

154 

. 000000E+00 

-. 207394E+01 

155 

. 000000E+00 

. 231483E+01 

156 

.  000000E+00 

-. 257492E+01 

157 

. 000000E+00 

. 285909E+01 

158 

. OO00OOE+00 

-. 318488E+01 

159 

. 000000E+00 

. 354945E+01 

160 

. 000000E+00 

-.394512E+01 

161 

. 000000E+00 

.438677E+01 

162 

. 000000E+00 

-.488765E+01 

163 

. 000000E+00 

.544063E+01 

164 

. 000000E+00 

-.604826E+01 

165 

. 000000E+00 

. 673134E+01 

166 

. 000000E+00 

-. 749681E+01 
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167 

. 000000E+00 

. 834016E+01 

168 

.  00000'  E+00 

-.927626E+01 

169 

. OOOOOCZ+00 

. 103271E+02 

170 

. 0O0000E+00 

-.  114959E+02 

171 

.  000000E+00 

. 127878E+02 

172 

. OO0000E+OO 

-. 142290E+02 

173 

. 000000E+00 

. 158402E+02 

174 

. 000000E+00 

-.  176274E+02 

175 

. 000000E+00 

. 196109E+02 

176 

. 000000E+00 

-.218257E+02 

177 

. 000OO0E+OO 

. 242932E+02 

178 

. 000O00E+0O 

-. 270308E+02 

179 

.000000E+00 

.300773E+02 

180 

. O000OOE+00 

-. 334757E+02 

181 

. 000000E+00 

. 372549E+02 

182 

. 000000E+00 

-.414534E+02 

183 

. 000000E+00 

.461307E+02 

184 

. 000000E+00 

-. 513410E+02 

185 

. 0C0000E+00 

.571326E+02 

186 

. 000000E+00 

-. 635747E+02 

187 

. 000000E+00 

. 707512E+02 

188 

. 000000E+00 

-. 787379E+02 

189 

. 000000E+00 

. 876183E+02 

190 

. 000000E+00 

-. 975028E+02 

191 

. 000000E+00 

. 108509E+03 

192 

. 000000E+00 

-. 120753E+03 

193 

. 000000E+00 

. 134374E+03 

194 

. 000000E+00 

-. 149537E+03 

195 

. 0000OOE+O0 

. 166415E+03 

196 

. 0OO00OE+0O 

-. 185190E+03 

197 

. 000000E+00 

.  206083E+03 

198 

. 000000E+00 

-. 229340E+03 

199 

. 000000E+00 

. 255219E+03 

200 

. 000000E+00 

-. 284013E+03 

201 

. 000000E+00 

. 316059E+03 

202 

.000000E+00 

-. 351727E+03 

203 

. 000000E+00 

. 391413E+03 

204 

. 000000E+00 

-.435575E+03 

205 

. 000000E+00 

.484725E+03 

206 

. 000000E+00 

-.539423E+03 

207 

. 000000E+00 

. 600286E+03 

208 

. O0OOOOE+00 

-.668018E+03 

209 

. 000000E+00 

. 743398E+03 

210 

. 000000E+00 

-.827281E+03 

211 

. 000000E+00 

. 920624E+03 

212 

.000000E+00 

-. 102450E+04 

213 

. 000000E+00 

. 114011E+04 

214 

.OOOOOOE+00 

-. 126875E+04 

215 

. OOOOOOE+00 

. 141191E+04 

216 

. 000000E+00 

-.  157123E+04 

217 

. OOOOOOE+00 

. 174852E+04 

218 

.000000E+00 

-. 194581E+04 

219 

. 000000E+00 

. 216537E+04 

220 

. 000000E+00 

-. 240970E+04 

221 

. 000000E+00 

.268160E+04 

222 

. 000000E+00 

-. 298418E+04 
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223 

. 000000E+00 

. 332090E+04 

224 

. 000000E+00 

-. 369562E+04 

225 

. 000000E+00 

.411262E+04 

226 

. 000000E+00 

-.457667E+04 

227 

. 000000E+00 

. 509308E+04 

228 

. 000000E+00 

-.566777E+04 

229 

.O0O000E+0O 

. 630729E+04 

230 

. 000000E+00 

-. 701898E+04 

231 

.000000E+00 

. 781098E+04 

232 

. OOOOOOE+OO 

-. 869233E+04 

233 

. OOOOOOE+OO 

. 967314E+04 

234 

. OOOOOOE+OO 

-. 107646E+05 

235 

. OOOOOfE+00 

. 119793E+05 

236 

. OOOOOGE+OO 

-. 133309E+05 

237 

. OOOOOOE+OO 

. 148351E+05 

238 

. OOOOOOE+OO 

-. 165091E+05 

239 

. OOOOOOE+OO 

. 183719E+05 

240 

. OOOOOOE+OO 

-. 204449E+05 

241 

. OOOOOOE+OO 

. 227518E+05 

242 

.OOOOOOE+OO 

-.  253191E+05 

243 

. OOOOOOE+OO 

. 281760E+05 

244 

. OOOOOOE+OO 

-. 313552E+05 

245 

. OOOOOOE+OO 

. 348932E+05 

246 

. OOOOOOE+OO 

-. 388304E+05 

247 

. OOOOOOE+OO 

. 432119E+05 

248 

. OOOOOOE+OO 

-.480877E+05 

249 

. OOOOOOE+OO 

.535138E+05 

250 

. OOOOOOE+OO 

-.595520E+05 

251 

. OOOOOOE+OO 

. 662716E+05 

252 

. OOOOOOE+OO 

-. 737494E+05 

253 

. OOOOOOE+OO 

. 820710E+05 

254 

. OOOOOOE+OO 

-. 913316E+05 

255 

. OOOOOOE+OO 

. 101637E+06 

256 

.OOOOOOE+OO 

-. 113105E+06 

257 

. OOOOOOE+OO 

. 125868E+06 

258 

. OOOOOOE+OO 

-. 140070E+06 

259 

. OOOOOOE+OO 

. 155875E+06 

260 

. OOOOOOE+OO 

-. 173463E+06 

261 

. OOOOOOE+OO 

. 193036E+06 

262 

. OOOOOOE+OO 

-. 214818E+06 

263 

.  OOOOOOE+OO 

. 239057E+06 

264 

. OOOOOOE+OO 

-. 266031E+06 

265 

. OOOOOOE+OO 

. 296049E+06 

266 

.OOOOOOE+OO 

-. 329454E+06 

267 

. OOOOOOE+OO 

. 366628E+06 

268 

.OOOOOOE+OO 

-.407997E+06 

269 

.  OOOOOOE+OO 

.454033E+06 

270 

. OOOOOOE+OO 

-.505265E+06 

271 

. OOOOOOE+OO 

.562277E+06 

272 

. OOOOOOE+OO 

-. 625722E+06 

273 

. OOOOOOE+OO 

. 696325E+06 

274 

. OOOOOOE+OO 

-. 774896E+06 

275 

. OOOOOOE+OO 

. 862332E+06 

276 

.OOOOOOE+OO 

-. 959634E+06 

277 

. OOOOOOE+OO 

. 106792E+07 

278 

. OOOOOOE+OO 

-.  118841E+07 
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279 

. 000000E+00 

. 132251E+07 

280 

. 000000E+00 

-. 147174E+07 

281 

. 000000E+00 

. 163780E+07 

282 

. 000000E+00 

-. 182260E+07 

283 

. 000000E+00 

. 202826E+07 

284 

. 000000E+00 

-.225712E+07 

285 

. 000000E+00 

. 251180E+07 

286 

. 000000E+00 

-. 279523E+07 

287 

. 000000E+00 

. 311063E+07 

288 

. 000000E+00 

-. 346162E+07 

289 

. 000000E+00 

. 385221E+07 

290 

.000000E+00 

-.428688E+07 

291 

. 000000E+00 

.477059E+07 

292 

. 000000E+00 

-.530889E+07 

293 

. 00000OE+OO 

.590792E+07 

294 

. 000000E+00 

-.657455E+07 

295 

. 000000E+00 

. 731639E+07 

296 

. 000000E+00 

-. 814194E+07 

297 

. O0O000E+0O 

. 906065E+07 

298 

. 000000E+00 

-. 100830E+08 

299 

. 000000E+00 

. 112207E+08 

300 

. 000000E+00 

-.  124868E+08 

301 

. 000000E+00 

. 138958E+08 

302 

.  OOOOO'.fE+OO 

-. 154637E+08 

303 

. OOOOOOE+00 

. 172086E+08 

304 

. OOOOOOE+OO 

-. 191504E+08 

305 

. OOOOOOE+00 

. 213112E+08 

306 

. 0^0C00£+00 

-. 237159E+08 

307 

. OOOOOOE+OO 

. 263919E+08 

308 

. OOOOOOE+00 

-. 293698E+08 

309 

. OOOOOOE+00 

. 326838E+08 

310 

. OOOOOOE+OO 

-. 363717E+08 

311 

. OOOOOOE+00 

.404757E+08 

312 

. OOOOOOE+OO 

-.450429E+08 

313 

. OOOOOOE+00 

.501253E+08 

314 

.OOOOOOE+OO 

-.557812E+08 

315 

. OOOOOOE+00 

. 620754E+08 

316 

. OOOOOOE+OO 

-. 690797E+08 

317 

. OOOOOOE+00 

. 768744E+08 

318 

. OOOOOOE+OO 

-. 855485E+08 

319 

. OOOOOOE+00 

. 952015E+08 

320 

. OOOOOOE+00 

-. 105944E+09 

321 

. OOOOOOE+00 

. 117898E+09 

322 

. OOOOOOE+00 

-. 131201E+09 

323 

. OOOOOOE+00 

. 146005E+09 

324 

. OOOOOOE+00 

-. 162480E+09 

325 

. OOOOOOE+00 

. 180813E+09 

326 

. OOOOOOE+00 

-. 201216E+09 

327 

. OOOOOOE+00 

. 223920E+09 

328 

. OOOOOOE+00 

-. 249186E+09 

329 

. OOOOOOE+00 

. 277303E+09 

330 

. OOOOOOE+00 

-. 308593E+09 

331 

. OOOOOOE+OO 

. 343413E+09 

332 

. OOOOOOE+00 

-. 382163E+09 

333 

. OOOOOOE+OO 

.425284E+09 

334 

. OOOOOOE+00 

-.473272E+09 
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335 

. 0OO0O0E+OO 

.526674E+09 

336 

. 000000E+00 

-.586101E+09 

337 

. OOOOOOE+OO 

. 652235E+09 

338 

. OOOOOOE+OO 

-. 725830E+09 

339 

. OOOOOOE+OO 

. 807730E+09 

340 

. OOOOOOE+OO 

-.898871E+09 

341 

. OOOOOOE+OO 

. 100030E+10 

342 

. OOOOOOE+OO 

-. 111316E+10 

343 

. OOOOOOE+OO 

. 123877E+10 

344 

. OOOOOOE+OO 

-. 137855E+10 

345 

.OOOOOOE+OO 

. 153410E+10 

346 

. OOOOOOE+OO 

-. 170720E+10 

347 

. OOOOOOE+OO 

. 189983E+10 

348 

. OOOOOOE+OO 

-. 211420E+10 

349 

. OOOOOOE+OO 

. 235276E+10 

350 

. OOOOOOE+OO 

-.261823E+10 

351 

.  OOOOOOE+OO 

.291366E+10 

352 

. OOOOOOE+OO 

-. 324243E+10 

353 

.  OOOOOOE+OO 

.360829E+10 

354 

. OOOOOOE+OO 

-.401544E+10 

355 

. OOOOOOE+OO 

. 446852E+10 

356 

. OOOOOOE+OO 

-.497273E+10 

357 

. OOOOOOE+OO 

.553384E+10 

358 

. OOOOOOE+OO 

-. 615825E+10 

359 

. OOOOOOE+OO 

. 685312E+10 

360 

. OOOOOOE+OO 

-. 762640E+10 

361 

. OOOOOOE+OO 

. 848693E+10 

362 

. OOOOOOE+OO 

-. 944456E+10 

363 

. OOOOOOE+OO 

. 105102E+11 

364 

. OOOOOOE+OO 

-. 116962E+11 

365 

. OOOOOOE+OO 

. 130159E+11 

366 

. OOOOOOE+OO 

-. 144846E+11 

367 

. OOOOOOE+OO 

. 161190E+11 

368 

. OOOOOOE+OO 

-. 179378E+11 

369 

.  00000r,E+00 

. 199618E+11 

370 

. OOOOOOE+OO 

-. 222142E+11 

371 

. OOOOOOE+OO 

. 247208E+11 

372 

.  oooooorj-oo 

-. 275102E+11 

373 

. OCOuOOE+00 

. 306143E+11 

374 

. OOOOOOE+OO 

-. 340687E+11 

375 

. OOOOOOE+OO 

. 379129E+11 

376 

. OOOOOOE+OO 

-.421908E+11 

377 

. OOOOOOE+OO 

.469514E+11 

378 

. OOOOOOE+OO 

-. 522492E+11 

379 

. OOOOOOE+OO 

.581448E+11 

380 

. OOOOOOE+OO 

-. 647056E+11 

381 

. OOOOOOE+OO 

. 720067E+11 

382 

. OOOOOOE+OO 

-.801317E+11 

383 

. OOOOOOE+OO 

. 891734E+11 

384 

. OOOOOOE+OO 

-. 992354E+11 

385 

. OOOOOOE+OO 

. 110433E+12 

386 

.OOOOCOE+00 

-. 122893E+12 

387 

. OOOOOOE+OO 

. 136760E+12 

388 

.  OOOOOOE+OO 

-. 152192E+12 

389 

. OOOOOOE+OO 

. 169364E+12 

390 

.OOOOOOE+OO 

-. 188475E+12 
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391 

. 000000E+00 

. 209741E+12 

392 

. 000000E+00 

-. 233408E+12 

393 

. 000000E+00 

. 259745E+12 

394 

. 000000E+00 

-. 289053E+12 

395 

. 000000E+00 

. 321669E+12 

396 

. 000000E+00 

-.357964E+12 

397 

. 000000E+00 

. 398356E+12 

398 

. 000000E+00 

-.443305E+12 

399 

. 000000E+00 

.493325E+12 

400 

. 000000E+00 

-.548990E+12 

E.     INPUT,  OUTPUT  DATA  OF  THE  AUGMENTED  IIR  NOTCH  FILTER:  P  =  3 


INPUT  DATA  FOR  PROBLEM  #  1 


PROBLEM  #  1    INPUT  DATA  SOURCEFILE:  dif feq.  tst 

THE  NUMBER  OF   INPUT  DELAYS:  L  =   5 

THE  NUMBER  OF  OUTPUT  DELAYS:  N  =   5 

THE  VALUE  OF  nstop  IS:  400 

THE  COEFFICIENTS  b(0),  b(l),  ...,  b(L)  ARE: 

990000E+00  -.329900E-02   . 186124E-01   . 937440E-02  -. 661007E+00  -. 604706E+00 


THE  COEFFICIENTS  a(l),  ...,  a(N)  ARE: 

000000E+00   .  000000E+00   . 000000E+00   . 652382E+00   . 598659E+00 


OUTPUT  DATA  FOR  PROBLEM  #  1 


ns 

x(.ns) 

y(ns) 

-5 

. 000000E+00 

. 000000E+00 

-4 

. 000000E+00 

. 000000E+00 

-3 

. 000000E+00 

. 000000E+00 

-2 

. OOOOOOE+00 

. 000000E+00 

-1 

. 000000E+00 

. 000000E+00 

0 

. 100000E+01 

. 990000E+00 

1 

. 000000E+00 

-.329900E-02 
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2 

. 000000E+00 

.  186124E-01 

3 

. O00O0OE+0O 

.937440E-02 

4 

. 000000E+00 

-.  151485E-01 

5 

. 0OO00OE+00 

-.  141854E-01 

6 

. 00O00OE+00 

. 101674E-01 

7 

. 000000E+00 

. 172582E-01 

8 

. 000000E+00 

-.427053E-02 

9 

. 000000E+00 

-.  183231E-01 

10 

. 000000E+00 

-. 185917E-02 

11 

. 000000E+00 

. 173458E-01 

12 

. 000000E+00 

. 754575E-02 

13 

.000000E+00 

-.  145103E-01 

14 

. 000000E+00 

-.  121822E-01 

15 

. 000000E+00 

. 102031E-01 

16 

. 000000E+00 

. 153069E-01 

17 

. 000000E+00 

-.494890E-02 

18 

.000000E+00 

-.  166341E-01 

19 

. 000000E+00 

-. 636683E-03 

20 

. 00C000E+00 

. 160941E-01 

21 

. 000000E+00 

.593505E-02 

22 

. 000000E+00 

-. 138145E-01 

23 

.  O00O00E+00 

-.  103735E-01 

24 

. 000000E+00 

. 101184E-01 

25 

. OOOOOOE+00 

.  135068E-01 

26 

. 000000E+00 

-.545927E-02 

27 

. 000000E+00 

-. 150377E-01 

28 

. 00OOOOE+O0 

. 390820E-03 

29 

. 000000E+00 

. 148691E-01 

30 

. 000000E+00 

.452445E-02 

31 

. 000000E+00 

-. 130786E-01 

32 

.  000000E+00 

-. 874750E-02 

33 

. 000000E+00 

. 993428E-02 

34 

. 000000E+00 

. 118532E-01 

35 

. 000000E+00 

-.582363E-02 

36 

. 000000E+00 

-. 135363E-01 

37 

. 000000E+00 

. 124418E-02 

38 

. 000000E+00 

. 136801E-01 

39 

. 000000E+00 

. 329678E-02 

40 

. 000000E+00 
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